An automated modeling method for space assembly process of variable topology aerospace structures

Abstract

The embodiment of the application discloses an automatic modeling method of a variable topology aerospace structure space assembly process, and the application can update and solve the whole structure dynamics equation by executing multiple assembly modeling processes to establish a dynamic model, realize automatic expansion and update of the dynamics model in the space assembly process, and represent the dynamics evolution law of the variable topology structure in the assembly process, so that the space assembly process of the variable topology aerospace structure can be accurately simulated, the modeling efficiency and accuracy can be improved, a large amount of data and complex assembly logic can be effectively processed, strong theoretical support is provided for the on-orbit assembly of the variable topology aerospace structure, a theoretical basis is provided for on-orbit assembly dynamics evolution analysis, and the application has important engineering application value. The application is widely applied in the technical field of computer software.

Description

Technical Field

[0001] The embodiments of the present invention relate to the fields of computer software technology and aerospace dynamics, and in particular to an automated modeling method for the space assembly process of variable topology aerospace structures. Background Technology

[0002] With the continuous development of aerospace technology, variable topology aerospace structures are increasingly widely used in space exploration, astronomical observation, and other fields. Variable topology aerospace structures refer to spacecraft systems whose structural form or connection relationships can be actively or passively changed during on-orbit operation, such as space solar power stations and ultra-large artificial gravity spacecraft. These structures can reach dimensions of hundreds or even thousands of meters. Due to their enormous size and complex structure, their scale far exceeds the payload limit of a single rocket launch. Therefore, "modular design and manufacturing, and on-orbit assembly" has become one of the main methods for constructing variable topology aerospace structures. Space assembly refers to the technical process of combining multiple independent modules or components into a complete functional structure according to a predetermined plan in outer space, through the operation of astronauts or space robots. During the space assembly process of variable topology aerospace structures, accurate modeling is crucial for predicting structural performance, guiding assembly operations, and ensuring structural stability.

[0003] However, traditional aerospace structure modeling methods are mostly based on overall static models, which are difficult to adapt to the problems of dynamic addition and subtraction of modules, real-time changes in connection relationships and complex environmental disturbances during space assembly. They are also difficult to meet the modeling and analysis needs of complex space assembly processes for variable topology aerospace structures. Summary of the Invention

[0004] To address the technical challenges of current aerospace structure modeling methods in meeting the modeling requirements of complex space assembly processes for variable topology aerospace structures, the present invention aims to provide an automated modeling method for the space assembly process of variable topology aerospace structures.

[0005] On one hand, embodiments of the present invention include an automated modeling method for the space assembly process of variable topology aerospace structures, the automated modeling method for the space assembly process of variable topology aerospace structures comprising: Obtain the model database; the model database stores parameter information of multiple assembly modules of the variable topology aerospace structure; Obtain the assembly database; the assembly database stores the assembly sequence information of each assembly module, and the assembly sequence is obtained according to intelligent planning; Execute multiple rounds of assembly modeling process and obtain the solution results of the overall structural dynamic equations after each round of assembly modeling process; The assembly modeling process in any round includes the following steps: An overall structural dynamic equation is established, which includes orbital dynamic equation, attitude dynamic equation, and structural dynamic equation, and there is mutual coupling between orbit, attitude, and structure. In this round of assembly modeling process, the overall structural dynamic equation is obtained through initialization. Otherwise, the overall structural dynamic equation updated in the previous round of assembly modeling process is obtained as the overall structural dynamic equation for this round of assembly modeling process. Based on the assembly sequence information, the parameter information of the assembly module corresponding to the assembly modeling process in this round is obtained from the model database; The overall structural dynamic equations are updated based on the assembly sequence information and the parameter information. Solve the updated overall structural dynamic equations.

[0006] Furthermore, the acquisition of the model database includes: Obtain structural information of the variable topology aerospace structure to be modeled; To obtain information on the carrying capacity of space launch vehicles; Based on the carrying capacity information, the structural information is subjected to fractal processing to obtain parameter information of multiple assembly modules; The parameter information of each assembly module is stored in the model database.

[0007] Further, updating the overall structural dynamic equations based on the assembly sequence information and the parameter information includes: Based on the assembly sequence information, constraint equations are added to the overall structural dynamics equations; Update the mass parameters in the overall structural dynamics equations based on the parameter information; Based on the updated mass parameters, the interaction coupling calculations between the various assembly modules in the overall structure are performed, and the coupling parameters in the overall structure dynamic equations are updated based on the calculation results; wherein, the overall structure is the structure described by the overall structure dynamic equations; Obtain space environment condition information, and update the external force parameters in the overall structural dynamic equations based on the space environment condition information.

[0008] Furthermore, solving the updated global structural dynamics equations includes: The updated overall structural dynamics equations were solved using numerical calculation methods. Based on the solution results, the dynamic response of the overall structure is obtained.

[0009] Furthermore, the automated modeling method for the space assembly process of the variable topology aerospace structure also includes: The parameter information read from the model database is activated.

[0010] Furthermore, the automated modeling method for the space assembly process of the variable topology aerospace structure also includes: The parameter information that has not been read from the model database is deactivated.

[0011] Furthermore, the multi-round assembly modeling process includes: The total number of rounds in the assembly modeling process is determined based on the total number of assembly modules in the model database. Establish the correspondence between each assembly module and each round of the assembly modeling process; Based on the assembly sequence information, determine the execution order of each round of the assembly modeling process; The assembly modeling process is executed sequentially in each round according to the execution order.

[0012] On the other hand, embodiments of the present invention also include an automated modeling system for the space assembly process of variable topology aerospace structures, the automated modeling system for the space assembly process of variable topology aerospace structures comprising: The first module is used to acquire a model database; the model database stores parameter information of multiple assembly modules of a variable topology aerospace structure; The second module is used to obtain the assembly database; the assembly database stores the assembly sequence information of each assembly module; The third module is used to execute multiple rounds of assembly modeling process and obtain the solution results of the overall structural dynamic equations after each round of assembly modeling process. The assembly modeling process in any round includes the following steps: Establish the overall structural dynamics equation; wherein, when the assembly modeling process described in this round is the first round, the overall structural dynamics equation is obtained through initialization; otherwise, the overall structural dynamics equation updated in the previous round of assembly modeling process is obtained as the overall structural dynamics equation for the assembly modeling process described in this round. Based on the assembly sequence information, the parameter information of the assembly module corresponding to the assembly modeling process in this round is obtained from the model database; The overall structural dynamic equations are updated based on the assembly sequence information and the parameter information. Solve the updated overall structural dynamic equations.

[0013] On the other hand, embodiments of the present invention also include a computer device, including a memory and a processor, the memory for storing at least one program, and the processor for loading at least one program to execute the automated modeling method for the space assembly process of variable topology aerospace structures in the embodiments.

[0014] On the other hand, embodiments of the present invention also include a computer-readable storage medium storing a processor-executable program, which, when executed by a processor, is used to perform the automated modeling method for the space assembly process of variable topology aerospace structures in the embodiments.

[0015] The beneficial effects of the embodiments of the present invention are as follows: The automated modeling method for the space assembly process of variable topology aerospace structures in the embodiments, by automatically modeling the space assembly process of variable topology aerospace structures and executing multiple assembly modeling processes to update and solve the dynamic equations of the overall structure, can establish a dynamic model, realize the automated expansion and updating of the dynamic model during the space assembly process, characterize the dynamic evolution law of the variable topology structure during the assembly process, and thus accurately simulate the space assembly process of variable topology aerospace structures. It can improve modeling efficiency and accuracy, effectively handle a large amount of data and complex assembly logic, provide strong theoretical support for the on-orbit assembly of variable topology aerospace structures, provide a theoretical basis for the dynamic evolution analysis of on-orbit assembly, and has important engineering application value. Attached Figure Description

[0016] Figure 1 This is a schematic diagram of the configuration of the multi-rotational joint space solar power station (MR-SPS) in the embodiment; Figure 2 This is a schematic diagram of the assembly process of the multi-rotational joint space solar power station (MR-SPS) in the embodiment; Figure 3 This is a schematic diagram of the main support structure of the multi-rotational joint space solar power station (MR-SPS) in the embodiment; Figure 4 This is a schematic diagram illustrating the steps of an automated modeling method for the space assembly process of variable topology aerospace structures in an embodiment. Figure 5 This is a flowchart and principle diagram of the automated modeling method for the space assembly process of variable topology aerospace structures in the embodiment; Figure 6 This is a schematic diagram of the assembly modeling process for each round in the embodiment; Figure 7 This is a schematic diagram showing the connection relationship between the assembly modules in the overall structure of the embodiment; Figure 8 This is a schematic diagram illustrating the principle of constraint processing in the embodiment; Figure 9This is a schematic diagram of the displacement curves of the space truss structure during the 6th and 12th assembly stages when the simulation verification was performed in the embodiment; Figure 10 This is a schematic diagram of the position error variation curve during simulation verification in the embodiment; Figure 11 This is a schematic diagram of the speed error variation curve during simulation verification in the embodiment; Figure 12 This is a schematic diagram of the attitude angle change curve during simulation verification in the embodiment; Figure 13 This is a schematic diagram of the attitude angular velocity change curve during simulation verification in the embodiment; Figure 14 This is a schematic diagram of the displacement response curve of the first motor during simulation verification in the embodiment; Figure 15 This is a schematic diagram of the displacement response curve of the second motor during simulation verification in the embodiment; Figure 16 This is a schematic diagram of the displacement response curve of the third motor during simulation verification in the embodiment; Figure 17 This is a schematic diagram of the displacement response curve of the fourth motor during simulation verification in the embodiment; Figure 18 This is a schematic diagram illustrating the parallel execution of each round of assembly modeling in the embodiment. Detailed Implementation

[0017] Traditional aerospace structure modeling methods are mostly based on overall static models, which are difficult to adapt to the problems of dynamic addition and subtraction of modules, real-time changes in connection relationships and complex environmental disturbances during space assembly. They are also difficult to meet the modeling requirements of complex space assembly processes for variable topology aerospace structures. Therefore, there is an urgent need for an automated modeling method to improve modeling efficiency and accuracy and ensure the smooth progress of space assembly of variable topology aerospace structures.

[0018] This embodiment provides an automated modeling method for the space assembly process of variable topology aerospace structures. Specifically, a space solar power satellite (SPS) is used as an example of a variable topology aerospace structure to illustrate the automated modeling method for the space assembly process of variable topology aerospace structures.

[0019] Space solar power stations (SPS) are a typical variable-topology space structure and an important solution for future clean energy acquisition. Traditional space solar power station designs face several technical challenges, such as the high technical difficulty of ultra-high-power conductive rotating joints in non-concentrated centralized power supply schemes, the inability of non-concentrated distributed power supply schemes to demonstrate continuous power supply advantages, and the challenges of designing high-precision concentrating systems in concentrating schemes. To simplify the technical difficulties of ultra-high-power conductive rotating joints in non-concentrated space solar power stations, the China Academy of Space Technology proposed a concept scheme for multiple rotationjoint space solar power stations (MR-SPS). The core of the MR-SPS is modular design, mainly composed of a central body (central truss structure), microwave transmitting antenna, connecting trusses, and a battery array [including solar panels and motor modules], with the following configuration: Figure 1 As shown.

[0020] Specifically, Figure 1 The Multi-Rotating Space Solar Power Station (MR-SPS) shown has a "one-line" supporting truss structure, and its assembly process is as follows: Figure 2 As shown. (Refer to...) Figure 2 The space solar power station (SPS) has a total truss structure length of 351 m, with each truss assembly module measuring 19.5 m in length. The entire assembly process is carried out in 18 stages. Compared to the structure in the first assembly stage, the final truss structure is 18 times larger.

[0021] Figure 1 The main support structure of the multi-rotational joint space solar power station (MR-SPS) shown is as follows: Figure 3 As shown, the automated modeling method for the assembly process proposed in this paper is validated. In this embodiment, using... Figure 3 Taking the modeling and simulation of the assembly process of the main support structure of the multi-rotational-joint space solar power station (MR-SPS) as an example, this paper illustrates the automated modeling method for the space assembly process of variable topology aerospace structures. The finite element modeling and numerical simulation steps used in the modeling process can all be implemented using custom-written programs in MATLAB.

[0022] Reference Figure 4 The automated modeling method for the space assembly process of variable topology aerospace structures includes the following steps: S1. Obtain the model database; S2. Obtain the assembly database; S3. Execute multiple rounds of assembly modeling process and obtain the solution results of the overall structural dynamic equations after each round of assembly modeling process.

[0023] Each round of assembly modeling process includes the following steps: S301. Establish the overall structural dynamics equation; wherein, when the current assembly modeling process is the first round, the overall structural dynamics equation is obtained through initialization; otherwise, the overall structural dynamics equation updated in the previous round of assembly modeling process is obtained as the overall structural dynamics equation for the current round of assembly modeling process. S302. Based on the assembly sequence information, obtain the parameter information of the assembly module corresponding to this round of assembly modeling process from the model database; S303. Update the overall structural dynamic equations based on the assembly sequence information and parameter information; S304. Solve the updated global structural dynamics equations.

[0024] Figure 4 The process and principle of automated modeling methods for the space assembly process of variable topology aerospace structures are as follows: Figure 5 As shown.

[0025] Reference Figure 4 and Figure 5 In step S1, the structural information of the variable topology space structure to be modeled [i.e., multi-rotational joint space solar power station (MR-SPS)] is obtained, that is... Figure 2 and Figure 3 The structural information is shown; next, the carrying capacity information of the space launch vehicle is obtained, specifically including the envelope size and carrying capacity of the launch vehicle used to transport the Multi-Rotating Joint Space Solar Power Station (MR-SPS). The carrying capacity information indicates the size of the object that the space launch vehicle can carry at one time, forming a constraint on the MR-SPS. If the MR-SPS exceeds the range of the carrying capacity information, then the MR-SPS needs to be disassembled into multiple independent and assembleable modules, i.e., assembly modules.

[0026] In this embodiment, based on the constraints of the space launch vehicle's carrying capacity information, the main support structure of the multi-rotating joint space solar power station (MR-SPS) is subjected to fractal processing, thereby disassembling it into assembly modules such as the central body, connecting trusses, and motor modules. For Figure 3 The main support structure of the multi-rotational joint space solar power station (MR-SPS) shown has the same number of connecting trusses and motor modules obtained by fractal processing (e.g., N connecting trusses and N motor modules).

[0027] In this embodiment, for each assembly module, its mass parameters (such as the mass and center of mass of the motor module), material properties (such as the elastic modulus, Poisson's ratio, and density of the material used in the truss module), and geometric dimensional information are obtained through measurement and material performance testing. Assembly modules of the same type can have the same geometric and material parameter configurations. In step S1, the obtained parameter information of these assembly modules is stored in the model database. Some of the parameter information stored in the model database is shown in Table 1.

[0028] Table 1. Parameter information of assembly modules stored in the model database.

[0029] In step S2, assembly sequence information is generated based on the assembly order of each assembly module in the model database. For example, for Figure 3 The main support structure of the Multi-Rotating Joint Space Solar Power Station (MR-SPS) shown can be assembled using a symmetrical assembly sequence, where connecting trusses and motor modules are assembled on both sides of the power station at different times. Therefore, with N connecting trusses and N motor modules, the entire assembly process can be divided into N stages. In the first stage, a set of connecting trusses and motor modules is installed on the foundation of the central body to form the overall structure. Then, in the second stage, a new set of connecting trusses and motor modules is installed on the overall structure of the first stage to form a new overall structure, and so on, achieving an assembly process that cycles through "connecting trusses-motor modules." The assembly sequence information records the module numbers, positions, and connection relationships with other modules to be installed at each stage of the assembly process, saved in a time-series format, and stored in the assembly database.

[0030] The principle behind step S2 is that the geometric configuration of the space structure gradually changes as modules are assembled throughout the on-orbit assembly phase. A reasonable assembly sequence is planned to reduce structural vibration and attitude changes during assembly. Taking into account the on-orbit state of the space structure (position, attitude, etc.), inherent structural properties (such as fundamental frequency), the motion trajectory of the assembly robot / manipulator, and fuel consumption, an optimization algorithm is used to rationally plan the assembly sequence, clarifying the order of each module during space assembly. Each module is assigned a unique number, its initial position in the overall structure is determined, and the connection relationships between modules, including connection methods and positions, are recorded. The number, position, connection relationships, and assembly sequence of each module are determined, and this data is stored in the assembly database. Through the establishment of the model database and the assembly database, systematic management of information on each module of the variable topology space structure is achieved.

[0031] In this embodiment, refer to Figure 6Based on the assembly sequence information, it is determined that phase 1, phase 2, ... phases need to be executed. N Multiple assembly stages with sequential relationships can be executed during step S3. N The assembly modeling process is carried out in rounds, and each round of assembly modeling corresponds to an assembly stage, thereby simulating the assembly stage.

[0032] In this embodiment, one round of assembly modeling process (the first) is used. i Taking the wheel assembly modeling process as an example, the following steps are described: S301. Establish the overall structural dynamics equations; S302. Based on the assembly sequence information, obtain the parameter information of the assembly module corresponding to this round of assembly modeling process from the model database; S303. Update the overall structural dynamic equations based on the assembly sequence information and parameter information; S304. Solve the updated global structural dynamics equations.

[0033] In this embodiment, the overall structural dynamic equation is a dynamic equation describing the overall structure, where the overall structure refers to the current execution of the... i Wheel assembly modeling process (corresponding stages) i At that time, the first round of assembly modeling process (corresponding to stage 1), the second round of assembly modeling process (corresponding to stage 2), ... have already been executed. i -1 round of assembly modeling process (corresponding stage) i -1) Assembled structure. For example, when i When =1, the overall structure refers to Figure 3 The central body in; when i When =2, the overall structure refers to the whole formed by the central body plus the assembly modules (connecting truss and motor modules) assembled on the central body in the first round of assembly modeling process (corresponding to stage 1).

[0034] In this embodiment, for Figure 3 The overall structural dynamic equations of the main support structure of the multi-rotational joint space solar power station (MR-SPS) shown can be derived through the following process: A dynamic model of the overall structure was established using finite element analysis and the first-order Lagrange method. Considering the connection relationship of the truss modules, appropriate element types and meshing parameters were set to accurately simulate the mass distribution and stiffness characteristics of the structure. Factors such as microgravity in the space environment were added to the model as load conditions.

[0035] a. Calculate kinetic energy and potential energy First, calculate the kinetic and potential energy of the overall structure. The kinetic energy mainly includes the kinetic energy of the central body, the kinetic energy of the connecting trusses, the kinetic energy of the motor modules, and the strain energy of the connecting trusses.

[0036] ① Kinetic energy of the central rigid body arbitrary mass element on the central body The position vector is represented as (1) in Mass element in the body coordinate system Relative center of mass position vector, This is the coordinate transformation matrix from the inertial coordinate system to the body coordinate system.

[0037] Differentiating the position vector yields (2) in Let be the angular velocity vector of the body at the center of the body coordinate system relative to the inertial coordinate system.

[0038] The kinetic energy of the central rigid body is expressed as (3) in For the mass of the central body, Let be the rotational inertia matrix of the central body.

[0039] ② Connecting the kinetic and potential energy of the truss The connecting truss module is discretized using Euler-Bernoulli beam elements, defining the first... The first truss module The generalized coordinates of each node are (4) in, For nodes Displacement vectors along each direction of the assembly module coordinate system. For nodes The rotation vector of the section around the three coordinate axes.

[0040] No. The generalized coordinates of each truss module node are (5) No. Any mass element on a truss module The position vector is represented as (6) in, This represents the position vector of the mass element relative to the connection point when no deformation occurs. The position vector, which is the position vector relative to the undeformed state, can be obtained through calculation using finite element theory. Differentiating the position vector yields (7) No. The kinetic energy of each truss module is expressed as follows: (8) in For the mass of the connecting truss.

[0041] No. The potential energy of each truss module is expressed as follows: (9) in Here is the stiffness matrix of the truss module.

[0042] ③ Motor module kinetic energy Definition of the first The generalized coordinates of each motor module are (10) No. The position vector of the centroid of each motor module in the inertial coordinate system is: (11) Differentiating the position vector yields (12) The kinetic energy of the motor module can be expressed as (13) in For the quality of the motor module.

[0043] The kinetic energy of the overall structure is expressed as (14) The potential energy of the overall structure is expressed as (15) The Lagrangian function of the overall structure is expressed as: (16) b. Constraint Handling During assembly, the constraints between modules consist of two parts: geometric constraints and boundary condition constraints. The connection diagram is as follows: Figure 7 As shown.

[0044] Here, constraint equations are used to characterize the connection relationships between the modules. The constraint equations are expressed as follows: (17) c. Dynamic equations Define generalized coordinates as (18) The first kind of Lagrange equation is expressed as: (19) Substituting equations (16) and (17) into the equations, we can obtain... (20) (twenty one) (twenty two) (twenty three) in, For the overall structural quality, For the rotational inertia of the structure, For the mass matrix of the truss module, For the truss module damping matrix, , These are the net external force and external torque acting on the structure, respectively. , These represent the external forces acting on the motor module, the external forces acting on the truss module, and the resultant external force, respectively. , , , , These are the corresponding coupling coefficient matrices.

[0045] Therefore, equations (20)-(23) constitute the overall structural dynamics equations to be established in step S301, where equation (20) represents the orbital motion equation, equation (21) represents the overall structural attitude dynamics model, equation (22) represents the motor module motion equation, and equation (23) represents the truss module vibration equation. The overall structural dynamics equations composed of equations (20)-(23) establish a dynamic model that includes the overall mechanical properties of the structure, including orbital dynamics, attitude dynamics, and structural dynamics, reflecting the design requirements and structural characteristics of variable topology aerospace structures. This equation comprehensively considers the influence of the structure's mass distribution, stiffness characteristics, damping characteristics, and external space environment factors on the structural dynamics behavior, providing a basic framework for subsequent assembly process modeling and a foundation for dynamic modeling during the assembly process. The dynamics equations can be gradually updated as assembly progresses.

[0046] In this embodiment, when executing step S3, the parameter information that has been read from the model database is activated. That is, the overall structural dynamics equations established in step S301 are considered to take into account the parameter information that has been read. For example, when executing the first step... iDuring the assembly modeling process, the overall structural dynamics equations established in step S301 consider the first round of assembly modeling, the second round of assembly modeling, and so on. i -1. Parameter information already read during the assembly modeling process; however, parameter information not read in the model database is deactivated, meaning that the unread parameter information is not considered in the overall structural dynamics equations established in step S301 (equivalent to the unread parameter information being invalid in the overall structural dynamics equations). For example, in the execution of the... i During the wheel assembly modeling process, the first wheel is not considered. i +1 round of assembly modeling process, the first i +2 rounds of assembly modeling process... Part N The parameter information is read only during the wheel assembly modeling process.

[0047] By deactivating unread parameter information, the amount of data that needs to be processed can be reduced, which helps to speed up the execution of each round of assembly modeling and reduce the demand for computing resources.

[0048] The dynamic characteristics of variable topology aerospace structures are analyzed from the perspective of structural dynamics. Assuming that during assembly, other influencing factors such as the connection gaps between assembly modules are not considered, only the truss structure before and after assembly is considered, and that the truss assembly modules are rigidly locked to the assembled parts immediately after assembly, forming an integral structure. Using rod elements as the basic units of the truss structure, and considering only the influence of structural dynamics, the structural dynamic equations of the overall space truss are derived based on (23): (twenty four) in, , , Let represent the mass matrix, damping matrix, and stiffness matrix of the entire space truss structure, respectively, where the damping is Rayleigh damping; Represents the displacement vector of a node; and These represent the control and disturbance forces, respectively. and These represent the position matrices of the control force and the disturbance force, respectively.

[0049] according to Figure 2 The assembly process shown proceeds step by step, updating the "activated" structural dynamic parameters. N The structural dynamics equations of the space truss structure during each assembly stage: (25) This analysis examines the dynamic response of variable topology spacecraft structures from the perspective of coupled dynamics involving orbit, attitude, and structure. Therefore, for the currently executing [number missing] [section missing]...i In the wheel assembly modeling process, if i =1, meaning that if the first round of assembly modeling is currently being performed, then since only the central body assembly module is currently available, and no parameter information of any truss modules, motor modules, or other assembly modules has been read, the overall structural dynamics equation in step 301 only retains the mechanical properties and other parameter information of the basic central body for model calculation. The parameter information related to the truss modules, motor modules, and other assembly modules in the overall structural dynamics equation is set to invalid. At this time, the overall structural dynamics equations shown in equations (20)-(23) can be simplified to (26) (27) Equations (26) and (27) represent orbital motion and attitude motion, respectively, which can be simply represented as rigid body motion.

[0050] In step S302, based on the assembly sequence information, the first... i Wheel assembly modeling process corresponding stages i Obtain the stage from the model database i The parameter information of the assembly module to be installed.

[0051] In step S303, the parameter information obtained in step S302 is activated and used to update the overall structural dynamic equations. The process of updating the overall structural dynamic equations using the parameter information specifically includes four steps: constraint handling, parameter updating, coupling calculation, and external disturbance loading. a. Constraint handling: Assume the first... i Wheel assembly modeling process (corresponding stages) i The assembly modules to be assembled into the overall structure are rigidly connected to the existing assembly modules within the overall structure. Based on the connection relationships between the assembly modules represented by the assembly sequence information in the assembly database, appropriate constraints are accurately applied in the model to ensure their correct position and orientation during assembly. Specifically, this can be achieved as follows: Figure 8 As shown, constraints such as "nodal degree of freedom loading" and "constraint equation connection" are used to constrain different assembly modules. For example, for the connection of truss modules, the relative displacement and rotation between truss modules are restricted to ensure that the connection method of the modules in the model is consistent with the actual assembly situation. By performing the above operations, the dynamic equations of the corresponding modules are activated, and the first dynamic equation is introduced into the overall structural dynamic equation before the update. i Wheel assembly modeling process (corresponding stages) i The constraint equations of the assembly module to be assembled.

[0052] b. Parameter Update: As assembly progresses and new assembly modules are added, the overall structural quality parameters will change. , The parameter changes abruptly, therefore according to the first... i Wheel assembly modeling process (corresponding stages) i The parameter information of the assembly modules to be assembled is used to update the mass parameters in the overall structural dynamic equation. Similarly, the material properties and other information of the new assembly modules can be extracted from the model database and integrated into the overall structural dynamic equation. At the same time, the parameters such as the moment of inertia and stiffness matrix in the overall structural dynamic equation are updated according to the actual assembly position and connection state of the new assembly modules, so that the updated overall structural dynamic equation model can accurately reflect the mechanical characteristics of the structure under the current assembly state.

[0053] c. Coupling Calculation: Based on the updates of parameters a and c, the interaction coupling of the structure is calculated. Specifically, consider the first... i Wheel assembly modeling process (corresponding stages) i The impact of the interactions between the assembled modules and the orbit-attitude-structure interactions on the overall structure's dynamic behavior, and their influence on the overall structural dynamic equations, are discussed. , , , , , , , , The coupling parameter matrix is ​​updated.

[0054] d. External disturbance loading: Based on actual space environment conditions (such as gravity gradient, space radiation, etc.), apply corresponding external disturbance loads to the overall structural dynamic equations, such as assembly contact collision force, solar radiation pressure, etc. For example, for equation (20) (orbit equation) in the overall structural dynamic equations, the resultant external force... Increase the gravity of the assembly module; for equation (21) (attitude equation) in the overall structural dynamics equation, change the net external torque. The gravitational gradient torque experienced by the structure; for equations (22) (vibration equation of the motor module) and (23) (vibration equation of the truss module) in the overall structural dynamics equations, the change in contact collision force , .

[0055] By executing the above process ad, step S303 is completed, and the overall structural dynamic equations are updated. Next, step S304 is executed, and the updated overall structural dynamic equations are solved by numerical calculation, thereby obtaining the dynamic responses of the overall structure, such as displacement, velocity, and acceleration, under the current assembly state, and analyzing the stability and mechanical properties of the structure.

[0056] Reference Figure 1 or Figure 2 After completing the first i Wheel assembly modeling process (corresponding stages) i After steps S301-S304 in ), if i < N If this indicates that the entire assembly modeling process has not yet been completed, then proceed to the next step. i +1 round of assembly modeling process (corresponding stage) i +1) steps S301-S304, wherein the execution of step i During step S301 of the +1 round assembly modeling process, the first step can be obtained. i Step S303 in the wheel assembly modeling process: Updated overall structural dynamic equations. Repeat this process until all steps are completed. N The wheel assembly modeling process enables automated modeling during spatial assembly.

[0057] In summary, the automated modeling method for the space assembly process of variable topology aerospace structures in this embodiment has the following technical effects: 1. Improved modeling efficiency: Through automated processes, the repetitive and tedious dynamic modeling work at different assembly stages is avoided. It can quickly complete the modeling of the space assembly process of variable topology aerospace structures, greatly improving modeling efficiency and meeting the requirements for modeling speed in engineering practice.

[0058] 2. Accurate description of the assembly process: Based on detailed structural fractals and database establishment, this method can accurately record and process the information of each module of the variable topology aerospace structure. In addition, it performs precise constraint processing, parameter updates and coupling calculations during the "activation" of modules, and fully considers various factors in the structural assembly process as well as the influence of the external space environment, so that the model can more realistically reflect the actual situation of the space assembly process of the variable topology aerospace structure.

[0059] 3. Optimize assembly schemes: By automating the assembly process of variable topology aerospace structures, different assembly sequences and external environmental conditions can be simulated in the model. The dynamic response and performance of the structure under various conditions can be analyzed, providing strong theoretical support for optimizing assembly schemes and helping to improve the success rate of space assembly of variable topology aerospace structures and the stability of the structures.

[0060] 4. Reduce costs and risks: Accurate and efficient modeling can identify potential problems in the space assembly process of variable topology aerospace structures in advance, avoiding design defects and safety hazards during actual assembly, thereby reducing R&D costs and engineering risks, and improving the economic benefits and reliability of aerospace engineering.

[0061] Simulation verification Efficiency analysis of automated modeling during the assembly process was conducted. In the simulation, the structural dynamics modeling in the automated modeling method for the space assembly process of the variable topology aerospace structure in this embodiment was numerically simulated and compared with the dynamic model established using the traditional finite element method. Here, the traditional modeling method refers to remodeling the already assembled overall structure after each new assembly module is assembled. The specific process of model building in both methods starts from the rod element, including the establishment of the structural dynamics model of the assembled parts at each assembly stage, until the establishment of the overall structural model at the final assembly stage. In this study, the model building and updating were implemented using a self-written program in MATLAB.

[0062] The structural dynamics models of the space truss at each assembly stage were established using the automated modeling method for the variable topology aerospace structure space assembly process in this embodiment, and numerical simulations were performed. The response curves of the truss structure at assembly stages 6 and 12 are shown below. Figure 9 As shown. After the sixth truss substructure is assembled, the length of the space truss structure is 117 m. Assuming the driving force at the end node is 0.5 N, the corresponding structural response is as follows. Figure 9 As shown in part (a); after the 12th truss substructure is assembled, the length of the space truss structure is 234 m. Assuming the driving force at the end node is 0.05 N, the corresponding structural response is as follows: Figure 9 As shown in part (b) of the document.

[0063] The advantages of automated modeling will be discussed below in terms of efficiency improvement, model scalability, functional scalability, automation level, and modeling continuity. Table 2 shows that, compared with traditional modeling methods, the automated modeling method for the variable topology aerospace structure space assembly process in this embodiment reduces the modeling time from 564.084 s to 145.93 s when establishing the dynamic model of the truss structure with 18 assembly stages, improving efficiency by approximately 74.13%. It should be noted that the space truss structure involved in this simulation only covers 18 substructures. If the overall structure contains hundreds or thousands of substructures, the advantage of this method in assembly modeling efficiency will become even more prominent. This efficiency improvement is due to the scalability of the proposed method. Using this method, the overall structural model is pre-built, and then only the corresponding structural model needs to be "activated" at each assembly stage, simplifying the modeling process and avoiding repetitive modeling work. In contrast, traditional methods require remodeling the assembled part of the overall structure or modeling the new type of substructure to be assembled before integration, which increases workload and complexity.

[0064] Furthermore, the automated modeling method for the space assembly process of variable topology aerospace structures in this embodiment also demonstrates advantages in terms of functional scalability. By pre-considering the relationship between the model of the substructure to be assembled and the model of the assembled parts, it facilitates the subsequent expansion and update design of the structure controller during the on-orbit assembly phase. Throughout the assembly phase, the dynamic modeling is coherent, reducing manual intervention and improving the degree of automation. In contrast, traditional modeling methods require manual analysis of new or new types of assembled substructures multiple times when establishing the structural dynamic model for the next assembly phase, which affects the continuity of modeling. Therefore, the proposed method not only improves the efficiency of dynamic modeling during the space assembly phase but also enhances the scalability and functional extensibility of the model, while ensuring the continuity and automation of the modeling process.

[0065] Table 2 Comparative Analysis of Automated Modeling Methods

[0066] To address the problem of coupled dynamics modeling during the on-orbit assembly of variable topology aerospace structures, this paper presents an automated modeling method for the space assembly process of variable topology aerospace structures. Simulation verification of the dynamics and control during the space assembly process of the MR-SPS main support structure is conducted. The establishment and updating of the dynamic model in this study are implemented by a self-written program in MATLAB.

[0067] Assume the initial attitude angle of the main structure is The attitude angular velocity is The time interval between two adjacent assembly is And a sliding mode controller is used to achieve stable control. For example... Figures 10-17 The curves for position error, velocity error, attitude angle, attitude angular velocity, and displacement response of the four motor modules are displayed respectively.

[0068] Figure 10 and Figure 11 The position and velocity error curves shown reflect the change in the deviation between the actual position and the desired commanded position of the central body over time. The curves indicate that the assembly process meets the stringent requirements for precise positioning. However, slight fluctuations occur in the velocity error curve during each assembly step, due to track-structure coupling. Effective suppression of velocity error is crucial for ensuring smooth movement and reducing impact, especially in precision docking operations.

[0069] Figure 12 , Figure 13 The diagram shows the attitude angle and angular velocity changes of the central body during assembly. Initially exhibiting attitude deviations, active control quickly reduced these deviations to within ±0.1°, maintaining high-precision attitude stability throughout the assembly process. Stable attitude is a prerequisite for ensuring accurate positioning and avoiding collisions during assembly.

[0070] Figure 14 , Figure 15 , Figure 16 and Figure 17 The displacement response curves of the first through fourth motor modules during the assembly process are shown. As key moving components connecting the satellite body to the large solar array, the dynamic characteristics of the motors directly affect the overall satellite stiffness and pointing accuracy. The displacement response curves of all four motor modules show vibration characteristics synchronized with the assembly process. Comparing the responses of the four modules, some differences can be observed in vibration amplitude, frequency, and decay time. This difference likely stems primarily from the different positions of the motor modules and the varying coupling effects they experience. This results in the fourth module having a relatively large structural amplitude, but under the constraint of the connecting truss module, it quickly decays and tends to a new equilibrium position, without exhibiting sustained divergence or constant-amplitude oscillations.

[0071] Based on the proposed automated modeling method for the spatial assembly process, we conducted assembly dynamics modeling and control simulation verification of the MR-SPS main support structure, which helps to reveal the dynamic evolution law during the assembly process.

[0072] In this embodiment, after executing all N During the wheel assembly modeling process, in addition to, Figure 6 The diagram shows the assembly modeling process executed sequentially in a serial manner. Alternatively, it can be implemented using... Figure 18 The parallel format shown. (Refer to...) Figure 18 It can be all NThe assembly modeling process is divided into multiple process groups. The assembly modeling processes within the same process group have a sequential execution order. For example, process group 1 includes consecutive rounds of assembly modeling processes, such as round 1, round 2, and so on. L The assembly modeling process is performed in rounds, and the number of rounds in the assembly modeling process within different process groups decreases sequentially according to the execution order. For example... Figure 18 You can set it in the middle. L > M - L > L - N That is, the number of rounds of the assembly modeling process included in process group 1 ( L The highest number of assembly modeling rounds is found in process group 2. M - L The number of rounds in the assembly modeling process included in process group 3 is moderate. L - N (Minimum). After partitioning, with sufficient computing resources, process group 1, process group 2, and process group 3 can be executed in parallel. During the execution of the same process group, each round of assembly modeling is executed sequentially. For example, in process group 1, the first round of assembly modeling, the second round of assembly modeling, and so on are executed sequentially. L Wheel assembly modeling process.

[0073] Figure 18 In this context, process group 1 corresponds to stage 1, stage 2, ... stage 1. L Simulation modeling was performed; process group 2 is equivalent to stage... L +1, Phase 2... Phase M Simulation modeling was performed; process group 3 is equivalent to stage... M +1, Phase 2... Phase N Simulation modeling was performed, and parallel execution can improve the modeling speed.

[0074] After the parallel execution of process group 1, process group 2 and process group 3, the execution results of each process group 1, process group 2 and process group 3 can be serially executed multiple times to obtain the final result.

[0075] For example, the execution result of process group 1 (updating the obtained overall structural dynamic equation) can be used as the initial overall structural dynamic equation obtained in the first round of assembly modeling, and the execution result of process group 2 can be used as the assembly module that needs to be assembled into the overall structure in the first round of assembly modeling, and the first round of assembly modeling can be executed; then, the execution result of process group 3 can be used as the assembly module that needs to be assembled into the overall structure in the third round of assembly modeling, and the second round of assembly modeling can be executed.

[0076] because Figure 6 The serial execution process shown is closer to the actual assembly process of variable topology aerospace structures, while Figure 18 Parallel execution processes can cause deviations from serial execution processes. By setting the number of rounds of assembly modeling processes in different process groups to decrease sequentially in the direction of execution order, it is possible to simulate the process of the overall structure of the variable topology aerospace structure increasing sequentially relative to the newly added assembly modules in the actual assembly process while achieving parallel execution. This reduces the deviation from the serial execution process and achieves a balance between modeling speed and modeling accuracy.

[0077] An automated modeling system for the space assembly process of variable topology aerospace structures can be run to execute the automated modeling method for the space assembly process of variable topology aerospace structures. Specifically, step S1 is executed by the first module of the automated modeling system, step S2 by the second module, and step S3 by the third module, thereby achieving the same technical effect as the automated modeling method for the space assembly process of variable topology aerospace structures.

[0078] A computer program can be written to execute the automated modeling method for the space assembly process of the variable topology aerospace structure in this embodiment. The computer program can be written into a computer device or storage medium. When the computer program is read out and run, the automated modeling method for the space assembly process of the variable topology aerospace structure in this embodiment can be executed, thereby achieving the same technical effect as the automated modeling method for the space assembly process of the variable topology aerospace structure in this embodiment.

[0079] It should be noted that, unless otherwise specified, when a feature is referred to as "fixed" or "connected" to another feature, it can be directly fixed or connected to the other feature, or indirectly fixed or connected to the other feature. Furthermore, the descriptions of "upper," "lower," "left," and "right" used in this disclosure are only relative to the relative positional relationships of the components of this disclosure in the accompanying drawings. The singular forms "a," "an," and "the" used in this disclosure are also intended to include the plural forms, unless the context clearly indicates otherwise. Moreover, unless otherwise defined, all technical and scientific terms used in this embodiment have the same meaning as commonly understood by one of ordinary skill in the art. The terminology used in this embodiment specification is only for describing specific embodiments and is not intended to limit the embodiments of the invention. The term "and / or" as used in this embodiment includes any combination of one or more of the associated listed items.

[0080] It should be understood that although the terms first, second, third, etc., may be used to describe various elements in this disclosure, these elements should not be limited to these terms. These terms are only used to distinguish elements of the same type from each other. For example, a first element may also be referred to as a second element without departing from the scope of this disclosure, and similarly, a second element may also be referred to as a first element. The use of any and all instances or exemplary language (“e.g.,” “such as,” etc.) provided in this embodiment is intended only to better illustrate embodiments of the invention and, unless otherwise required, does not impose a limitation on the scope of embodiments of the invention.

[0081] It should be recognized that embodiments of the present invention can be implemented or carried out by computer hardware, a combination of hardware and software, or by computer instructions stored in a non-transitory computer-readable storage medium. The method can be implemented using standard programming techniques—including a non-transitory computer-readable storage medium configured with a computer program, wherein such a storage medium causes the computer to operate in a specific and predefined manner—according to the methods and drawings described in the specific embodiments. Each program can be implemented in a high-level procedural or object-oriented programming language to communicate with the computer system. However, if desired, the program can be implemented in assembly or machine language. In any case, the language can be a compiled or interpreted language. Furthermore, for this purpose, the program can run on a programmed application-specific integrated circuit (ASIC).

[0082] Furthermore, the procedures described in this embodiment can be performed in any suitable order unless otherwise indicated by this embodiment or otherwise obviously contradict the context. The procedures (or variations and / or combinations thereof) described in this embodiment can be executed under the control of one or more computer systems configured with executable instructions, and can be implemented by hardware or a combination thereof as code (e.g., executable instructions, one or more computer programs, or one or more applications) that commonly executes on one or more processors. A computer program includes a plurality of instructions executable by one or more processors.

[0083] Furthermore, the method can be implemented in any suitable type of computing platform, including but not limited to personal computers, minicomputers, mainframes, workstations, networked or distributed computing environments, standalone or integrated computer platforms, or in communication with charged particle tools or other imaging devices, etc. Aspects of embodiments of the invention can be implemented as machine-readable code stored on a non-transitory storage medium or device, whether removable or integrated into a computing platform, such as a hard disk, optical read and / or write storage medium, RAM, ROM, etc., such that it is readable by a programmable computer, and when the storage medium or device is read by the computer, it can be used to configure and operate the computer to perform the processes described herein. Furthermore, the machine-readable code, or portions thereof, can be transmitted via wired or wireless networks. The invention of this embodiment includes these and other different types of non-transitory computer-readable storage media when such media comprises instructions or programs that implement the steps above in conjunction with a microprocessor or other data processor. Embodiments of the invention also include the computer itself when programmed according to the methods and techniques of embodiments of the invention.

[0084] A computer program can be applied to input data to perform the functions of this embodiment, thereby transforming the input data to generate output data stored in non-volatile memory. The output information can also be applied to one or more output devices, such as a display. In a preferred embodiment of the invention, the transformed data represents physical and tangible objects, including a specific visual depiction of physical and tangible objects generated on the display.

[0085] The above are merely preferred embodiments of the present invention. The embodiments of the present invention are not limited to the above-described implementations. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the embodiments of the present invention, as long as they achieve the same technical effects, should be included within the scope of protection of the embodiments of the present invention. Within the scope of protection of the embodiments of the present invention, the technical solutions and / or implementation methods can have various modifications and variations.

Claims

1. A method for automated modeling of a variable topology aerospace structure space assembly process, characterized by, The automated modeling method for the space assembly process of the variable topology aerospace structure includes: Obtain structural information of the variable topology aerospace structure to be modeled; To obtain information on the carrying capacity of space launch vehicles; Based on the carrying capacity information, the structural information is subjected to fractal processing to obtain parameter information of multiple assembly modules; The parameter information of each assembly module is stored in the model database; the model database stores the parameter information of multiple assembly modules of the variable topology aerospace structure. Obtain the assembly database; the assembly database stores the assembly sequence information of each assembly module; The assembly modeling process is executed in multiple rounds, and the solution results of the overall structural dynamic equations updated in each round of the assembly modeling process are obtained, thereby obtaining the dynamic response of the overall structure under the current assembly state; the dynamic response includes displacement, velocity and acceleration. The assembly modeling process in any round includes the following steps: Establish the overall structural dynamics equations; wherein, when the assembly modeling process described in this round is the first round, the overall structural dynamics equations are obtained through initialization; otherwise, the overall structural dynamics equations updated in the previous round of assembly modeling process are obtained as the overall structural dynamics equations for the assembly modeling process described in this round. The overall structural dynamics equations include orbital dynamics equations, attitude dynamics equations, and structural dynamics equations, and there is mutual coupling between orbit, attitude, and structure. Based on the assembly sequence information, the parameter information of the assembly module corresponding to the current assembly modeling process is obtained from the model database; the parameter information read from the model database is activated, and the parameter information not read from the model database is deactivated; the deactivation process is used to invalidate the unread parameter information in the overall structural dynamic equation. Based on the assembly sequence information and the parameter information, the overall structural dynamic equation is updated so that the updated overall structural dynamic equation model reflects the mechanical properties of the structure in the current assembly state. Solve the updated overall structural dynamics equations; The overall structural dynamic equation is: in, For the overall structural quality, For the rotational inertia of the structure, For the mass matrix of the truss module, For the truss module damping matrix, The net external force acting on the structure. The external torque acting on the structure The external force acting on the motor module The external forces acting on the truss module For the net external force received, , , , , , , , The coupling coefficient matrix, With the center of mass, Let be the angular velocity vector of the central body relative to the inertial coordinate system in the body coordinate system. For the first Generalized coordinates of each motor module For the first The generalized coordinates of each truss module node.

2. The automated modeling method for the space assembly process of variable topology aerospace structures according to claim 1, characterized in that, The step of updating the overall structural dynamic equations based on the assembly sequence information and the parameter information includes: Based on the assembly sequence information, constraint equations are added to the overall structural dynamics equations; Update the mass parameters in the overall structural dynamics equations based on the parameter information; Based on the updated mass parameters, the interaction coupling calculations between the various assembly modules in the overall structure are performed, and the coupling parameters in the overall structure dynamic equations are updated based on the calculation results; wherein, the overall structure is the structure described by the overall structure dynamic equations; Obtain space environment condition information, and update the external force parameters in the overall structural dynamic equations based on the space environment condition information.

3. The automated modeling method for the space assembly process of variable topology aerospace structures according to claim 2, characterized in that, Solving the updated global structural dynamics equations includes: The updated overall structural dynamics equations were solved using numerical calculation methods. Based on the solution results, the dynamic response of the overall structure is obtained.

Images