Equivalent mechanical modeling method and system for variable-thickness composite material component

Abstract

The invention belongs to the field of mechanical simulation analysis of composite material structural parts, and particularly discloses an equivalent mechanical modeling method of a variable-thickness composite material component, which comprises the following steps: acquiring a laying deformation measure, and determining deformation distribution of fiber cloth by using the laying deformation measure through a fishing net method; mapping the deformation distribution of the fiber cloth to a target structure grid by using a grid mapping method; constructing a numerical model of equivalent mechanical characteristics of the target structure grid by using a multilayer structure equivalent calculation formula; the equivalent mechanical characteristics comprise equivalent density, an equivalent stiffness matrix and an equivalent damping matrix; and determining a numerical value implementation program of the numerical value model, and performing equivalent mechanical analysis simulation on the variable-thickness composite material component by the numerical value model through the numerical value implementation program by adopting a finite element method. According to the method, the modeling accuracy of the variable-thickness composite material component is improved.

Description

Technical Field

[0001] The present application belongs to the field of mechanical simulation analysis of composite structural parts, and more specifically, relates to an equivalent mechanical modeling method and system for thick composite structural parts. Background Art

[0002] Fiber-reinforced composite materials are usually applied in the form of laminates, which are stacked in a certain order to form different composite components such as plates, columns, shells or beams. Due to the anisotropy of its single-layer material, single layers with different fiber orientations are combined at different angles to meet various required performance requirements. The ply design of composite components is divided into two types: constant thickness and variable thickness. When designing variable thickness composite components, a certain number of plies will be added or discarded at specific locations. By changing the local ply thickness or fiber orientation, the structure can achieve different stiffness changes at different locations. It generally has the characteristics of uneven thickness, decreasing plies, and complex orientation. In the area where the thickness of the variable thickness component changes, the orientation of adjacent plies is not parallel, and a resin-filled area will be formed.

[0003] Variable-thickness composite components are commonly used in large structures such as aerospace, ship propulsion, and automotive manufacturing. Differential loads between different regions of a component often lead to different thickness designs, known as "drop-off" structures. This variable-thickness design not only offers design diversity and idealization, but also presents challenges such as difficult modeling and analysis, and increased manufacturing costs.

[0004] At present, the existing modeling methods cannot fully consider the complexity of the shape of variable thickness composite components and the factors affecting the draping process, making it difficult to accurately perform mechanical modeling of variable thickness composite components.

[0005] Therefore, how to improve the accuracy of modeling of variable thickness composite components is a technical problem that needs to be solved urgently. Summary of the Invention

[0006] In response to the defects of the existing technology, the purpose of this application is to provide an equivalent mechanical modeling method and system for thickness composite components, aiming to solve the problem in the existing technology that the modeling accuracy is low due to the inability to fully consider the complexity of the shape of variable thickness composite components and the influencing factors of the draping process.

[0007] To achieve the above objectives, in a first aspect, the present application provides an equivalent mechanical modeling method for a variable thickness composite material component, comprising: Obtaining draping deformation metrics, and using the draping deformation metrics to determine the deformation distribution of the fiber cloth using the fishnet method, and mapping the deformation distribution of the fiber cloth to the target structure grid using the grid mapping method; A numerical model of the equivalent mechanical characteristics of the target structure grid is constructed using an equivalent calculation formula for a multi-layer structure; the equivalent mechanical characteristics include an equivalent density, an equivalent stiffness matrix, and an equivalent damping matrix; the equivalent calculation formula for the multi-layer structure is constructed based on the parameters of the equivalent mechanical characteristics, which are obtained by weighted calculation of the stiffness, damping, and density of each ply by combining ply information with coordinate system transformation.

[0008] A numerical implementation program of the numerical model is determined, and the numerical model is passed through the numerical implementation program to perform equivalent mechanical analysis simulation on a variable thickness composite material component using a finite element method.

[0009] Optionally, obtaining the draping deformation metric, determining the fiber cloth deformation distribution by using the draping deformation metric through a fishnet method, and mapping the fiber cloth deformation distribution to a target structure grid by using a grid mapping method includes: Using the fishing net method to perform draping simulation, determine the degree of deformation of the fiber cloth during the draping process, and obtain a deformation measure of the fiber cloth draping according to the deformation degree; Determine the draping simulation result based on the deformation metric, and determine the deformation distribution of the fiber cloth in the variable thickness composite material component; the fiber cloth deformation distribution includes the spatial shell unit grid and the fiber deformation metric distribution after the fiber cloth is draped; The shell-solid mesh mapping method is used to map the deformation distribution of the fiber cloth from the spatial shell element mesh to the target structural mesh of the variable thickness composite component for the hexahedral mesh used in mechanical simulation.

[0010] Optionally, the draping deformation metric includes: a warp shear angle, a weft shear angle, a normal unit vector, a shear deformation reference unit vector, and a fiber cloth thickness.

[0011] Optionally, the implementation process of the shell-solid mesh mapping method includes: Constructing an initialization matrix of size m×n, wherein the initialization matrix is ​​used to store the intersection relationship between plies and solid elements; m represents the number of solid elements in the target structure grid, and n represents the number of plies in the space shell grid; For each solid element in the target structural mesh, traverse the triangular mesh of the target ply of the spatial shell element mesh; Using a spatial triangle-spatial hexahedron collision detection algorithm to determine whether the solid unit and the target ply intersect, marking the empty matrix according to the determination result to obtain a marking matrix; For each entity cell in the target structure grid, searching for non-zero entries in the target row of the label matrix; According to the search results, the warp shear angle, weft shear angle, normal unit vector, shear deformation reference unit vector, and fiber cloth thickness of the target ply at the solid element are written into the solid element.

[0012] Optionally, the numerical model for determining the equivalent mechanical properties of the target structure grid using a multi-layer structure equivalent calculation formula includes: Determining the ply order of the solid elements according to the normal element vectors of the solid elements of the target structural network; Determine the warp unit vector and the weft unit vector of each layer of fiber cloth based on the ply order, the warp shear angle, the weft shear angle, the normal unit vector, and the shear deformation reference unit vector; Determine the fiber cloth density, elastic stiffness matrix and material damping matrix of each layer of fiber cloth; Based on the ply sequence, the equivalent mechanical properties corresponding to the solid unit are determined in combination with the warp unit vector, weft unit vector, fiber cloth thickness and density, elastic stiffness matrix and material damping matrix to construct the numerical model.

[0013] Optionally, based on the ply sequence, the equivalent mechanical properties corresponding to the solid element are determined in combination with the warp unit vector, the weft unit vector, the fiber cloth thickness, the fiber cloth density, the elastic stiffness matrix, and the material damping matrix to construct the numerical model, including: Establish a local coordinate system corresponding to the number of plies; Obtaining a normal unit vector of the entity element, and using the first local coordinate in the unit direction of the normal vector as the local coordinate system of the entity element; Determine the transverse coordinate axis of the local coordinate system according to the warp direction unit vector and the weft direction unit vector, use the normal vector as the height coordinate axis of the local coordinate system, and rotate the height coordinate axis according to the right-hand rule to obtain the longitudinal coordinate axis of the local coordinate system; Establishing a transformation matrix between each ply and the solid element, wherein the transformation matrix is ​​used to represent the transformation relationship between stress or strain between the local coordinate system and the global coordinate system; Convert the elastic stiffness matrix and material damping matrix of each ply from the local coordinate system to the global coordinate system to obtain the stiffness matrix and damping matrix of each ply in the global coordinate system; The equivalent stiffness matrix of the solid element is calculated by weighted average based on the stress continuity condition, the strain continuity condition, the volume fraction of each ply, and the converted stiffness matrix. The volume fraction is determined based on the ply material thickness and the total thickness of the solid element. Based on the converted damping matrix, the equivalent stiffness matrix and the equivalent flexibility matrix are combined to calculate and obtain the equivalent damping matrix of the solid element; the equivalent flexibility matrix is ​​the inverse matrix of the equivalent stiffness matrix; Based on the fiber cloth thickness and fiber cloth density of each layer, the equivalent density of the solid unit is calculated by volume weighting; The equivalent stiffness matrix, equivalent damping matrix and equivalent density are used as equivalent mechanical characteristic parameters of the solid element to construct a numerical model.

[0014] Optionally, the implementation process of the numerical model includes: Perform a fishnet draping simulation and read the simulation results to determine the deformation distribution of the fiber cloth and the deformation measurement data of each draping; Using a programming interface of finite element analysis software to read the target structure grid data and ply material performance data; Using the shell-solid mesh mapping method, the entity elements of the target structure mesh are traversed to complete the mesh data mapping; The equivalent calculation formula of multi-layer structure is used to calculate the equivalent mechanical properties of the target structure grid; The programming interface function of the finite element analysis software is used to assign attributes to the target structure mesh.

[0015] In a second aspect, the present application further provides an equivalent mechanical modeling system for a variable thickness composite material component, comprising: a deformation distribution determination module for obtaining a draping deformation metric and determining a fiber cloth deformation distribution using the draping deformation metric through a fishnet method; A mapping module, used for mapping the deformation distribution of the fiber cloth to the target structure grid using a grid mapping method; A model construction module is used to construct a numerical model of the equivalent mechanical characteristics of the target structure grid using an equivalent calculation formula for a multi-layer structure; the equivalent mechanical characteristics include an equivalent density, an equivalent stiffness matrix, and an equivalent damping matrix; the equivalent calculation formula for a multi-layer structure is constructed based on the parameters of the equivalent mechanical characteristics, which are obtained by weighted calculation of the stiffness, damping, and density of each ply by combining ply information with coordinate system transformation.

[0016] The simulation module is used to determine the numerical implementation program of the numerical model, and use the finite element method to perform equivalent mechanical analysis simulation on the variable thickness composite material component through the numerical model.

[0017] In a third aspect, the present application provides an electronic device comprising: at least one memory for storing programs; and at least one processor for executing the programs stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method described in the first aspect or any possible implementation of the first aspect.

[0018] In a fourth aspect, the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method described in the first aspect or any possible implementation of the first aspect.

[0019] In a fifth aspect, the present application provides a computer program product, which, when executed on a processor, enables the processor to execute the method described in the first aspect or any possible implementation of the first aspect.

[0020] It can be understood that the beneficial effects of the second to fifth aspects mentioned above can be found in the relevant description of the first aspect mentioned above, and will not be repeated here.

[0021] In general, the above technical solutions conceived by this application have the following beneficial effects compared with the existing technologies: (1) This application combines the draping deformation measurement with the fishing net method to accurately capture the deformation distribution of the fiber cloth during the draping process, and uses the grid mapping method to accurately transfer it to the target structure grid, ensuring the fidelity of the initial deformation information. A numerical model containing equivalent density, stiffness and damping is constructed using the multi-layer structure equivalent calculation formula, which can effectively simulate the complex mechanical behavior of composite laminated structures. Finally, the analysis and simulation of variable thickness composite components by combining the numerical implementation program with the finite element method can not only overcome the difficulties of traditional methods in dealing with complex geometry and variable thickness problems, and comprehensively consider the influencing factors in the draping process, but also more accurately reflect the response of variable thickness composite components under actual working conditions, thereby significantly improving the accuracy of modeling.

[0022] (2) This application uses a fishing net method to simulate the deformation distribution of fiber cloth draping, which can fully consider the influence of the forming process of variable thickness composite components. In addition, the use of solid meshing can preserve the geometric shape during the simulation process, and the solid element supports three-dimensional stress-strain simulation, which greatly improves the accuracy of mechanical modeling of variable thickness composite components.

[0023] (3) This application establishes a corresponding numerical model from the perspective of equivalent mechanical modeling, and uses the program development interface of the finite element framework to complete the numerical implementation program. It is compatible with the general finite element framework and automates the calculation of equivalent mechanical characteristics, reducing the difficulty of equivalent mechanical modeling and analysis. It has the advantages of high computational efficiency and accuracy, and accelerates the design of variable thickness composite components. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] Figure 1 This is a flow chart of an equivalent mechanical modeling method for a variable thickness composite material component provided in an embodiment of the present application; Figure 2 This is a flow chart of the numerical model of the variable thickness composite material component constructed in the embodiment of the present application, which is numerically realized through the Abaqus & Python secondary development interface; Figure 3 is a schematic diagram of an equivalent mechanical finite element model of a composite propeller constructed according to a preferred embodiment of the present application; Figure 4 It is a vibration shape cloud diagram and modal frequency distribution diagram of the wet modal analysis of the composite propeller constructed according to the preferred embodiment of the present application; Figure 5 This is a schematic structural diagram of an equivalent mechanical modeling system for a variable thickness composite material component provided in an embodiment of the present application; Figure 6 It is a structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0025] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.

[0026] The term "and / or" as used herein describes an association between related objects, indicating that three possible relationships exist. For example, "A and / or B" can represent: A exists alone, A and B exist simultaneously, or B exists alone. The symbol " / " as used herein indicates that the related objects are in an "or" relationship, for example, A / B means either A or B.

[0027] The terms "first" and "second" in this specification and claims are used to distinguish different objects rather than to describe a specific order of objects. For example, "first response message" and "second response message" are used to distinguish different response messages rather than to describe a specific order of response messages.

[0028] In the embodiments of this application, words such as "exemplary" or "for example" are used to indicate examples, illustrations, or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of this application should not be interpreted as being preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary" or "for example" is intended to present the relevant concepts in a concrete manner.

[0029] In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more, for example, multiple processing units means two or more processing units, etc.; multiple elements means two or more elements, etc.

[0030] The embodiments of the present application are described below in conjunction with the drawings in the embodiments of the present application.

[0031] Reference Figure 1 The present application provides an equivalent mechanical modeling method for a variable thickness composite material component, comprising: S101. Obtaining a draping deformation metric, and using the draping deformation metric to determine the fiber cloth deformation distribution by the fishnet method; S102. Mapping the deformation distribution of the fiber cloth to the target structure grid using a grid mapping method; S103. Construct a numerical model of equivalent mechanical properties of the target structure grid using a multi-layer structure equivalent calculation formula; the equivalent mechanical properties include equivalent density, equivalent stiffness matrix, and equivalent damping matrix; S104. Determine a numerical implementation program for the numerical model, and use the finite element method to perform equivalent mechanical analysis simulation on the variable thickness composite material component using the numerical model through the numerical implementation program.

[0032] Specifically, the present application obtains draping deformation metrics, uses the draping deformation metrics to determine the fiber cloth deformation distribution through the fishnet method, and uses the grid mapping method to map the fiber cloth deformation distribution to the target structure grid.

[0033] First, the local deformation distribution of the fiber fabric is quantified by accurately measuring or calculating the deformation of the fiber fabric due to bending, stretching, or shearing during the draping process, combined with the fishnet method. Optionally, the draping deformation metric includes: warp shear angle, weft shear angle, normal unit vector, shear deformation reference unit vector, and fiber fabric thickness.

[0034] It should be noted that the fishing net method in the embodiment of the present application is a numerical simulation method for analyzing the deformation distribution of fiber cloth during the draping process. Its core idea is to discretize the fiber cloth into a mesh structure composed of nodes and line segments, and quantify the local deformation of the fiber cloth by simulating the deformation of the grid under force or process constraints.

[0035] After obtaining the deformation distribution of the fiber cloth, the mesh mapping method is used to accurately transfer this deformation information from the original mesh of the fiber cloth to the finite element mesh of the target structure, ensuring that the target structure mesh can reflect the actual deformation state of the fiber cloth and provide accurate geometric initial conditions for subsequent mechanical analysis.

[0036] Among them, the target structure in the embodiment of the present application is a hexahedral grid system of a variable thickness composite material component used for subsequent mechanical simulation.

[0037] Furthermore, a numerical model of the equivalent mechanical properties of the target structure grid is constructed using the multi-layer structure equivalent calculation formula. After obtaining the deformed target structure grid, this step uses the multi-layer structure equivalent calculation formula to perform equivalent processing to obtain the equivalent mechanical properties of the entire structure or local area. Specifically, it includes calculating the equivalent density, constructing the equivalent stiffness matrix, and determining the equivalent damping matrix. These equivalent parameters can comprehensively reflect the overall mechanical response of the multi-layer composite material. In this way, the complex anisotropic multi-layer structure is simplified into a numerical model with equivalent mechanical properties, which not only retains the key mechanical behavior of the material but also reduces the complexity of the model.

[0038] Finally, a numerical implementation program for the numerical model is determined, and the numerical model is passed through the numerical implementation program to perform equivalent mechanical analysis simulation on the variable thickness composite component using the finite element method. Optionally, the numerical implementation program is implemented using the Python interface of Abaqus.

[0039] A numerical model containing equivalent density, stiffness matrix, and damping matrix was imported into the program, and the finite element method (FEM) was used to perform mechanical analysis and simulation of the variable-thickness composite component. This FEM analysis yielded equivalent mechanical analysis results, leveraging the precise deformation information and equivalent mechanical properties provided by the previous steps to ensure that the simulation results accurately reflect the actual mechanical behavior of the variable-thickness composite component.

[0040] The implementation tool of this embodiment is Abaqus' Python secondary development interface. The program implements the numerical model by programming a series of Python module functions (including a fishnet draping simulation function, a 3D collision detection function, a shell-to-solid mesh mapping function, a mechanical property equivalence function, and an equivalent attribute assignment function) to facilitate integration with the finite element solver, including consideration of fiber draping deformation and component equivalent modeling. The fishnet method is used to calculate fiber draping deformation data; a hierarchical bounding box data structure is combined with a collision detection algorithm to efficiently map draping deformation data to a mesh; a multilayer structure equivalent calculation formula is used to solve the equivalent mechanical properties of variable-thickness composite components; and the Python secondary development interface is used to implement automated finite element modeling.

[0041] Through the established numerical model and finite element software, the deformation, stress, mode and other results of variable thickness composite components under different loading / boundary conditions are obtained.

[0042] Using established numerical implementation methods, the numerical model was placed within the Python and finite element framework to perform equivalent mechanical modeling of variable-thickness composite components, enabling simulation analysis of deformation, stress, and modal characteristics. By considering the fiber deformation distribution during the draping process and retaining the three-dimensional mesh system, the accuracy and efficiency of modeling of variable-thickness composite components were effectively improved.

[0043] Optionally, obtaining the draping deformation metric, determining the fiber cloth deformation distribution by using the draping deformation metric through a fishnet method, and mapping the fiber cloth deformation distribution to a target structure grid by using a grid mapping method includes: Using the fishing net method to perform draping simulation, determine the degree of deformation of the fiber cloth during the draping process, and obtain a deformation measure of the fiber cloth draping according to the deformation degree; Determine the draping simulation result based on the deformation metric, and determine the deformation distribution of the fiber cloth in the variable thickness composite material component; the fiber cloth deformation distribution includes the spatial shell unit grid and the fiber deformation metric distribution after the fiber cloth is draped; The shell-solid mesh mapping method is used to map the deformation distribution of the fiber cloth from the spatial shell element mesh to the target structural mesh of the variable thickness composite component for the hexahedral mesh used in mechanical simulation.

[0044] Specifically, the fishing net method of the embodiment of the present application abstracts the fiber cloth into a grid model composed of fiber paths, simulates the deformations such as stretching, bending and shearing caused by geometric constraints in the process of laying the fiber cloth on a complex curved surface component, thereby quantifying the degree of deformation of the fiber cloth in different areas and obtaining the deformation measurement of the fiber cloth laying accordingly.

[0045] Based on the obtained deformation metrics, detailed draping simulation results can be further determined to determine the deformation distribution of the fiber cloth within the variable-thickness composite component. This distribution is specifically manifested as the spatial shell unit grid shape formed by the fiber cloth within the component after draping, and the corresponding fiber deformation metric distribution.

[0046] Using a shell-solid mesh mapping method, the deformation distribution of the fiber cloth is accurately mapped from its own spatial shell element mesh to the target structural mesh, typically composed of hexahedral elements, prepared for the mechanical analysis of variable-thickness composite components. This mapping ensures that the mesh used in the mechanical simulation is not only geometrically accurate but also contains true information about the initial deformation of the fiber cloth.

[0047] Optionally, the implementation process of the shell-solid mesh mapping method includes: Constructing an initialization matrix of size m×n, wherein the initialization matrix is ​​used to store the intersection relationship between plies and solid elements; m represents the number of solid elements in the target structure grid, and n represents the number of plies in the space shell grid; For each solid element in the target structural mesh, traverse the triangular mesh of the target ply of the spatial shell element mesh; Using a spatial triangle-spatial hexahedron collision detection algorithm to determine whether the solid unit and the target ply intersect, marking the empty matrix according to the determination result to obtain a marking matrix; For each entity cell in the target structure grid, searching for non-zero entries in the target row of the label matrix; According to the search results, the warp shear angle, weft shear angle, normal unit vector, shear deformation reference unit vector, and fiber cloth thickness of the target ply at the solid element are written into the solid element.

[0048] Specifically, in this embodiment, first initialize an empty matrix D with a matrix size of (m, n); For the solid element H in the target structure mesh i The position information of the space shell grid is traversed j The triangular mesh of the layer is judged whether it intersects. If it intersects, D ij Marked as 1, otherwise D ij is 0; wherein the intersection judgment is realized by a spatial triangle-space hexahedron collision detection algorithm; For the solid element H in the target structure mesh i , find the non-zero entries in the i-th row of the D matrix. If D ij =1, then it is a solid element H i Write layer P j Warp / weft shear angle, normal unit vector, shear deformation reference unit vector, and fabric thickness at the solid element.

[0049] Optionally, the numerical model for determining the equivalent mechanical properties of the target structure grid using a multi-layer structure equivalent calculation formula includes: Determining the ply order of the solid elements according to the normal element vectors of the solid elements of the target structural network; Determine the warp unit vector and the weft unit vector of each layer of fiber cloth based on the ply order, the warp shear angle, the weft shear angle, the normal unit vector, and the shear deformation reference unit vector; Determine the fiber cloth density, elastic stiffness matrix and material damping matrix of each layer of fiber cloth; Based on the ply sequence, the equivalent mechanical properties corresponding to the solid unit are determined in combination with the warp unit vector, weft unit vector, fiber cloth thickness and density, elastic stiffness matrix and material damping matrix to construct the numerical model.

[0050] Specifically, the present embodiment first determines the order of composite material plies within a solid unit of the target structural network by analyzing the normal unit vectors of the solid unit. Specifically, the normal vector indicates the orientation of the unit, and the stacking order of the various layers of fiber cloth contained in the unit can be inferred from the normal vector.

[0051] After the ply order is determined, the actual warp and weft unit vectors of each layer of fiber cloth in the current unit are calculated based on the order, the initial shear angles of the warp and weft yarns of the fiber cloth, the normal unit vector of the solid unit, and the shear deformation reference unit vector.

[0052] Finally, based on the determined ply order, the warp unit vector, weft unit vector, fiber cloth thickness, density, elastic stiffness matrix, and material damping matrix of each layer are comprehensively calculated to determine the overall equivalent mechanical properties of the solid element. The equivalent properties usually include equivalent density, equivalent stiffness matrix, and equivalent damping matrix.

[0053] Optionally, based on the ply sequence, the equivalent mechanical properties corresponding to the solid element are determined in combination with the warp unit vector, the weft unit vector, the fiber cloth thickness, the fiber cloth density, the elastic stiffness matrix, and the material damping matrix to construct the numerical model, including: Establish a local coordinate system corresponding to the number of plies; Obtaining a normal unit vector of the entity element, and using the first local coordinate in the unit direction of the normal vector as the local coordinate system of the entity element; Determine the transverse coordinate axis of the local coordinate system according to the warp direction unit vector and the weft direction unit vector, use the normal vector as the height coordinate axis of the local coordinate system, and rotate the height coordinate axis according to the right-hand rule to obtain the longitudinal coordinate axis of the local coordinate system; Establishing a transformation matrix between each ply and the solid element, wherein the transformation matrix is ​​used to represent the transformation relationship between stress or strain between the local coordinate system and the global coordinate system; Convert the elastic stiffness matrix and material damping matrix of each ply from the local coordinate system to the global coordinate system to obtain the stiffness matrix and damping matrix of each ply in the global coordinate system; The equivalent stiffness matrix of the solid element is calculated by weighted average based on the stress continuity condition, the strain continuity condition, the volume fraction of each ply, and the converted stiffness matrix. The volume fraction is determined based on the ply material thickness and the total thickness of the solid element. Based on the converted damping matrix, the equivalent stiffness matrix and the equivalent flexibility matrix are combined to calculate and obtain the equivalent damping matrix of the solid element; the equivalent flexibility matrix is ​​the inverse matrix of the equivalent stiffness matrix; Based on the fiber cloth thickness and fiber cloth density of each layer, the equivalent density of the solid unit is calculated by volume weighting; The equivalent stiffness matrix, equivalent damping matrix and equivalent density are used as equivalent mechanical characteristic parameters of the solid element to construct a numerical model.

[0054] Specifically, in this embodiment, a local coordinate system independent of the overall structure coordinate system is first established for each ply.

[0055] Next, the warp yarn direction unit vector and the weft yarn direction unit vector of the ply determined in the previous step are used to determine the x-axis of the local coordinate system, the normal unit vector is used as the z-axis, and the warp yarn is rotated 90° around the z-axis according to the right-hand rule to obtain the y-axis.

[0056] Then, based on the established relationship between the local coordinate system and the global structure coordinate system, a transformation matrix is ​​constructed for each ply. The transformation matrix accurately describes the transformation relationship between the stress or strain tensor in the local coordinate system and the global coordinate system, ensuring that mechanical quantities in different coordinate systems can be mapped and calculated.

[0057] The original elastic stiffness matrix and material damping matrix of each ply defined in the local coordinate system are converted to the global coordinate system using the conversion matrix to obtain the stiffness matrix and damping matrix of each ply in the global coordinate system.

[0058] After obtaining the mechanical parameters of all plies in the global coordinate system, based on the stress continuity conditions and strain continuity conditions, combined with the volume fraction of each ply (the volume fraction is determined by the ratio of the thickness of the ply to the total thickness of the solid unit), and the converted stiffness matrix, the equivalent stiffness matrix of the entire solid unit is calculated by the weighted average method.

[0059] Similarly, the converted damping matrix is ​​combined with the just-calculated equivalent stiffness matrix and its inverse (i.e., the equivalent flexibility matrix) to determine the equivalent damping matrix of the solid element using a specific formula. The equivalent density of the solid element is calculated using a volume-weighted average based on the fiber cloth thickness and density of each ply.

[0060] Finally, the three key parameters of the calculated equivalent stiffness matrix, equivalent damping matrix and equivalent density are used as the equivalent mechanical characteristic parameters of the solid element.

[0061] Specifically, the equivalent calculation formula process of the multi-layer structure described in this embodiment is as follows: For solid element Hi , assuming that it corresponds to N layers, there are a series of local coordinate systems ; Take H i The normal unit vector of the ground element is ,Will The first local coordinate in the direction As unit H i The local coordinate system of 1) The transformation relationship between the stiffness / damping matrix and different coordinate systems includes: Define the R matrix: ; The T matrix is ​​used to represent the transformation relationship between stress / strain in the local coordinate system and the global coordinate system. The shear strain in is the engineering shear strain:

[0062] in:

[0063] Global coordinate system and local coordinate system The direction cosines of each coordinate axis between them are shown in the following table:

[0064] Based on the above, the expression of the material stiffness matrix in the global coordinate system can be obtained:

[0065] Also, the stiffness matrix is ​​converted from the local coordinate system To the local coordinate system The conversion relationship:

[0066] in, are the stiffness matrices of the corresponding coordinate systems, are the T matrices of the corresponding coordinate systems respectively.

[0067] In one vibration cycle, the dissipated energy of the system is:

[0068] in, is the material damping matrix, which represents the damping capacity of the material itself. This matrix is ​​determined by experiment. Similar to the stiffness matrix, it is easy to obtain the damping matrix from the local coordinate system. To the local coordinate system The conversion relationship:

[0069] in, are the damping matrices of the corresponding coordinate systems.

[0070] 2) The calculation process of the equivalent stiffness matrix is ​​as follows: No. Layers in the material coordinate system The stress-strain relationship under satisfies the orthotropic anisotropy:

[0071] Definition of effective stress and effective strain of multilayer structure:

[0072] Where V is the unit volume of the sublayer , the stress and strain of each single layer in the sublayer unit remain uniform in the in-plane direction. Further, the effective stress and effective strain are expressed as:

[0073] and It is Layer stress and strain, represents the volume fraction. Considering the continuity conditions of stress and strain of the laminate, we have:

[0074] After complex derivation, the equivalent stiffness matrix satisfy:

[0075]

[0076]

[0077]

[0078] 3) The calculation process of the equivalent damping matrix is ​​as follows: For a multilayer structure with N layers, the dissipated energy per unit volume during one vibration cycle can be expressed as:

[0079] Dissipated energy density used 、 、 Express:

[0080] Damping Matrix Convert to After the coordinate system, the equivalent damping matrix can be expressed as:

[0081] in, is the equivalent flexibility matrix and the equivalent stiffness matrix, and are as follows:

[0082]

[0083]

[0084]

[0085]

[0086]

[0087] 4) The equivalent density calculation process is as follows: The equivalent density calculation formula is:

[0088] Optionally, the implementation process of the numerical model includes: Perform a fishnet draping simulation and read the simulation results to determine the deformation distribution of the fiber cloth and the deformation measurement data of each draping; Using a programming interface of finite element analysis software to read the target structure grid data and ply material performance data; Using the shell-solid mesh mapping method, the entity elements of the target structure mesh are traversed to complete the mesh data mapping; The equivalent calculation formula of multi-layer structure is used to calculate the equivalent mechanical properties of the target structure grid; The programming interface function of the finite element analysis software is used to assign attributes to the target structure mesh.

[0089] The implementation process is implemented through Abaqus, specifically: The deformation distribution of the fiber cloth and various draping deformation metrics are read from the draping simulation results using the fishing net method. The target structure mesh data and ply material performance data are read using the Abaqus Python interface. The shell-solid mesh mapping method is used to traverse the solid elements of the target structure mesh and complete the mesh data mapping. The multi-layer structure equivalent calculation formula is used to calculate the equivalent mechanical properties of the target structure mesh. The Abaqus Python interface function is used to assign attributes to the target structure mesh.

[0090] Specifically, the numerical model is implemented in the framework of the Abaqus&Python interface function, and the program implementation process is as follows: Reference Figure 2 The numerical model has the following numerical implementation process: the fishing net method is used to calculate the deformation data of the fiber cloth draping, and then the ply material performance data, draping deformation results, and target structure mesh data are read respectively; the three types of data are passed into the shell-solid mesh mapping function for mesh mapping; the equivalent mechanical properties of the target structure mesh are calculated, and finally the Abaqus Python secondary development interface is called for automated equivalent mechanical modeling.

[0091] The present application is described in detail below with reference to specific embodiments.

[0092] Example 1: The equivalent mechanical modeling method of variable thickness composite components proposed in this application was used to perform wet modal simulation of the composite propeller DTMB4119 component.

[0093] like Figure 3 As shown in Figure 1, the finite element model of the propeller component is established using the equivalent mechanical modeling method in Abaqus software. Figure 3 As shown, the wet modal analysis includes the structural domain (propeller) and the fluid domain (water). As described above, the numerical model and numerical implementation program proposed in this application are used, the fishing net method is used for draping simulation, and the numerical implementation program is used to assign properties of equivalent mechanical characteristics to the propeller mesh. A spherical fluid domain with a diameter of 3000mm is used, the acoustic element type is AC3D10, and the acoustic infinite element ACIN3D6 is used to simulate the infinite fluid domain. ASI elements are used to couple the displacement and pressure degrees of freedom between the solid and fluid. The Abaqus acoustic-solid coupling module is used to solve the first 40 modes. Fixed displacement boundary conditions (U1=U2=U3=UR1=UR2=UR3=0) are applied to the bottom node of the propeller hub.

[0094] Reference Figure 4 , Figure 4 The vibration mode cloud diagrams and modal frequency distributions of the first four wet modes completed using the equivalent mechanical modeling method proposed in this application are shown. Figure 4 In the propeller, due to its rotationally symmetrical structure, there are multiple vibration modes at certain frequencies. However, compared with the dry mode vibration mode, the vibration modes between the blades are not completely consistent. This is caused by the added mass effect of the fluid. In the first-order vibration mode, all blades twist around their respective neutral lines; in the second-order vibration mode, the blades deform around the hub; in the third-order vibration mode, the lower blade undergoes obvious bending deformation, with an amplitude of about half that of the other two blades; the fourth-order vibration mode is the opposite of the third-order vibration mode, with the lower blade almost not vibrating, while the upper blade undergoes bending deformation. There are also multiple identical vibration modes in the subsequent vibration modes, but the difference is that the blades are different.

[0095] Reference Figure 5The present application also provides an equivalent mechanical modeling system for a variable thickness composite material component, comprising: a deformation distribution determination module 510 for obtaining a draping deformation metric and determining a fiber cloth deformation distribution using the draping deformation metric by a fishnet method; A mapping module 520 is used to map the deformation distribution of the fiber cloth to the target structure grid using a grid mapping method; Model construction module 530 is used to construct a numerical model of the equivalent mechanical characteristics of the target structure grid using a multi-layer structure equivalent calculation formula; the equivalent mechanical characteristics include equivalent density, equivalent stiffness matrix, and equivalent damping matrix; the multi-layer structure equivalent calculation formula is constructed based on the parameters of the equivalent mechanical characteristics, and the parameters of the equivalent mechanical characteristics are obtained by combining ply information with coordinate system transformation and performing a weighted calculation of the stiffness, damping, and density of each ply.

[0096] The simulation module 540 is used to determine a numerical implementation program of the numerical model, and to perform equivalent mechanical analysis simulation on a variable thickness composite material component by using the finite element method through the numerical model.

[0097] Optionally, obtaining the draping deformation metric, determining the fiber cloth deformation distribution by using the draping deformation metric through a fishnet method, and mapping the fiber cloth deformation distribution to a target structure grid by using a grid mapping method includes: Using the fishing net method to perform draping simulation, determine the degree of deformation of the fiber cloth during the draping process, and obtain a deformation measure of the fiber cloth draping according to the deformation degree; Determine the draping simulation result based on the deformation metric, and determine the deformation distribution of the fiber cloth in the variable thickness composite material component; the fiber cloth deformation distribution includes the spatial shell unit grid and the fiber deformation metric distribution after the fiber cloth is draped; The shell-solid mesh mapping method is used to map the deformation distribution of the fiber cloth from the spatial shell element mesh to the target structural mesh of the variable thickness composite component for the hexahedral mesh used in mechanical simulation.

[0098] Optionally, the draping deformation metric includes: a warp shear angle, a weft shear angle, a normal unit vector, a shear deformation reference unit vector, and a fiber cloth thickness.

[0099] Optionally, the implementation process of the shell-solid mesh mapping method includes: Constructing an initialization matrix of size m×n, wherein the initialization matrix is ​​used to store the intersection relationship between plies and solid elements; m represents the number of solid elements in the target structure grid, and n represents the number of plies in the space shell grid; For each solid element in the target structural mesh, traverse the triangular mesh of the target ply of the spatial shell element mesh; Using a spatial triangle-spatial hexahedron collision detection algorithm to determine whether the solid unit and the target ply intersect, marking the empty matrix according to the determination result to obtain a marking matrix; For each entity cell in the target structure grid, searching for non-zero entries in the target row of the label matrix; According to the search results, the warp shear angle, weft shear angle, normal unit vector, shear deformation reference unit vector, and fiber cloth thickness of the target ply at the solid element are written into the solid element.

[0100] Optionally, the numerical model for determining the equivalent mechanical properties of the target structure grid using a multi-layer structure equivalent calculation formula includes: Determining the ply order of the solid elements according to the normal element vectors of the solid elements of the target structural network; Determine the warp unit vector and the weft unit vector of each layer of fiber cloth based on the ply order, the warp shear angle, the weft shear angle, the normal unit vector, and the shear deformation reference unit vector; Determine the fiber cloth density, elastic stiffness matrix and material damping matrix of each layer of fiber cloth; Based on the ply sequence, the equivalent mechanical properties corresponding to the solid unit are determined in combination with the warp unit vector, weft unit vector, fiber cloth thickness and density, elastic stiffness matrix and material damping matrix to construct the numerical model.

[0101] Optionally, based on the ply sequence, the equivalent mechanical properties corresponding to the solid element are determined in combination with the warp unit vector, the weft unit vector, the fiber cloth thickness, the fiber cloth density, the elastic stiffness matrix, and the material damping matrix to construct the numerical model, including: Establish a local coordinate system corresponding to the number of plies; Obtaining a normal unit vector of the entity element, and using the first local coordinate in the unit direction of the normal vector as the local coordinate system of the entity element; Determine the transverse coordinate axis of the local coordinate system based on the warp direction unit vector and the weft direction unit vector, use the normal vector as the height coordinate axis of the local coordinate system, rotate the height coordinate axis according to the right-hand rule to obtain the longitudinal coordinate axis of the local coordinate system, and establish a transformation matrix between each ply and the solid element, wherein the transformation matrix is ​​used to represent the transformation relationship between stress or strain in the local coordinate system and the global coordinate system; Convert the elastic stiffness matrix and material damping matrix of each ply from the local coordinate system to the global coordinate system to obtain the stiffness matrix and damping matrix of each ply in the global coordinate system; The equivalent stiffness matrix of the solid element is calculated by weighted average based on the stress continuity condition, the strain continuity condition, the volume fraction of each ply, and the converted stiffness matrix. The volume fraction is determined based on the ply material thickness and the total thickness of the solid element. Based on the converted damping matrix, the equivalent stiffness matrix and the equivalent flexibility matrix are combined to calculate and obtain the equivalent damping matrix of the solid element; the equivalent flexibility matrix is ​​the inverse matrix of the equivalent stiffness matrix; Based on the fiber cloth thickness and fiber cloth density of each layer, the equivalent density of the solid unit is calculated by volume weighting; The equivalent stiffness matrix, equivalent damping matrix and equivalent density are used as equivalent mechanical characteristic parameters of the solid element to construct a numerical model.

[0102] Optionally, the implementation process of the numerical model includes: Perform a fishnet draping simulation and read the simulation results to determine the deformation distribution of the fiber cloth and the deformation measurement data of each draping; Using a programming interface of finite element analysis software to read the target structure grid data and ply material performance data; Using the shell-solid mesh mapping method, the entity elements of the target structure mesh are traversed to complete the mesh data mapping; The equivalent calculation formula of multi-layer structure is used to calculate the equivalent mechanical properties of the target structure grid; The programming interface function of the finite element analysis software is used to assign attributes to the target structure mesh.

[0103] It is understandable that the detailed functional implementation of each of the above units / modules can be found in the introduction of the aforementioned method embodiment, and will not be repeated here.

[0104] It should be understood that the above-mentioned device is used to execute the method in the above-mentioned embodiment. The implementation principle and technical effect of the corresponding program module in the device are similar to those described in the above-mentioned method. The working process of the device can refer to the corresponding process in the above-mentioned method and will not be repeated here.

[0105] Reference Figure 6 Based on the method in the above embodiment, an embodiment of the present application provides an electronic device, which may include: a processor (Processor) 610, a communication interface (Communications Interface) 620, a memory (Memory) 630 and a communication bus 640, wherein the processor 610, the communication interface 620, and the memory 630 communicate with each other via the communication bus 640. The processor 610 can call the logic instructions in the memory 630 to execute the method in the above embodiment.

[0106] In addition, the logic instructions in the aforementioned memory 630 can be implemented in the form of a software functional unit and, when sold or used as an independent product, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the portion that contributes to the prior art, or the portion of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application.

[0107] Based on the method in the above embodiment, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method in the above embodiment.

[0108] Based on the method in the above embodiment, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method in the above embodiment.

[0109] It is understood that the processor in the embodiments of the present application may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field programmable gate arrays (FPGA), other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. The general-purpose processor may be a microprocessor or any conventional processor.

[0110] The method steps in the embodiments of the present application can be implemented by hardware or by a processor executing software instructions. The software instructions can be composed of corresponding software modules, which can be stored in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), registers, hard disks, mobile hard disks, CD-ROMs, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be an integral part of the processor. The processor and storage medium can be located in an ASIC.

[0111] The above embodiments can be implemented in whole or in part using software, hardware, firmware, or any combination thereof. When implemented using software, they can be implemented in whole or in part in the form of a computer program product. The computer program product comprises one or more computer instructions. When loaded and executed on a computer, the computer program instructions fully or partially produce the processes or functions described in the embodiments of this application. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted via the computer-readable storage medium. The computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium accessible by a computer or a data storage device such as a server or data center that integrates one or more available media. The available medium can be magnetic media (e.g., floppy disk, hard disk, tape), optical media (e.g., DVD), or semiconductor media (e.g., solid-state drive (SSD)).

[0112] It will be understood that the various numerical numbers involved in the embodiments of the present application are merely distinctions for the convenience of description and are not intended to limit the scope of the embodiments of the present application.

[0113] It is easy for those skilled in the art to understand that the above is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application should be included in the scope of protection of the present application.

Claims

1. An equivalent mechanical modeling method for a variable thickness composite material component, characterized in that: include: Obtaining a draping deformation metric, and using the draping deformation metric to determine a fiber cloth deformation distribution using a fishnet method; The deformation distribution of the fiber cloth is mapped to the target structure grid using the grid mapping method; A numerical model of the equivalent mechanical characteristics of the target structure grid is constructed using a multi-layer structure equivalent calculation formula; the equivalent mechanical characteristics include equivalent density, equivalent stiffness matrix, and equivalent damping matrix; the parameters of the equivalent mechanical characteristics are obtained by combining ply information with coordinate system transformation and performing a weighted calculation of the stiffness, damping, and density of each ply; A numerical implementation program of the numerical model is determined, and the numerical model is passed through the numerical implementation program to perform equivalent mechanical analysis simulation on a variable thickness composite material component using a finite element method.

2. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 1, characterized in that: The obtaining of the draping deformation metric, determining the fiber cloth deformation distribution by using the draping deformation metric through a fishnet method, and mapping the fiber cloth deformation distribution to a target structure grid by using a grid mapping method includes: Using the fishing net method to perform draping simulation, determine the degree of deformation of the fiber cloth during the draping process, and obtain a deformation measure of the fiber cloth draping according to the deformation degree; Determine the draping simulation result based on the deformation metric, and determine the deformation distribution of the fiber cloth in the variable thickness composite material component; the fiber cloth deformation distribution includes the spatial shell unit grid and the fiber deformation metric distribution after the fiber cloth is draped; The shell-solid mesh mapping method is used to map the deformation distribution of the fiber cloth from the spatial shell element mesh to the target structural mesh of the variable thickness composite component for the hexahedral mesh used in mechanical simulation.

3. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 1, characterized in that: The draping deformation metric includes: warp yarn shear angle, weft yarn shear angle, normal unit vector, shear deformation reference unit vector and fiber cloth thickness.

4. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 2, characterized in that: The implementation process of the shell-solid mesh mapping method includes: Constructing an initialization matrix of size m×n, wherein the initialization matrix is ​​used to store the intersection relationship between plies and solid elements; m represents the number of solid elements in the target structure grid, and n represents the number of plies in the space shell grid; For each solid element in the target structural mesh, traverse the triangular mesh of the target ply of the spatial shell element mesh; Using a spatial triangle-spatial hexahedron collision detection algorithm to determine whether the solid unit and the target ply intersect, marking the empty matrix according to the determination result to obtain a marking matrix; For each entity cell in the target structure grid, searching for non-zero entries in the target row of the label matrix; According to the search results, the warp shear angle, weft shear angle, normal unit vector, shear deformation reference unit vector, and fiber cloth thickness of the target ply at the solid element are written into the solid element.

5. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 2, characterized in that: The numerical model for determining the equivalent mechanical properties of the target structure grid using the multi-layer structure equivalent calculation formula includes: Determining the ply order of the solid elements according to the normal element vectors of the solid elements of the target structural network; Determine the warp unit vector and the weft unit vector of each layer of fiber cloth based on the ply order, the warp shear angle, the weft shear angle, the normal unit vector, and the shear deformation reference unit vector; Determine the fiber cloth density, elastic stiffness matrix and material damping matrix of each layer of fiber cloth; Based on the ply sequence, the equivalent mechanical properties corresponding to the solid unit are determined in combination with the warp unit vector, weft unit vector, fiber cloth thickness and density, elastic stiffness matrix and material damping matrix to construct the numerical model.

6. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 5, characterized in that: Based on the ply sequence, the equivalent mechanical properties corresponding to the solid element are determined in combination with the warp unit vector, the weft unit vector, the fiber cloth thickness, the fiber cloth density, the elastic stiffness matrix, and the material damping matrix to construct the numerical model, including: Establish a local coordinate system corresponding to the number of plies; Obtaining a normal unit vector of the entity element, and using the first local coordinate in the unit direction of the normal vector as the local coordinate system of the entity element; Determine the transverse coordinate axis of the local coordinate system according to the warp direction unit vector and the weft direction unit vector, use the normal vector as the height coordinate axis of the local coordinate system, and rotate the height coordinate axis according to the right-hand rule to obtain the longitudinal coordinate axis of the local coordinate system; Establishing a transformation matrix between each ply and the solid element, wherein the transformation matrix is ​​used to represent the transformation relationship between stress or strain between the local coordinate system and the global coordinate system; Convert the elastic stiffness matrix and material damping matrix of each ply from the local coordinate system to the global coordinate system to obtain the stiffness matrix and damping matrix of each ply in the global coordinate system; The equivalent stiffness matrix of the solid element is calculated by weighted average based on the stress continuity condition, the strain continuity condition, the volume fraction of each ply, and the converted stiffness matrix. The volume fraction is determined based on the ply material thickness and the total thickness of the solid element. Based on the converted damping matrix, the equivalent stiffness matrix and the equivalent flexibility matrix are combined to calculate and obtain the equivalent damping matrix of the solid element; the equivalent flexibility matrix is ​​the inverse matrix of the equivalent stiffness matrix; Based on the fiber cloth thickness and fiber cloth density of each layer, the equivalent density of the solid unit is calculated by volume weighting; The equivalent stiffness matrix, equivalent damping matrix and equivalent density are used as equivalent mechanical characteristic parameters of the solid element to construct a numerical model.

7. The equivalent mechanical modeling method for a variable thickness composite material component according to claim 1, characterized in that: The implementation process of the numerical model includes: Perform a fishnet draping simulation and read the simulation results to determine the deformation distribution of the fiber cloth and the deformation measurement data of each draping; Using a programming interface of finite element analysis software to read the target structure grid data and ply material performance data; Using the shell-solid mesh mapping method, the entity elements of the target structure mesh are traversed to complete the mesh data mapping; The equivalent calculation formula of multi-layer structure is used to calculate the equivalent mechanical properties of the target structure grid; The programming interface function of the finite element analysis software is used to assign attributes to the target structure mesh.

8. An equivalent mechanical modeling system for variable thickness composite material components, characterized in that: include: a deformation distribution determination module for obtaining a draping deformation metric and determining a fiber cloth deformation distribution using the draping deformation metric through a fishnet method; A mapping module, used for mapping the deformation distribution of the fiber cloth to the target structure grid using a grid mapping method; A model construction module is used to construct a numerical model of the equivalent mechanical characteristics of the target structure grid using an equivalent calculation formula for a multi-layer structure; the equivalent mechanical characteristics include an equivalent density, an equivalent stiffness matrix, and an equivalent damping matrix; the equivalent calculation formula for a multi-layer structure is constructed based on the parameters of the equivalent mechanical characteristics, which are obtained by weighted calculation of the stiffness, damping, and density of each ply by combining ply information with coordinate system transformation. The simulation module is used to determine the numerical implementation program of the numerical model, and use the finite element method to perform equivalent mechanical analysis simulation on the variable thickness composite material component through the numerical model.

9. An electronic device, characterized in that: include: at least one memory for storing a computer program; At least one processor is used to execute the program stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed on a processor, the processor is caused to execute the method according to any one of claims 1 to 7.

Images