A method and apparatus for impact damage analysis of an attitude control power system

Abstract

This invention discloses a method and apparatus for impact damage analysis of an attitude control power system, relating to the field of attitude control power system technology. It addresses the problems of over- or under-testing in existing impact damage analysis methods, as well as high computational costs and low efficiency. The method includes: constructing a whole-machine impact response model of the target attitude control power system; using an explicit solver to determine the impact response of the whole-machine impact response model, obtaining the impact response results; identifying the target weak component based on the stress response of each component in the whole-machine impact response model; constructing an impact damage model of the target weak component; using an explicit solver to determine the impact damage of the impact damage model, obtaining the impact damage data of the target weak component. The impact damage analysis method for an attitude control power system provided by this invention can avoid over- or under-testing in the impact damage analysis of attitude control power systems and reduce computational costs.

Description

Technical Field

[0001] This invention relates to the field of attitude control power system technology, and in particular to an impact damage analysis method and apparatus for attitude control power systems. Background Technology

[0002] An attitude control propulsion system is generally composed of components such as a frame, thrust chamber, gas generator, gas pipelines, liquid pipelines, storage tanks, and high-pressure gas cylinders, all organically connected together. During operation, this high-dimensional dynamic system may be subjected to various intense environmental loads, including vibrations, shocks, and pressure pulsations. Connections such as duct roots, lugs, support plates, and clamps are prone to damage under dynamic loads, thus threatening the safe operation of the attitude control engine.

[0003] Existing methods for impact damage analysis that use local critical structure models generally apply impact loads directly to the local critical structures. However, external load excitation is usually transmitted from the docking plane between the power system and the compartment to the power system structure. During the transmission process, due to the damping or resonance effects of various components, the load may attenuate or amplify. Therefore, it is easy to apply incorrect load boundaries, resulting in over- or under-testing. On the other hand, the method of establishing a whole engine model for damage analysis has huge computational costs and low efficiency, and is not suitable for actual engineering requirements. Summary of the Invention

[0004] The purpose of this invention is to provide a method and apparatus for impact damage analysis of attitude control power systems, which solves the problems of over- or under-testing of existing impact damage analysis methods, as well as high computational cost and low efficiency.

[0005] To achieve the above objectives, the present invention provides the following technical solution:

[0006] In a first aspect, the present invention provides a method for impact damage analysis of an attitude control power system, comprising:

[0007] A whole-machine impact response model of the target attitude control power system is constructed; the whole-machine impact response model includes finite element models of each component;

[0008] An explicit solver is used to determine the impact response of the whole machine impact response model, and the impact response results are obtained; the impact response results include the stress response and displacement response of each component;

[0009] The target weak component is determined based on the stress response of each component in the overall impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine.

[0010] Construct an impact damage model for the target weak component;

[0011] An explicit solver is used to determine the impact damage of the impact damage model, thereby obtaining the impact damage data of the target weak component; the boundary condition of the impact damage model is the displacement response corresponding to the impact damage model in the impact response results.

[0012] Optionally, the overall impact response model of the target attitude control power system includes:

[0013] A three-dimensional model of the target attitude control power system is constructed based on the geometric structure of the target attitude control power system; the target attitude control power system includes a thin-walled structure, a first force transmission structure, a second force transmission structure, and a heavy-duty structure;

[0014] Finite element models are constructed for the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure; the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure are connected and assembled into a whole machine finite element model by binding connection or multi-point constraint.

[0015] The material properties of the whole machine finite element model are set to a non-destructive nonlinear constitutive model to obtain the whole machine impact response model.

[0016] Optionally, constructing finite element models for the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure includes:

[0017] In the three-dimensional model, the first force transmission structure is equivalent to the central axis, the thin-walled structure is equivalent to the central plane, and the heavy structure is equivalent to the center of mass, thus obtaining a simplified three-dimensional model.

[0018] The simplified 3D model is meshed, the central axis is simulated using beam elements, the center surface is simulated using shell elements, the center of mass is simulated using mass points and assigned a corresponding mass, and the second force transmission structure is simulated using solid elements.

[0019] Optionally, constructing the impact damage model of the target weak component includes:

[0020] Based on Saint-Venant's principle, the target weak component and the model within the target connection area of ​​the target weak component are extracted from the overall impact response model to obtain the first model;

[0021] If the first model includes beam elements and shell elements, then import the original geometric model corresponding to the beam element at the position of the beam element, import the original geometric model corresponding to the shell element at the position of the shell element, and use solid elements to simulate to obtain the second model;

[0022] By changing the material properties in the second model to a strain rate-dependent nonlinear dynamic damage constitutive model, a third model is obtained.

[0023] The contact between the target weak component and the surrounding connecting components is established in the third model to obtain the fourth model;

[0024] The displacement response of the nodes or surfaces on the boundary of the fourth model in the impact response results is applied to the fourth model to simulate the boundary conditions, thereby obtaining the impact damage model of the target weak component.

[0025] Optionally, when the target weak component is a liquid-bearing structure, the contact between the target weak component and the surrounding connecting components is established in the third model to obtain the fourth model, which further includes:

[0026] A fluid domain geometry model is imported into the target weak component of the third model, and the fluid domain geometry model is meshed to obtain fluid domain elements;

[0027] The fluid domain element is set as an arbitrary Lagrange-Euler region, and the properties of the fluid domain element are set as Euler properties;

[0028] The outer wall constraint of the target weak component is set as a Lagrangian constraint, and the constraint at the free liquid surface of the target weak component is set as an Euler boundary.

[0029] The contact surface between the target weak component and the fluid domain unit is coupled using the coupled Eulerian-Lagrange method;

[0030] The material properties of the fluid domain element are set using state equations to obtain the third model.

[0031] Optionally, determining the target weak component based on the stress response of each component in the overall impact response model includes:

[0032] Components whose stress response exceeds the yield strength are identified as target weak components.

[0033] Optionally, the analysis steps for determining the impact response of the whole machine impact response model and determining the impact damage of the impact damage model are the same.

[0034] Optionally, the material constitutive model of the impact damage model is a nonlinear dynamic damage constitutive model related to strain rate.

[0035] Optionally, when determining impact damage in the impact damage model, if the damage variable of a unit in the impact damage model is 1, then the unit with a damage variable of 1 is deleted.

[0036] Compared with existing technologies, the present invention provides an impact damage analysis method for an attitude control power system, comprising: constructing a whole-machine impact response model of the target attitude control power system; the whole-machine impact response model includes finite element models of each component; using an explicit solver to determine the impact response of the whole-machine impact response model and obtain impact response results; the impact response results include the stress response and displacement response of each component; determining the target weak component based on the stress response of each component in the whole-machine impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under whole-machine impact; constructing an impact damage model of the target weak component; using an explicit solver to determine the impact damage of the impact damage model and obtain impact damage data of the target weak component. This application, by determining the impact response of the entire machine, can accurately and efficiently identify the target weak components that undergo plastic deformation under impact, i.e., the damaged components. Then, only the target weak components need to be analyzed for impact damage, without having to determine which components are susceptible to damage by performing impact damage analysis on each component of the entire machine. This greatly reduces computational costs and improves analysis efficiency. In addition, the boundary conditions of the impact damage model of the target weak components in this application are the displacement responses corresponding to the impact damage model in the impact response results obtained from determining the impact response of the entire machine. This can ensure the accuracy of the boundary conditions applied to the impact damage model, avoid over-testing or under-testing, and improve the accuracy of impact damage analysis.

[0037] Secondly, the present invention also provides an impact damage analysis device for an attitude control power system, comprising:

[0038] The whole machine shock response model construction module is used to construct the whole machine shock response model of the target attitude control power system; the whole machine shock response model includes the finite element model of each component;

[0039] The impact response determination module is used to determine the impact response of the whole machine impact response model using an explicit solver, and obtain the impact response results; the impact response results include the stress response and displacement response of each component;

[0040] The target weak component determination module is used to determine the target weak component based on the stress response of each component in the whole machine impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine;

[0041] The impact damage model construction module is used to construct the impact damage model of the target weak component;

[0042] The impact damage determination module is used to determine the impact damage of the impact damage model using an explicit solver to obtain the impact damage data of the target weak component; the boundary condition of the impact damage model is the displacement response corresponding to the impact damage model in the impact response results.

[0043] Compared with the prior art, the beneficial effects of the impact damage analysis device for attitude control power system provided by the present invention are the same as the beneficial effects of the impact damage analysis method for attitude control power system described in the above technical solution, and will not be repeated here. Attached Figure Description

[0044] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this invention, illustrate exemplary embodiments of the invention and are used to explain the invention, but do not constitute an undue limitation of the invention. In the drawings:

[0045] Figure 1 A flowchart of an impact damage analysis method for an attitude control power system provided by the present invention;

[0046] Figure 2 This is a schematic diagram of the impact damage analysis device for an attitude control power system provided by the present invention. Detailed Implementation

[0047] To facilitate a clear description of the technical solutions in the embodiments of the present invention, the terms "first" and "second" are used to distinguish identical or similar items with essentially the same function and effect. For example, the first threshold and the second threshold are merely used to distinguish different thresholds and do not limit their order. Those skilled in the art will understand that the terms "first" and "second" do not limit the quantity or execution order, and that the terms "first" and "second" are not necessarily different.

[0048] It should be noted that in this invention, the terms "exemplary" or "for example" are used to indicate examples, illustrations, or descriptions. Any embodiment or design described as "exemplary" or "for example" in this invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of terms such as "exemplary" or "for example" is intended to present the relevant concepts in a concrete manner.

[0049] In this invention, "at least one" refers to one or more, and "more than one" refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone, where A and B can be singular or plural. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "At least one of the following" or similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one of a, b, or c can represent: a, b, c, a combination of a and b, a combination of a and c, a combination of b and c, or a, b, and c, where a, b, and c can be single or multiple.

[0050] Before introducing the embodiments of the present invention, the relevant terms involved in the embodiments of the present invention are first defined as follows:

[0051] ALE, short for Arbitrary Lagrangian-Eulerian, is a numerical computation method that combines the advantages of the Lagrangian and Eulerian descriptions. It is primarily used to handle complex problems involving large deformations, moving boundaries, or fluid-structure interactions. Its core idea is to construct a mesh framework independent of material motion, allowing for local reconstruction of the mesh in time and space to adapt to changes in boundary conditions while maintaining the quality of the computational mesh.

[0052] CLE, short for Coupled Eulerian-Lagrangian, is a computational method in finite element method software used to simulate multiphase flow and solid-fluid interactions. It combines the advantages of the Eulerian and Lagrangian methods, enabling more accurate simulation of complex engineering problems.

[0053] ABAQUS: A finite element analysis software widely used in engineering simulation and design verification. Its full name is "Advanced Simulation for Engineering and Sciences." It focuses on solving problems ranging from simple linear analysis to complex nonlinear problems, using the finite element method to numerically simulate various engineering problems, helping engineers and researchers optimize designs and improve product performance and safety.

[0054] SPH particles: In smooth particle hydrodynamics, particles are the basic units used in meshless numerical computation methods to simulate the behavior of continuous media.

[0055] During the design phase of attitude control engines, when evaluating key local structures, simulation calculations often result in over- or under-evaluation due to inaccurate application of load boundaries, necessitating further verification of the engine's dynamic strength through experiments. Furthermore, traditional engineering calculation models typically use static strength verification theories, such as tensile strength and elongation at fracture, to verify impact strength, neglecting the strain rate effect and damage caused by impact loads. This leads to significant errors between the calculated results and actual conditions.

[0056] To address the aforementioned problems, this invention provides a method and apparatus for impact damage analysis of an attitude control power system. The method uses a first-level whole-machine impact response model to determine the overall impact response, identifies the target weak component, obtains the load boundary of the target weak component, and uses a second-level impact damage model to obtain the structural damage under transient impact loads. The following description is in conjunction with the accompanying drawings.

[0057] See Figure 1 The impact damage analysis method for an attitude control dynamic system provided by the present invention includes the following steps:

[0058] Step 100: Construct the overall impact response model of the target attitude control power system; the overall impact response model includes the finite element models of each component.

[0059] The overall impact response model is a finite element model, including beam elements, shell elements, mass points, and second force transmission structure solid elements. The overall impact response model can be constructed using finite element software, including combining shape, material properties, boundary conditions, etc.

[0060] Specifically, step 100 can be implemented based on the following steps:

[0061] Step 110: Construct a three-dimensional model of the target attitude control power system based on its geometric structure; the target attitude control power system includes a thin-walled structure, a first force transmission structure, a second force transmission structure, and a heavy-duty structure;

[0062] For example, thin-walled structures include gas cylinders, hoops, etc.; first force transmission structures include conduits, bolts, etc., which transmit force in a single direction; second force transmission structures include joints, supports, etc., which are key force transmission structures; and heavy-duty structures include thrust devices, valves, etc., which are large-mass non-force transmission structures.

[0063] Step 120: Construct finite element models of the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure; the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure are connected and assembled into a whole machine finite element model by binding connection or multi-point constraint.

[0064] Specifically, step 120 may include the following steps:

[0065] Step 121: In the three-dimensional model, the first force transmission structure is equivalent to the central axis, the thin-walled structure is equivalent to the central plane, and the heavy structure is equivalent to the center of mass, to obtain a simplified three-dimensional model;

[0066] In practical applications, first delete non-critical parts such as chamfers and holes; extract and retain the central axis of conduits, bolts, etc., and delete the geometric models of conduits and bolts; extract the center surface of thin-walled structures such as storage tanks, gas cylinders, and hoop, delete the geometric model and retain the center surface; extract the center of mass of large-mass non-force-transmitting structures such as thrust devices and valves, delete the geometric model of large-mass non-force-transmitting structures and retain the center of mass; retain the key force-transmitting structures.

[0067] Step 122: Mesh the simplified 3D model, simulate the central axis using beam elements, simulate the center surface using shell elements, simulate the center of mass using mass points and assign corresponding mass, and simulate the second force transmission structure using solid elements.

[0068] In practical applications, a relatively new mesh is used, ensuring the mesh count does not exceed 300,000. After modeling, the model quality is checked, including aspects such as element aspect ratio, warping, mesh distortion, coincident nodes, normal direction, element size, model boundaries, and penetration. The established whole-machine impact response model must ensure that the simplified model is consistent with the original structure's mass characteristics and stiffness, and that the load transfer relationship is consistent, in order to accurately and efficiently obtain the whole-machine impact response, identify target weak components, and accurately obtain the load boundaries of the target weak components.

[0069] Step 130: Set the material properties of the whole machine finite element model to a non-destructive nonlinear constitutive model to obtain the whole machine impact response model.

[0070] For liquid-filled structures such as liquid-filled tanks and liquid-carrying conduits, the liquid mass needs to be added. The simplified, high-mass non-load-bearing structures are assigned a mass point of equal mass at their center of mass. The model mass is checked to ensure that the simplified model mass matches the actual structure.

[0071] Step 200: Use an explicit solver to determine the impact response of the whole machine impact response model and obtain the impact response results; the impact response results include the stress response and displacement response of each component.

[0072] For example, in ABAQUS software, first set up the analysis step. In the analysis step settings, select explicit dynamic analysis, set the initial time step and maximum time step, and set the total time length of the analysis step. Then apply boundary conditions to the model: apply load boundaries to the model according to the actual working conditions. Output settings: output stress response, acceleration response, and displacement response. Select explicit solver type for solving. Finally, submit the analysis job in the job module of the software and wait for the solver to complete the calculation. Use the post-processor to view the simulation results, i.e., the impact response results.

[0073] Step 300: Determine the target weak component based on the stress response of each component in the overall impact response model.

[0074] The weakest component of the target is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine;

[0075] Specifically, the stress response of each component is compared with the yield stress of the component material. Components whose stress response exceeds the yield strength of the component material are identified as target weak components. When the stress response exceeds the yield stress of the material, the material begins to undergo plastic deformation. Plastic deformation is an irreversible permanent deformation; once plastic deformation occurs, the material cannot recover its original shape and mechanical properties even if the external force is removed. Therefore, components that undergo plastic deformation can be considered target weak components.

[0076] Step 400: Construct an impact damage model for the target weak component; the material constitutive model of the impact damage model is a strain rate-dependent nonlinear dynamic damage constitutive model.

[0077] Specifically, step 400 can be implemented based on the following steps:

[0078] Step 410: Based on Saint-Venant's principle, extract the target weak component and the model within the target connection area of ​​the target weak component from the overall impact response model to obtain the first model; the geometric boundary of the first model corresponds to the overall impact response model.

[0079] The core of Saint-Venant's principle is that the stress distribution caused by a load applied to an elastic body, far from the load application area, is essentially related only to the resultant force and resultant moment of the load, and not to the specific distribution of the load. This means that the specific distribution of the load only affects the stress distribution near the load application area, while far away from this area, the stress distribution is uniform.

[0080] In practical applications, the target connection area of ​​the target weak component is first randomly or determined based on experience. It is then determined whether the stress distribution in this area is non-uniform. If so, the target connection area is expanded; if not, the target connection area is reduced until the stress distribution in the area outside the target connection area is uniform.

[0081] Step 420: If the first model includes beam elements and shell elements, then import the original geometric model corresponding to the beam element onto the beam element, import the original geometric model corresponding to the shell element at the position of the shell element, and use solid elements for simulation, that is, generate a solid network to obtain the second model;

[0082] Step 430: Change the material properties in the second model to a strain rate-dependent nonlinear dynamic damage constitutive model to obtain the third model;

[0083] For example, the Jonson-Cook damage constitutive model can be selected as the nonlinear dynamic damage constitutive model. Choosing the nonlinear dynamic damage constitutive model can solve the problem in traditional engineering calculation models, which usually use static strength verification theories, such as tensile strength and elongation after fracture, to verify impact strength, without considering the strain rate effect and damage caused by impact load, resulting in a large error between the calculation results and the actual situation.

[0084] Step 440: Establish the contact between the target weak component and the surrounding connecting components in the third model to obtain the fourth model;

[0085] In practical applications, if the target weak component is in contact with the surrounding connected structure, a contact pair needs to be established.

[0086] Step 450: Apply the displacement response of the nodes or surfaces on the boundary of the fourth model in the impact response results to the fourth model to simulate the boundary conditions, and obtain the impact damage model of the target weak component.

[0087] In practical applications, the displacement response of nodes or surfaces on the boundary of the fourth model is selected, the file path of the result database file of the whole machine impact response model is specified, the corresponding analysis step and frame are selected, and the displacement or acceleration response obtained from the impact response calculation is selected as the load boundary of the impact damage model at the key weak position, so as to ensure that the load boundary of the target weak component is accurate during the load transfer process in the first-level and multi-level models.

[0088] Step 500: Use an explicit solver to determine the impact damage of the impact damage model and obtain the impact damage data of the target weak component; the boundary condition of the impact damage model is the displacement response corresponding to the impact damage model in the impact response results.

[0089] For example, in ABAQUS software, the output is set to output stress response, acceleration response, displacement response, and damage variable; element deletion is defined when determining impact damage in an impact damage model, deleting elements with a damage variable of 1 when the damage variable in the impact damage model is 1; the analysis step is set to be the same for determining impact response in the whole machine impact response model and determining impact damage in the impact damage model. If the whole machine impact response model is a multi-step analysis, multiple analysis steps are defined simultaneously in the impact damage model. Submit the analysis job in the job module of the software and wait for the solver to complete the calculation. Use the post-processor to view the simulation results, i.e., the impact damage data.

[0090] After modeling is completed, it is necessary to check whether the simulation results of the impact damage model and the simulation results of the whole machine impact response model are reasonably connected at the boundary. Compare the acceleration and displacement responses obtained from the impact damage model with those obtained from the whole machine impact response model to see if they correspond to each other, and analyze the damage of the target weak component.

[0091] As can be seen from the above method, this application can accurately and efficiently identify the target weak components that undergo plastic deformation under the impact of the whole machine by determining the impact response of the whole machine, that is, the components that are damaged. Then, only the target weak components need to be analyzed for impact damage, without having to determine the impact damage of each component of the whole machine to determine which components are susceptible to damage, which greatly reduces the computational cost and improves the analysis efficiency. In addition, the boundary conditions of the impact damage model of the target weak components in this application are the displacement response corresponding to the impact damage model in the impact response results obtained by determining the impact response of the whole machine. This can ensure the accuracy of the boundary conditions applied to the impact damage model, avoid over-testing or under-testing, and improve the accuracy of impact damage analysis.

[0092] The established impact damage model for the target weak component must ensure accurate calculation of dynamic response and damage mode under high strain rate impact loads. This is especially important for special liquid-bearing containers, such as liquid-filled tanks and liquid-bearing conduits, where the influence of liquid mass, stiffness, and self-sloshing on the structural impact response and damage under impact loads must be considered. Traditional acoustic-structure interaction methods are computationally efficient but cannot calculate nonlinear large deformation problems. While the SPH particle method can accurately simulate liquids within containers, its computational cost is too high and its stability is low, making it unsuitable for impact damage analysis. Therefore, a suitable method must be selected to simulate liquids within the structure. Specifically, when the target weak component is a liquid-bearing structure, the following steps are included before step 440:

[0093] Step 431: Import the fluid domain geometry model into the target weak component of the third model, and mesh the fluid domain geometry model to obtain fluid domain elements;

[0094] In this case, the mesh of the fluid domain element needs to be smaller than that of the solid element.

[0095] Step 432: Set the fluid domain unit to an arbitrary Lagrange-Euler region, and set the properties of the fluid domain unit to Euler properties;

[0096] For example, in the ABAQUS software mesh, an adaptive mesh domain is selected, the fluid domain element is specified as an arbitrary Lagrange-Euler region, and Euler properties are assigned.

[0097] Step 433: Set the outer wall constraint of the target weak component to a Lagrangian constraint, and set the constraint at the free liquid surface of the target weak component to an Eulerian boundary;

[0098] For example, the outer walls of the storage tank, liquid-filled conduit, etc. are set as Lagrangian constraints, and the free liquid surface is set as an Eulerian boundary.

[0099] Step 434: Couple the contact surface between the target weak component and the fluid domain unit using the coupled Eulerian-Lagrange method;

[0100] The Coupled Eulerian-Lagrange method (CEL) is a numerical simulation technique that combines Eulerian and Lagrange descriptions. It is primarily used to handle engineering problems involving complex physical phenomena such as large deformations, fluid-structure interaction, or multiphase flow. In practical applications, Eulerian meshes are used to describe fluid motion, while Lagrange meshes are used to describe solid deformation. The interaction between the fluid and solid is achieved through contact relationships, including pressure transfer and momentum exchange.

[0101] Step 435: Set the material properties of the fluid domain unit using the equation of state to obtain the third model.

[0102] Specifically, the EOS equation of state is used to define the liquid.

[0103] Step 436: Define fluid output: output fluid pressure; initial conditions of the fluid domain: set the hydrostatic pressure of the fluid. After the impact damage model is built, the pressure contour map of the fluid domain needs to be checked to verify whether the internal energy and kinetic energy are reasonable and to verify the accuracy of the fluid-structure interaction model.

[0104] This invention employs the ALE method to establish the fluid domain within the impact damage model. By dynamically adjusting the liquid mesh to adapt to boundary motion, it solves the problems of mesh distortion and numerical convergence caused by large deformations. It eliminates the need for mesh reconstruction, significantly improving computational efficiency while maintaining computational accuracy and stability. Furthermore, by combining nonlinear dynamic damage constitutive models and damage models, fluid-structure interaction analysis calculations under impact loads are performed with greater accuracy.

[0105] The embodiments of the present invention can divide functional modules according to the above method examples. For example, each function can be divided into its own functional module, or two or more functions can be integrated into one processing module. The integrated module can be implemented in hardware or as a software functional module. It should be noted that the module division in the embodiments of the present invention is illustrative and only represents one logical functional division; other division methods may be used in actual implementation.

[0106] When dividing each function into modules according to its corresponding function. Figure 2 A schematic diagram of the impact damage analysis device for an attitude control power system provided by the present invention is shown. Figure 2 As shown, the device includes:

[0107] The whole machine shock response model construction module 201 is used to construct the whole machine shock response model of the target attitude control power system; the whole machine shock response model includes the finite element model of each component;

[0108] The impact response determination module 202 is used to determine the impact response of the whole machine impact response model using an explicit solver, and obtain the impact response results; the impact response results include the stress response and displacement response of each component;

[0109] The target weak component determination module 203 is used to determine the target weak component based on the stress response of each component in the whole machine impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine;

[0110] The impact damage model construction module 204 is used to construct the impact damage model of the target weak component; the impact damage determination module 205 is used to determine the impact damage of the impact damage model using an explicit solver to obtain the impact damage data of the target weak component; the boundary condition of the impact damage model is the displacement response corresponding to the impact damage model in the impact response results.

[0111] Optionally, the whole-machine impact response model construction module 201 may specifically include:

[0112] A three-dimensional model building unit is used to construct a three-dimensional model of the target attitude control power system based on the geometric structure of the target attitude control power system; the target attitude control power system includes a thin-walled structure, a first force transmission structure, a second force transmission structure, and a heavy-duty structure;

[0113] The component finite element model building unit is used to build finite element models of the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure; the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure are connected and assembled into a whole machine finite element model by binding connection or multi-point constraint.

[0114] The constitutive model setting unit is used to set the material properties of the whole machine finite element model to a non-destructive nonlinear constitutive model, thereby obtaining the whole machine impact response model.

[0115] Optionally, the component finite element model building unit can be specifically used for:

[0116] In the three-dimensional model, the first force transmission structure is equivalent to the central axis, the thin-walled structure is equivalent to the central plane, and the heavy structure is equivalent to the center of mass, thus obtaining a simplified three-dimensional model.

[0117] The simplified 3D model is meshed, the central axis is simulated using beam elements, the center surface is simulated using shell elements, the center of mass is simulated using mass points and assigned a corresponding mass, and the second force transmission structure is simulated using solid elements.

[0118] Optionally, the impact damage model construction module 204 may specifically include:

[0119] The target connection area determination unit is used to extract the target weak component and the model within the target connection area of ​​the target weak component in the whole machine impact response model according to the Saint-Venant principle, so as to obtain the first model.

[0120] The model simplification unit is used to import the original geometric model corresponding to the beam element at the position of the beam element and the original geometric model corresponding to the shell element at the position of the shell element if the first model includes beam elements and shell elements, and simulate it using solid elements to obtain the second model.

[0121] The nonlinear dynamic damage constitutive model setting unit is used to change the material properties in the second model to a strain rate-dependent nonlinear dynamic damage constitutive model to obtain the third model.

[0122] A contact pair setting unit is used to establish contact between the target weak component and the surrounding connecting components in the third model to obtain a fourth model;

[0123] The boundary condition setting unit is used to apply the displacement response of the nodes or surfaces on the boundary of the fourth model in the impact response results to the fourth model to simulate the boundary conditions, so as to obtain the impact damage model of the target weak component.

[0124] Optionally, when the target weak component is a liquid-bearing structure, the impact damage model construction module 204 further includes:

[0125] Fluid domain meshing unit, used to import the fluid domain geometric model into the target weak component of the third model, and to mesh the fluid domain geometric model to obtain fluid domain unit;

[0126] The fluid domain unit attribute setting unit is used to set the fluid domain unit to an arbitrary Lagrange-Euler region and set the attributes of the fluid domain unit to Euler attributes.

[0127] The outer wall constraint setting unit is used to set the outer wall constraint of the target weak component as a Lagrangian constraint, and the constraint at the free liquid surface of the target weak component as an Euler boundary.

[0128] A coupling setting unit is used to couple the contact surface between the target weak component and the fluid domain unit using a coupled Euler-Lagrange method;

[0129] The material property setting unit is used to set the material properties of the fluid domain unit using state equations to obtain a third model.

[0130] Optionally, the target weak component determination module 203 can be specifically used for:

[0131] Components whose stress response exceeds the yield strength are identified as target weak components.

[0132] Optionally, the analysis steps for determining the impact response of the whole machine impact response model and determining the impact damage of the impact damage model are the same.

[0133] Optionally, the material constitutive model of the impact damage model is a nonlinear dynamic damage constitutive model related to strain rate.

[0134] Optionally, when determining impact damage in the impact damage model, if the damage variable of a unit in the impact damage model is 1, the unit with a damage variable of 1 is deleted.

[0135] The above mainly describes the solutions provided by the embodiments of the present invention from the perspective of the interaction between various modules. It is understood that, in order to achieve the above functions, it includes corresponding hardware structures and / or software modules for executing each function. Those skilled in the art should readily recognize that, in conjunction with the units and algorithm steps of the various examples described in the embodiments disclosed herein, the present invention can be implemented in hardware or a combination of hardware and computer software. Whether a function is executed in hardware or by computer software driving hardware depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0136] In the above embodiments, implementation can be achieved entirely or partially through software, hardware, firmware, or any combination thereof. When implemented using software, it can be implemented entirely or partially in the form of a computer program product. The computer program product includes one or more computer programs or instructions. When the computer program or instructions are loaded and executed on a computer, the processes or functions described in the embodiments of the present invention are performed entirely or partially. The computer can be a general-purpose computer, a special-purpose computer, a computer network, a terminal, a user equipment, or other programmable device. The computer program or instructions can be stored in a computer-readable storage medium or transferred from one computer-readable storage medium to another. For example, the computer program or instructions can be transferred from one website, computer, server, or data center to another website, computer, server, or data center via wired or wireless means. The computer-readable storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that integrates one or more available media. The available medium can be a magnetic medium, such as a floppy disk, hard disk, or magnetic tape; it can also be an optical medium, such as a digital video disc (DVD); or it can be a semiconductor medium, such as a solid-state drive (SSD).

[0137] Although the invention has been described herein in conjunction with various embodiments, those skilled in the art will understand and implement other variations of the disclosed embodiments by reviewing the accompanying drawings, the disclosure, and the appended claims in carrying out the claimed invention. In the claims, the word "comprising" does not exclude other components or steps, and "a" or "an" does not exclude a plurality. A single processor or other unit can implement several functions listed in the claims. While different dependent claims may recite certain measures, this does not mean that these measures cannot be combined to produce good results.

[0138] Although the invention has been described in conjunction with specific features and embodiments, it is obvious that various modifications and combinations can be made therein without departing from the spirit and scope of the invention. Accordingly, this specification and drawings are merely exemplary descriptions of the invention as defined by the appended claims, and are considered to cover any and all modifications, variations, combinations, or equivalents within the scope of the invention. Clearly, those skilled in the art can make various alterations and modifications to the invention without departing from its spirit and scope. Thus, if such modifications and modifications of the invention fall within the scope of the claims and their equivalents, the invention is also intended to include such modifications and modifications.

Claims

1. A method for impact damage analysis of an attitude control dynamic system, characterized in that, include: A whole-machine impact response model of the target attitude control power system is constructed; the whole-machine impact response model includes finite element models of each component; An explicit solver is used to determine the impact response of the whole machine impact response model, and the impact response results are obtained; the impact response results include the stress response and displacement response of each component; The target weak component is determined based on the stress response of each component in the overall impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine. Construct an impact damage model for the target weak component; An explicit solver is used to determine the impact damage of the impact damage model, thereby obtaining the impact damage data of the target weak component; the boundary conditions of the impact damage model are the displacement responses corresponding to the impact damage model in the impact response results. The overall impact response model of the target attitude control power system includes: A three-dimensional model of the target attitude control power system is constructed based on the geometric structure of the target attitude control power system; the target attitude control power system includes a thin-walled structure, a first force transmission structure, a second force transmission structure, and a heavy-duty structure; Finite element models are constructed for the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure; the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure are connected and assembled into a whole machine finite element model by binding connection or multi-point constraint. The material properties of the whole machine finite element model are set to a non-destructive nonlinear constitutive model to obtain the whole machine impact response model. The construction of the impact damage model for the target weak component includes: Based on Saint-Venant's principle, the target weak component and the model within the target connection area of ​​the target weak component are extracted from the overall impact response model to obtain the first model; If the first model includes beam elements and shell elements, then import the original geometric model corresponding to the beam element at the position of the beam element, import the original geometric model corresponding to the shell element at the position of the shell element, and use solid elements to simulate to obtain the second model; By changing the material properties in the second model to a strain rate-dependent nonlinear dynamic damage constitutive model, a third model is obtained. The contact between the target weak component and the surrounding connecting components is established in the third model to obtain the fourth model; The displacement response of the nodes or surfaces on the boundary of the fourth model in the impact response results is applied to the fourth model to simulate the boundary conditions, thereby obtaining the impact damage model of the target weak component.

2. The impact damage analysis method for the attitude control dynamic system according to claim 1, characterized in that, The finite element model construction of the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure includes: In the three-dimensional model, the first force transmission structure is equivalent to the central axis, the thin-walled structure is equivalent to the central plane, and the heavy structure is equivalent to the center of mass, thus obtaining a simplified three-dimensional model. The simplified 3D model is meshed, the central axis is simulated using beam elements, the center surface is simulated using shell elements, the center of mass is simulated using mass points and assigned a corresponding mass, and the second force transmission structure is simulated using solid elements.

3. The impact damage analysis method for the attitude control power system according to claim 1, characterized in that, When the target weak component is a liquid-containing structure, the contact between the target weak component and the surrounding connecting components is established in the third model to obtain the fourth model, which also includes the following: A fluid domain geometry model is imported into the target weak component of the third model, and the fluid domain geometry model is meshed to obtain fluid domain elements; The fluid domain element is set as an arbitrary Lagrange-Euler region, and the properties of the fluid domain element are set as Euler properties; The outer wall constraint of the target weak component is set as a Lagrangian constraint, and the constraint at the free liquid surface of the target weak component is set as an Euler boundary. The contact surface between the target weak component and the fluid domain unit is coupled using the coupled Eulerian-Lagrange method; The material properties of the fluid domain element are set using state equations to obtain the third model.

4. The impact damage analysis method for the attitude control dynamic system according to claim 1, characterized in that, Based on the stress response of each component in the overall impact response model, the target weak components are identified as follows: Components whose stress response exceeds the yield strength are identified as target weak components.

5. The impact damage analysis method for the attitude control power system according to claim 1, characterized in that, The analysis steps for determining the impact response of the whole machine impact response model and determining the impact damage of the impact damage model are the same.

6. The impact damage analysis method for the attitude control dynamic system according to claim 1, characterized in that, The material constitutive model of the impact damage model is a nonlinear dynamic damage constitutive model related to strain rate.

7. The impact damage analysis method for the attitude control power system according to claim 1, characterized in that, When determining impact damage in the impact damage model, if the damage variable of a unit in the impact damage model is 1, then the unit with a damage variable of 1 is deleted.

8. An impact damage analysis device for an attitude control power system, characterized in that, include: The whole machine shock response model construction module is used to construct the whole machine shock response model of the target attitude control power system; the whole machine shock response model includes the finite element model of each component; The impact response determination module is used to determine the impact response of the whole machine impact response model using an explicit solver, and obtain the impact response results; the impact response results include the stress response and displacement response of each component; The target weak component determination module is used to determine the target weak component based on the stress response of each component in the whole machine impact response model; the target weak component is the component in the target attitude control power system that undergoes plastic deformation under the impact of the whole machine; The impact damage model construction module is used to construct the impact damage model of the target weak component; the material constitutive model of the impact damage model is a nonlinear dynamic damage constitutive model related to strain rate. The impact damage determination module is used to determine the impact damage of the impact damage model using an explicit solver, thereby obtaining the impact damage data of the target weak component; the boundary conditions of the impact damage model are the displacement responses corresponding to the impact damage model in the impact response results. The whole machine shock response model construction module specifically includes: A three-dimensional model building unit is used to construct a three-dimensional model of the target attitude control power system based on the geometric structure of the target attitude control power system; the target attitude control power system includes a thin-walled structure, a first force transmission structure, a second force transmission structure, and a heavy-duty structure; The component finite element model building unit is used to build finite element models of the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure; the thin-walled structure, the first force transmission structure, the second force transmission structure, and the heavy structure are connected and assembled into a whole machine finite element model by binding connection or multi-point constraint. Constitutive model setting unit is used to set the material properties of the whole machine finite element model to a non-destructive nonlinear constitutive model to obtain the whole machine impact response model; The impact damage model construction module specifically includes: The target connection area determination unit is used to extract the target weak component and the model within the target connection area of ​​the target weak component in the whole machine impact response model according to the Saint-Venant principle, so as to obtain the first model. The model simplification unit is used to import the original geometric model corresponding to the beam element at the position of the beam element and the original geometric model corresponding to the shell element at the position of the shell element if the first model includes beam elements and shell elements, and simulate it using solid elements to obtain the second model. The nonlinear dynamic damage constitutive model setting unit is used to change the material properties in the second model to a strain rate-dependent nonlinear dynamic damage constitutive model to obtain the third model. A contact pair setting unit is used to establish contact between the target weak component and the surrounding connecting components in the third model to obtain a fourth model; The boundary condition setting unit is used to apply the displacement response of the nodes or surfaces on the boundary of the fourth model in the impact response results to the fourth model to simulate the boundary conditions, so as to obtain the impact damage model of the target weak component.

Images