Plate spring thickness determination method and device for multi-body dynamic model of plate spring

Abstract

The invention provides a leaf spring thickness determination method and device for a leaf spring multi-body dynamics model, and relates to the technical field of vehicle dynamics simulation. The method comprises the steps that geometric parameters and the number of discrete elements of the steel plate spring are obtained, and based on the geometric parameters and the number of discrete elements, the thickness value of each discrete element matched with the target mechanical property is determined through a publicly verifiable deterministic algorithm so as to form discrete element thickness distribution; and outputting the discrete unit thickness distribution as an input file of the steel plate spring multi-body dynamic model. According to the method, the publicly verifiable deterministic algorithm is introduced, and the thickness distribution matched with the target mechanical property is determined based on the geometrical parameters and the discrete element number of the steel plate spring, so that the problems of low model transparency and difficulty in ensuring the precision caused by dependence on artificial experience or a black box algorithm in the prior art are effectively improved; and the credibility and the prediction accuracy of the simulation model are obviously improved.

Description

Technical Field

[0001] This application relates to the field of vehicle dynamics simulation technology, and in particular to a method and apparatus for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring. Background Technology

[0002] As the core elastic and guiding element of commercial vehicle suspension systems, leaf springs directly affect the vehicle's handling stability, ride comfort, and reliability. In modern vehicle forward development processes, virtual simulation based on Multi-Body Dynamics (MBD) has become an indispensable part of predicting and optimizing suspension performance. To achieve efficient simulation, the continuous leaf spring is typically discretized into a series of multi-body models connected by beam elements. Among these challenges, assigning appropriate thickness parameters to each discrete element to ensure its lumped mechanical properties closely match the overall mechanical properties of the real leaf spring is fundamental to building a high-fidelity simulation model and has long been a key technical difficulty for engineers.

[0003] In related technologies, determining the leaf spring thickness in multibody dynamics models of leaf springs mainly relies on two existing methods. The first is the user input method, where engineers manually assign thickness values ​​to each discrete unit in the model based on design drawings and personal experience. The second is the software-internal calculation method, represented by mainstream commercial software such as Adams. Users only need to provide macroscopic parameters such as initial thickness and number of leaves, and the software's internal algorithm can automatically complete the thickness allocation. However, the former relies heavily on subjective experience, resulting in low modeling efficiency and difficulty in guaranteeing accuracy; the latter operates entirely in a "black box" mode, with its specific calculation formulas and algorithms not disclosed. Users cannot understand or verify the logic of the thickness distribution, and it is difficult to modify and optimize it according to actual needs.

[0004] Therefore, providing a transparent and controllable method for determining the thickness of leaf springs has become an urgent need to improve the accuracy and simulation reliability of multibody dynamics models of leaf springs. Summary of the Invention

[0005] This application provides a method and apparatus for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, in order to improve the problems of low model transparency and difficulty in guaranteeing accuracy caused by the reliance on manual experience or black-box algorithms in related technologies.

[0006] In a first aspect, this application provides a method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, including:

[0007] Obtain the geometric parameters and number of discrete elements of the leaf spring;

[0008] Based on geometric parameters and the number of discrete elements, a publicly verifiable deterministic algorithm is used to determine the thickness values ​​of each discrete element that match the target mechanical properties, thereby forming a discrete element thickness distribution.

[0009] The discrete element thickness distribution is output as the input file for the multibody dynamics model of the leaf spring.

[0010] In one possible implementation, the geometric parameters include the contour point coordinates and initial thickness of the leaf spring. Based on the geometric parameters and the number of discrete elements, the thickness value of each discrete element matching the target mechanical properties is determined by a publicly verifiable deterministic algorithm. This includes: performing curve fitting based on the contour point coordinates to obtain the working surface geometry of the leaf spring; and determining the thickness value of each discrete element in the normal direction at its corresponding position by coordinate transformation based on the working surface geometry and the initial thickness.

[0011] In one possible implementation, curve fitting is performed based on the coordinates of the contour points to obtain the working surface geometry of the leaf spring, including: performing a fourth-order polynomial fitting on the coordinates of the contour points using the least squares method to generate the contour curve of the working surface of the leaf spring.

[0012] In one possible implementation, based on the working surface geometry and initial thickness, the thickness value of each discrete element in the normal direction at its corresponding position is determined by coordinate transformation, including: determining the paired surface geometry corresponding to the working surface geometry based on the working surface geometry and initial thickness; and determining the thickness value of each discrete element in the normal direction based on the working surface geometry and the paired surface geometry.

[0013] In one possible implementation, the geometry of the paired surface corresponding to the working surface geometry is determined based on the working surface geometry and the initial thickness, including: determining the tangent angle at each discrete point according to the contour curve of the working surface geometry; determining the contour point coordinates of the paired surface based on the initial thickness and the tangent angle; and performing a fourth-order polynomial fitting on the contour point coordinates of the paired surface to generate a contour curve characterizing the geometry of the paired surface.

[0014] In one possible implementation, the coordinates of the contour points of the paired surfaces satisfy the following formula:

[0015]

[0016]

[0017] in, Let z be the z-coordinate of the i-th discrete point on the contour curve of the working surface. Let z be the z-coordinate of the corresponding point on the profile curve of the paired surface. For the initial thickness, Let be the tangent angle at the i-th discrete point; denoted as the vertical distance between the working surface geometry and the paired surface geometry at the i-th discrete point.

[0018] In one possible implementation, the thickness value of each discrete element in the normal direction is determined based on the working surface geometry and the paired surface geometry, including: determining the fitted vertical distance corresponding to each discrete point based on the contour curve of the working surface geometry and the contour curve of the paired surface geometry; and determining the thickness value of each discrete element in the normal direction through coordinate rotation transformation based on the fitted vertical distance corresponding to each discrete point and the tangent angle of the working surface geometry at each discrete point.

[0019] In one possible implementation, the thickness value of each discrete element in the normal direction is determined by coordinate rotation transformation, including: constructing an initial vector based on the fitted vertical distance corresponding to each discrete point; rotating the initial vector around the origin of the coordinate system clockwise according to the tangent angle of the working surface geometry at each discrete point to obtain the rotated vector; extracting the component of the rotated vector in the normal direction, and using the absolute value of the component as the thickness value of the corresponding discrete element in the normal direction.

[0020] Secondly, this application provides a leaf spring thickness determination device for a multibody dynamics model of a leaf spring, comprising:

[0021] The acquisition module is used to acquire the geometric parameters and the number of discrete elements of the leaf spring;

[0022] The determination module is used to determine the thickness value of each discrete element that matches the target mechanical properties based on geometric parameters and the number of discrete elements, using a publicly verifiable deterministic algorithm, so as to form a discrete element thickness distribution.

[0023] The output module is used to output the thickness distribution of discrete elements as the input file for the multibody dynamics model of the leaf spring.

[0024] Thirdly, this application provides an electronic device, including: a processor, and a memory communicatively connected to the processor;

[0025] Memory is used to store instructions executed by the computer;

[0026] A processor for executing computer-executable instructions stored in memory to implement any of the methods of the first aspect.

[0027] Fourthly, this application provides a computer-readable storage medium storing computer-executable instructions, which, when executed, are used to implement the method of any one of the first aspects.

[0028] Fifthly, this application provides a computer program product, including a computer program that, when executed, implements the method of any one of the first aspects.

[0029] The method and apparatus for determining the thickness of a leaf spring in a multibody dynamics model provided in this application obtain the geometric parameters and the number of discrete elements of the leaf spring. Based on the geometric parameters and the number of discrete elements, a publicly verifiable deterministic algorithm is used to determine the thickness value of each discrete element that matches the target mechanical properties, so as to form a discrete element thickness distribution. Furthermore, the discrete element thickness distribution is output as the input file of the multibody dynamics model of the leaf spring.

[0030] In this process, by introducing a publicly verifiable deterministic algorithm, a thickness distribution matching the target mechanical properties is determined based on the geometric parameters and the number of discrete elements of the leaf spring. This effectively improves the problems of low model transparency and difficulty in guaranteeing accuracy caused by reliance on manual experience or black-box algorithms in existing technologies. Specifically, the physical rationality of the thickness distribution is ensured through explicit mathematical principles, making the mechanical properties of the discrete beam model highly consistent with the real leaf spring, significantly improving the credibility and prediction accuracy of the simulation model. At the same time, the mechanism for automatically determining the thickness and outputting standardized documents achieves seamless integration with the multibody dynamics modeling environment, significantly reducing manual operation steps. While ensuring model quality, it significantly improves modeling efficiency and result consistency, thus providing a reliable technical guarantee for the refined design and performance optimization of vehicle suspension systems. Attached Figure Description

[0031] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.

[0032] Figure 1 A schematic diagram illustrating an application scenario of the leaf spring thickness determination method for a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application.

[0033] Figure 2 A flowchart illustrating a method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application;

[0034] Figure 3 A schematic diagram illustrating the positional relationship between the contour curves of the working surface and the contour curves of the mating surface, provided for an exemplary embodiment of this application;

[0035] Figure 4 Another schematic flowchart of a method for determining the thickness of a leaf spring for a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application;

[0036] Figure 5 A schematic diagram of a leaf spring thickness determination device for a multibody dynamics model of a leaf spring provided as an exemplary embodiment of this application;

[0037] Figure 6 A schematic diagram of the structure of an electronic device provided as an exemplary embodiment of this application.

[0038] The accompanying drawings illustrate specific embodiments of this application, which will be described in more detail below. These drawings and descriptions are not intended to limit the scope of the concept in any way, but rather to illustrate the concept of this application to those skilled in the art through reference to particular embodiments. Detailed Implementation

[0039] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0040] The terms “first,” “second,” etc., used in the specification and claims of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented, for example, in orders other than those illustrated or described herein. Furthermore, the terms “comprising” and “having,” and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, products, or apparatus.

[0041] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties. Furthermore, the collection, use and processing of the relevant data must comply with relevant laws, regulations and standards, and corresponding operation entry points are provided for users to choose to authorize or refuse.

[0042] In related technologies, existing multibody dynamics modeling methods for leaf springs have significant drawbacks in determining the thickness of discrete elements. Firstly, manual input methods relying on human experience are not only inefficient, but the rationality of thickness allocation also depends on the engineer's personal experience, leading to unstable model accuracy and poor reproducibility. Secondly, while using thickness allocation algorithms built into commercial software can achieve preliminary automation, their computational logic is completely closed. Users cannot understand the underlying mechanisms of thickness distribution, nor can they adjust and optimize them according to specific simulation needs. This "black box" operation mode makes the thickness allocation process opaque and uncontrollable, causing systematic deviations between the key mechanical parameters such as stiffness characteristics and stress distribution of the discrete beam model and the actual leaf spring. This severely weakens the credibility of simulation results, thereby restricting the design and optimization efficiency of simulation-based suspension systems.

[0043] To address the aforementioned issues, this application proposes a scheme for determining the leaf spring thickness in a multibody dynamics model of a leaf spring. By constructing a deterministic algorithm based on publicly available mathematical principles, the algorithm takes the geometric parameters of the leaf spring and the number of discrete elements as input, automatically generating a discrete element thickness distribution that matches the target mechanical properties. Specifically, verifiable mathematical methods such as curve fitting and coordinate transformation are used to accurately calculate the thickness value of each discrete element in the normal direction, ensuring the physical rationality of the thickness distribution and the accuracy of the model's mechanical properties. Then, through a standardized output interface, the thickness distribution data is directly converted into the input file for the multibody dynamics model of the leaf spring, achieving full automation from parameter input to model establishment. This effectively improves the limitations of traditional methods relying on manual experience or black-box algorithms, significantly enhancing the transparency, efficiency, and reproducibility of the modeling process while ensuring model accuracy.

[0044] Figure 1 This is a schematic diagram illustrating an application scenario of the leaf spring thickness determination method for a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application. For example... Figure 1 As shown, this application scenario includes a client 11 and a server 12, wherein the number of clients 11 can be at least one. In practical applications, when the server 12 detects the leaf spring thickness determination command submitted by the client 11, it executes the leaf spring thickness determination method for the multibody dynamics model of a leaf spring provided in this application, and outputs the discrete element thickness distribution as the input file of the multibody dynamics model of the leaf spring.

[0045] It should be noted that server 12 can also be replaced by a server cluster or other computing devices with a certain computing power. Both the first and second clients can be mobile phones, computers, laptops, or personal digital assistants (PDAs).

[0046] The following is combined with Figure 1 Application scenarios, refer to Figure 2 This application describes a method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, according to an exemplary embodiment of this application. It should be noted that the above application scenarios are shown only to facilitate understanding of the spirit and principles of this application, and the embodiments of this application are not limited to those described herein. Figure 1 The limitations of the application scenarios shown are explained. It should be noted that the following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments. The embodiments of this application will now be described with reference to the accompanying drawings.

[0047] Figure 2 This is a flowchart illustrating a method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application. Figure 2 As shown, the method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring includes the following steps:

[0048] S201. Obtain the geometric parameters and number of discrete elements of the leaf spring.

[0049] Among them, geometric parameters refer to a set of basic data describing the physical structure, shape, and size of the leaf spring; discrete elements refer to the basic mechanical elements into which a single leaf spring is divided along its length in the process of multibody dynamics modeling to transform a continuous elastic body into a computable discrete model. Each discrete element will be modeled as a beam element with independent cross-sectional properties in the simulation model; the number of discrete elements refers to the total number of discrete elements into which a single leaf spring is divided. This parameter is the core variable controlling the calculation accuracy and solution efficiency of the model. The more elements there are, the stronger the model's ability to characterize complex deformations, but the calculation cost also increases accordingly.

[0050] For example, in response to receiving a leaf spring thickness determination command, the geometric parameters and number of discrete elements of the leaf spring are read from a preset data file fore_1_part.dat.

[0051] S202. Based on geometric parameters and the number of discrete elements, the thickness values ​​of each discrete element that match the target mechanical properties are determined by a publicly verifiable deterministic algorithm to form a discrete element thickness distribution.

[0052] In some embodiments, the geometric parameters include the contour point coordinates and initial thickness of the leaf spring. Based on the geometric parameters and the number of discrete elements, the thickness value of each discrete element that matches the target mechanical properties is determined by a publicly verifiable deterministic algorithm. This includes: performing curve fitting based on the contour point coordinates to obtain the working surface geometry of the leaf spring; and determining the thickness value of each discrete element in the normal direction at its corresponding position by coordinate transformation based on the working surface geometry and the initial thickness.

[0053] The geometric parameters include the number of discrete points, leaf spring width, offset, and installation length. Target mechanical properties refer to the key mechanical characteristics that the leaf spring must achieve in its design, providing targets and constraints for the thickness distribution calculation. Specifically, target mechanical properties include, but are not limited to, the overall stiffness curve of the leaf spring, stress distribution under specific loads, rated load-bearing capacity, or specific deformation characteristics (such as arc height). A publicly verifiable deterministic algorithm is used: its core lies in converting the geometric parameters of the leaf spring into the thickness distribution of discrete elements through a series of explicit mathematical steps. This algorithm is completely public, logically deterministic, and does not rely on any undisclosed internal logic or empirical parameters. The following section uses the first half of the first leaf spring as an example to illustrate the specific implementation steps of this algorithm in detail.

[0054] For example, the geometric parameters and number of discrete elements of the leaf spring are read from the preset data file fore_1_part.dat, such as the number of discrete points being 11, the leaf spring width being 6.0000E+01, the initial thickness being 1.0000E+01, the offset being 0.0000E+00, the installation length being 1.1000E+02, and the number of discrete elements being 10; and the contour point coordinates of the leaf spring are read from the preset data file fore_1.dat, for example, Table 1 shows an example of the contour point coordinates of the leaf spring provided in the exemplary embodiment of this application.

[0055]

[0056] For example, curve fitting is performed based on the coordinates of the aforementioned contour points to obtain a mathematical model of the working surface representing the geometry of the leaf spring. Then, combined with the initial thickness parameters, the geometric data is converted into thickness values ​​in the normal direction of each discrete element at its corresponding position through coordinate transformation. Specifically, the algorithm constructs a geometric model of the working surface, calculates the local geometric features at each discrete point, derives the thickness distribution curve based on the initial thickness and geometric relationship, and finally obtains the normal thickness value of each discrete element through coordinate rotation transformation. The entire process ensures the physical rationality of the thickness distribution through rigorous mathematical derivation, enabling the final discrete element thickness distribution to accurately reflect the actual mechanical properties of the leaf spring and maintain a high degree of consistency with the target mechanical performance requirements. This deterministic algorithm-based calculation method guarantees the reliability and reproducibility of the results, providing accurate input data for subsequent multibody dynamics modeling.

[0057] By following the above steps, the thickness values ​​of all discrete elements that match the target mechanical properties can be determined, thereby forming a thickness distribution. The entire process is based on publicly available mathematical principles, with clear steps and definite results. Anyone skilled in the art can reproduce the implementation process based on this disclosure.

[0058] S203. Output the thickness distribution of discrete elements as the input file for the multibody dynamics model of the leaf spring.

[0059] For example, the normal thickness values ​​of the 10 discrete elements output in step S202 are sequentially packaged into a text file conforming to a specific format. This text file serves as direct input to the simulation model and can be parsed by multibody dynamics software such as Adams. It automatically completes the parameterized definition of the cross-sectional properties of each element in the leaf spring discrete beam model, ultimately driving the automatic generation of a high-precision model.

[0060] It should be noted that this method is universal. The thickness determination process described above, taking the first half of the first leaf spring as an example, is also applicable to the second half of the first leaf spring, as well as all leaf springs from the second to the nth leaf spring. That is, it is not only applicable to single leaf springs, but can also be extended to multi-leaf springs, gradually thickened springs and other composite structures. For these structures, the system will reuse the exact same deterministic algorithm. The only difference is that specific geometric parameters of the corresponding leaf spring or segment need to be input.

[0061] The leaf spring thickness determination method for multibody dynamics models of leaf springs provided in this application introduces a publicly verifiable deterministic algorithm. Based on the geometric parameters and number of discrete elements of the leaf spring, it determines a thickness distribution that matches the target mechanical properties. This effectively improves the problems of low model transparency and difficulty in guaranteeing accuracy caused by reliance on manual experience or black-box algorithms in the prior art. Specifically, it ensures the physical rationality of the thickness distribution through explicit mathematical principles, making the mechanical properties of the discrete beam model highly consistent with the real leaf spring, significantly improving the credibility and prediction accuracy of the simulation model. At the same time, the mechanism of automated thickness determination and standardized document output achieves seamless integration with the multibody dynamics modeling environment, significantly reducing manual operation steps. While ensuring model quality, it significantly improves modeling efficiency and result consistency, thus providing a reliable technical guarantee for the refined design and performance optimization of vehicle suspension systems.

[0062] In some embodiments, curve fitting is performed based on the coordinates of the contour points to obtain the working surface geometry of the leaf spring, including: performing fourth-order polynomial fitting on the coordinates of the contour points using the least squares method to generate the contour curve of the working surface of the leaf spring.

[0063] Among them, the least squares method is a mathematical optimization technique that finds the best function match for the data by minimizing the sum of squares of the errors. In this embodiment, it is used to ensure that the overall deviation between the fitted curve and the original contour points is minimized. The fourth-order polynomial fitting refers to the method of using a fourth-order polynomial function to fit the discrete data to the curve. This degree selection can ensure the fitting accuracy for the complex bending shape of the leaf spring, while avoiding overfitting that may be caused by too high-order polynomials. The contour curve of the working surface refers to a continuous mathematical function obtained by curve fitting that can accurately describe the geometry of the neutral surface or upper surface of the leaf spring. This curve transforms the discrete contour point data into a continuous geometric expression, providing an accurate mathematical basis for subsequent thickness calculation.

[0064] For example, by performing a fourth-order polynomial fitting on the x and z coordinates of a series of points using the least squares method, the contour curve of the working surface of the leaf spring is generated. Satisfy the following formula:

[0065]

[0066] in, , , , and The coefficients obtained by fitting the curve using the least squares method collectively determine the specific shape of the quartic curve, thus accurately matching the true contour of the leaf spring.

[0067] In this embodiment, the working surface geometry of the leaf spring is constructed by performing a fourth-order polynomial fitting on the contour point coordinates using the least squares method. Specifically, the contour curve that best represents the true geometric characteristics of the leaf spring is automatically generated through mathematical optimization, effectively eliminating the subjectivity of manual judgment and random errors in the original data. Furthermore, the fourth-order polynomial model can accurately capture the complex bending characteristics of the leaf spring (including asymmetry and local curvature changes) while reducing the risk of overfitting caused by excessively high orders, ensuring the accuracy and stability of the geometric model. The resulting high-precision working surface provides a reliable geometric benchmark for subsequent thickness value determination, ensuring the physical rationality of the thickness distribution of discrete units from the source. This results in the established multibody dynamics model of the leaf spring having higher simulation accuracy and prediction reliability.

[0068] In some embodiments, based on the working surface geometry and initial thickness, the thickness value of each discrete element in the normal direction at its corresponding position is determined by coordinate transformation, including: determining the paired surface geometry corresponding to the working surface geometry based on the working surface geometry and initial thickness; and determining the thickness value of each discrete element in the normal direction based on the working surface geometry and the paired surface geometry.

[0069] For example, based on the known working surface geometry (i.e., the mathematical model of the upper surface profile of the leaf spring) and the given initial thickness parameters, the initial spatial position data of the mating surfaces are derived through geometric relationships; the initial spatial position data are then reconstructed by curve fitting to generate a smooth and continuous mating surface geometry.

[0070] In some embodiments, the geometry of the paired surface corresponding to the working surface geometry is determined based on the working surface geometry and the initial thickness, including: determining the tangent angle at each discrete point according to the contour curve of the working surface geometry; determining the contour point coordinates of the paired surface based on the initial thickness and the tangent angle; and performing a fourth-order polynomial fitting on the contour point coordinates of the paired surface to generate a contour curve characterizing the geometry of the paired surface.

[0071] Among them, the tangent angle refers to the angle between the tangent line of the working surface profile curve of the leaf spring at a certain point and the horizontal direction, reflecting the local inclination of that point; the mating surface refers to the corresponding surface that maintains a constant normal distance from the working surface, which in this embodiment refers to the lower surface profile of the leaf spring.

[0072] In some embodiments, the coordinates of the contour points of the paired surfaces satisfy the following formula:

[0073]

[0074]

[0075] in, Let z be the z-coordinate of the i-th discrete point on the contour curve of the working surface. Let z be the z-coordinate of the corresponding point on the profile curve of the paired surface. For the initial thickness, Let be the tangent angle at the i-th discrete point; Let be the vertical distance between the working surface geometry and the paired surface geometry at the i-th discrete point. For example, Figure 3 A schematic diagram illustrating the positional relationship between the contour curves of the working surface and the contour curves of the mating surface, provided for an exemplary embodiment of this application.

[0076] Accordingly, the coordinates of the contour points of the paired surface are determined using the above formula; further, the x-coordinates of a series of points are determined using the least squares method. The coordinates are fitted with a fourth-order polynomial to generate profile curves representing the geometry of the paired surfaces. Satisfy the following formula:

[0077]

[0078] in, , , , and The coefficients obtained by fitting the curve using the least squares method collectively determine the specific shape of the quartic curve, thus accurately matching the true contour of the leaf spring.

[0079] In some embodiments, the thickness value of each discrete element in the normal direction is determined based on the working surface geometry and the paired surface geometry, including: determining the fitted vertical distance corresponding to each discrete point based on the contour curve of the working surface geometry and the contour curve of the paired surface geometry; and determining the thickness value of each discrete element in the normal direction through coordinate rotation transformation based on the fitted vertical distance corresponding to each discrete point and the tangent angle of the working surface geometry at each discrete point.

[0080] For example, for each discrete point, by fitting the profile curve based on the working surface geometry and the profile curve of the paired surface geometry, the fitted vertical distance corresponding to each discrete point is obtained. For example, the fitted vertical distance corresponding to each discrete point satisfies the following formula:

[0081]

[0082] in, The vertical distance corresponding to each discrete point after fitting is represented.

[0083] In some embodiments, the thickness value of each discrete unit in the normal direction is determined by coordinate rotation transformation, including: constructing an initial vector based on the fitted vertical distance corresponding to each discrete point; rotating the initial vector around the origin of the coordinate system clockwise according to the tangent angle of the working surface geometry at each discrete point to obtain the rotated vector; extracting the component of the rotated vector in the normal direction, and using the absolute value of the component as the thickness value of the corresponding discrete unit in the normal direction.

[0084] For example, based on the fitted vertical distances corresponding to each discrete point. Construct initial vectors This vector represents the vertical distance from the working surface to the paired surface in the global coordinate system.

[0085] Accordingly, the initial vector is rotated clockwise around the origin according to the tangent angle θ of the working surface geometry at each discrete point. For a point (x, y) on a two-dimensional plane, the formula for calculating its new coordinates (x', y') is as follows:

[0086]

[0087]

[0088] Represented in matrix form, it is:

[0089]

[0090] Correspondingly, the coordinates of the rotated vector: , Furthermore, the components of the rotated vector in the normal direction are extracted. The absolute value of this component This is used as the thickness value of the corresponding discrete element in the normal direction. The vertical distance is then transformed through coordinate rotation. Converted to the true thickness along the normal direction, when the leaf spring has curvature (θ≠0), the perpendicular distance... The thickness is always greater than the normal thickness, and is precisely corrected by multiplying by the cosθ factor; when the leaf spring is straight (θ=0), cosθ=1, and the thickness value is equal to the vertical distance, which matches the actual physical situation.

[0091] In this embodiment, a complete three-dimensional geometric representation of the leaf spring is established by accurately constructing corresponding paired surfaces based on the working surface geometry and initial thickness, ensuring the accuracy of thickness calculation from the outset. Subsequently, by analyzing the spatial positional relationship between the two surfaces at each discrete element, the true thickness value in the normal direction is determined, effectively eliminating the systematic error caused by the traditional method of mistakenly taking the vertical distance as thickness due to neglecting the surface curvature. This thickness determination method based on rigorous geometric principles ensures the physical accuracy of the discrete element cross-sectional properties, enabling the established leaf spring multibody dynamics model to accurately reflect the stiffness characteristics and mechanical behavior of the actual leaf spring. This significantly improves the reliability and prediction accuracy of simulation results, providing reliable technical support for the precise design of vehicle suspension systems.

[0092] Figure 4 Another schematic flowchart illustrating a method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application. Figure 4 As shown, the method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring includes the following steps:

[0093] S401. Obtain the geometric parameters and number of discrete elements of the leaf spring. The geometric parameters include the coordinates of the leaf spring's contour points and the initial thickness.

[0094] S402. The contour curve of the working surface of the leaf spring is generated by performing a fourth-order polynomial fitting on the coordinates of the contour points using the least squares method.

[0095] For example, by performing a fourth-order polynomial fitting on the x and z coordinates of a series of points using the least squares method, the contour curve of the working surface of the leaf spring is generated. Satisfy the following formula:

[0096]

[0097] in, , , , and The coefficients obtained by fitting the curve using the least squares method collectively determine the specific shape of the quartic curve, thus accurately matching the true contour of the leaf spring.

[0098] S403. Determine the tangent angle at each discrete point based on the contour curve of the working surface geometry.

[0099] S404. Based on the initial thickness and tangent angle, determine the coordinates of the contour points of the paired surfaces.

[0100] For example, the coordinates of the contour points of the paired surface satisfy the following formula:

[0101]

[0102]

[0103] in, Let z be the z-coordinate of the i-th discrete point on the contour curve of the working surface. Let z be the z-coordinate of the corresponding point on the profile curve of the paired surface. For the initial thickness, Let be the tangent angle at the i-th discrete point; denoted as the vertical distance between the working surface geometry and the paired surface geometry at the i-th discrete point.

[0104] S405. Perform a fourth-order polynomial fitting on the coordinates of the contour points of the paired surfaces to generate a contour curve that characterizes the geometry of the paired surfaces.

[0105] For example, the least squares method is used to determine the x-coordinates of a series of points. The coordinates are fitted with a fourth-order polynomial to generate profile curves representing the geometry of the paired surfaces. Satisfy the following formula:

[0106]

[0107] in, , , , and The coefficients obtained by fitting the curve using the least squares method collectively determine the specific shape of the quartic curve, thus accurately matching the true contour of the leaf spring.

[0108] S406. Based on the contour curve of the working surface geometry and the contour curve of the paired surface geometry, determine the fitted vertical distance of each discrete point.

[0109] For example, the fitted vertical distances for each discrete point satisfy the following formula:

[0110]

[0111] in, The vertical distance corresponding to each discrete point after fitting is represented.

[0112] S407. Construct an initial vector based on the fitted vertical distances corresponding to each discrete point.

[0113] For example, based on the fitted vertical distances corresponding to each discrete point. Construct initial vectors .

[0114] S408. Rotate the initial vector clockwise around the origin of the coordinate system according to the tangent angle of the working surface geometry at each discrete point to obtain the rotated vector.

[0115] S409. Extract the components of the rotated vector in the normal direction, and use the absolute value of the components as the thickness value of the corresponding discrete element in the normal direction to form the discrete element thickness distribution.

[0116] For example, the coordinates of the rotated vector: , Furthermore, the components of the rotated vector in the normal direction are extracted. The absolute value of this component This is the thickness value of the corresponding discrete element in the normal direction.

[0117] S410. Output the thickness distribution of discrete elements as the input file for the multibody dynamics model of the leaf spring.

[0118] In summary, this application has at least the following advantages:

[0119] I. By introducing a publicly verifiable deterministic algorithm, a thickness distribution matching the target mechanical properties is determined based on the geometric parameters and number of discrete elements of the leaf spring. This effectively improves the problems of low model transparency and difficulty in guaranteeing accuracy caused by reliance on manual experience or black-box algorithms in existing technologies. Specifically, the physical rationality of the thickness distribution is ensured through explicit mathematical principles, making the mechanical properties of the discrete beam model highly consistent with the real leaf spring, significantly improving the credibility and prediction accuracy of the simulation model. At the same time, the mechanism of automatically determining the thickness and outputting standardized documents achieves seamless integration with the multibody dynamics modeling environment, significantly reducing manual operation steps. While ensuring model quality, it significantly improves modeling efficiency and result consistency, thus providing a reliable technical guarantee for the refined design and performance optimization of vehicle suspension systems.

[0120] Second, the working surface geometry of the leaf spring is constructed by using the least squares method to fit the coordinates of the contour points with a fourth-order polynomial. Specifically, the contour curve that best represents the true geometric characteristics of the leaf spring is automatically generated through mathematical optimization, effectively eliminating the subjectivity of manual judgment and random errors in the original data. Moreover, the fourth-order polynomial model can accurately capture the complex bending characteristics of the leaf spring (including asymmetry and local curvature changes) and reduce the risk of overfitting caused by excessively high orders, ensuring the accuracy and stability of the geometric model. The resulting high-precision working surface provides a reliable geometric benchmark for subsequent thickness value determination, ensuring the physical rationality of the thickness distribution of discrete units from the source. This results in the established multibody dynamics model of the leaf spring having higher simulation accuracy and prediction reliability.

[0121] Third, by accurately constructing corresponding paired surfaces based on the working surface geometry and initial thickness, a complete three-dimensional geometric representation of the leaf spring was established, ensuring the accuracy of thickness calculation from the outset. Subsequently, by analyzing the spatial positional relationship between the two surfaces at each discrete element, the true thickness value in the normal direction was determined, effectively eliminating the systematic error caused by the traditional method of mistakenly taking the vertical distance as thickness due to neglecting the surface curvature. This thickness determination method based on rigorous geometric principles ensures the physical accuracy of the discrete element cross-sectional properties, enabling the established leaf spring multibody dynamics model to accurately reflect the stiffness characteristics and mechanical behavior of the actual leaf spring. This significantly improves the reliability and prediction accuracy of the simulation results, providing reliable technical support for the precise design of vehicle suspension systems.

[0122] IV. The thickness determination process provided in this application is not only applicable to the example single-leaf spring structure, but can also be seamlessly extended to multi-leaf springs, gradually thickening springs, and other composite leaf spring structures, demonstrating excellent versatility. Its core advantage lies in the fact that for different types of leaf springs, the system does not need to change the algorithm kernel; it only needs to input the specific geometric parameters of the corresponding leaf spring or segment to reuse the same publicly verified deterministic algorithm to automatically generate a high-precision thickness distribution. This "unified algorithm, parameter-driven" mode greatly improves the engineering applicability and modeling efficiency of the method, effectively reduces the technical threshold and development cost of simulating different types of leaf springs, and provides strong common technical support for the unified, standardized design and performance optimization of complex suspension systems.

[0123] The following are embodiments of the apparatus described in this application, which can be used to execute the embodiments of the method described in this application. For details not disclosed in the apparatus embodiments of this application, please refer to the embodiments of the method described in this application.

[0124] Figure 5 This is a schematic diagram of a leaf spring thickness determination device for a multibody dynamics model of a leaf spring, provided as an exemplary embodiment of this application. Figure 5 As shown, the leaf spring thickness determination device 50 for the multibody dynamics model of a leaf spring includes an acquisition module 51, a determination module 52, and an output module 53, wherein:

[0125] The acquisition module 51 is used to acquire the geometric parameters and the number of discrete units of the leaf spring;

[0126] The determination module 52 is used to determine the thickness value of each discrete element that matches the target mechanical properties based on geometric parameters and the number of discrete elements, through a publicly verifiable deterministic algorithm, so as to form a discrete element thickness distribution.

[0127] Output module 53 is used to output the thickness distribution of discrete elements as the input file for the multibody dynamics model of the leaf spring.

[0128] In one possible implementation, the geometric parameters include the contour point coordinates and initial thickness of the leaf spring. Based on the geometric parameters and the number of discrete units, the determining module 52 can be specifically used to: perform curve fitting based on the contour point coordinates to obtain the working surface geometry of the leaf spring; and determine the thickness value of each discrete unit in the normal direction at its corresponding position through coordinate transformation based on the working surface geometry and the initial thickness.

[0129] In one possible implementation, the determining module 52 can also be used to: perform a fourth-order polynomial fitting on the coordinates of the contour points using the least squares method to generate the contour curve of the working surface of the leaf spring.

[0130] In one possible implementation, the determining module 52 can also be used to: determine the paired surface geometry corresponding to the working surface geometry based on the working surface geometry and the initial thickness; and determine the thickness value of each discrete element in the normal direction based on the working surface geometry and the paired surface geometry.

[0131] In one possible implementation, the determining module 52 can also be used to: determine the tangent angle at each discrete point based on the contour curve of the working surface geometry; determine the contour point coordinates of the paired surface based on the initial thickness and the tangent angle; and perform a fourth-order polynomial fitting on the contour point coordinates of the paired surface to generate a contour curve characterizing the geometry of the paired surface.

[0132] In one possible implementation, the coordinates of the contour points of the paired surfaces satisfy the following formula:

[0133]

[0134]

[0135] in, Let z be the z-coordinate of the i-th discrete point on the contour curve of the working surface. Let z be the z-coordinate of the corresponding point on the profile curve of the paired surface. For the initial thickness, Let be the tangent angle at the i-th discrete point; denoted as the vertical distance between the working surface geometry and the paired surface geometry at the i-th discrete point.

[0136] In one possible implementation, the determining module 52 can also be used to: determine the fitted vertical distance corresponding to each discrete point based on the contour curve of the working surface geometry and the contour curve of the paired surface geometry; and determine the thickness value of each discrete unit in the normal direction through coordinate rotation transformation based on the fitted vertical distance corresponding to each discrete point and the tangent angle of the working surface geometry at each discrete point.

[0137] In one possible implementation, the determining module 52 can also be used to: construct an initial vector based on the fitted vertical distances corresponding to each discrete point; rotate the initial vector clockwise around the origin of the coordinate system according to the tangent angle of the working surface geometry at each discrete point to obtain the rotated vector; extract the component of the rotated vector in the normal direction, and use the absolute value of the component as the thickness value of the corresponding discrete unit in the normal direction.

[0138] The leaf spring thickness determination device for a multibody dynamics model of a leaf spring provided in this application embodiment can execute the technical solution shown in the above embodiment of the leaf spring thickness determination method for a multibody dynamics model of a leaf spring. Its implementation principle and beneficial effects are similar, and will not be repeated here.

[0139] It should be noted that, for the sake of simplicity, the foregoing method embodiments are all described as a series of actions. However, those skilled in the art should understand that this application is not limited to the described order of actions, as some steps may be performed in other orders or simultaneously according to this application. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily essential to this application.

[0140] It should be further noted that although the steps in the flowchart are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowchart may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the sub-steps or stages of other steps.

[0141] It should be noted that the above-described device embodiments are merely illustrative, and the device of this application can be implemented in other ways. Furthermore, it should be understood that the division of the various modules in the above device is only a logical functional division; in actual implementation, they can be fully or partially integrated into a single physical entity, or they can be physically separated. These modules can all be implemented in software via processing element calls; they can all be implemented in hardware; or some modules can be implemented by processing element calls to software, while others are implemented in hardware. For example, the determining module can be a separately established processing element, or it can be integrated into a chip in the above device. Alternatively, it can be stored as program code in the memory of the above device, and called and executed by a processing element of the above device. The implementation of other modules is similar. In addition, these modules can be fully or partially integrated together, or they can be implemented independently. The processing element here can be an integrated circuit with signal processing capabilities. During implementation, each step of the above method or each of the above modules can be completed through integrated logic circuits in the hardware of the processor element or through software instructions.

[0142] For example, these modules can be one or more integrated circuits configured to implement the above methods, such as one or more Application Specific Integrated Circuits (ASICs), one or more Digital Signal Processors (DSPs), or one or more Field Programmable Gate Arrays (FPGAs). As another example, when a module is implemented using processing element scheduler code, the processing element can be a general-purpose processor, such as a Central Processing Unit (CPU) or other processor capable of calling program code. Furthermore, these modules can be integrated together as a System-On-a-Chip (SOC).

[0143] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented in software, it can be implemented, in whole or in part, as a computer program product. A computer program product includes one or more computer instructions. When the computer instructions are loaded and executed on a computer, all or part of the flow or function according to the embodiments of this application is generated. 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 from one computer-readable storage medium to another. For example, 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, fiber optic, Digital Subscriber Line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) 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 media can be magnetic media (e.g., floppy disks, hard disks, magnetic tapes), optical media (e.g., Digital Video Discs, DVDs), or semiconductor media (e.g., solid-state disks (SSDs)).

[0144] Figure 6 A schematic diagram of the structure of an electronic device provided as an exemplary embodiment of this application. For example... Figure 6 As shown, the electronic device 60 in this embodiment includes:

[0145] At least one processor 61; and a memory 62 communicatively connected to the at least one processor;

[0146] The memory 62 stores instructions that can be executed by at least one processor 61 to cause the electronic device to perform the method as described in any of the above embodiments.

[0147] Alternatively, the memory 62 can be either standalone or integrated with the processor 61.

[0148] The memory 62 may include high-speed random access memory (RAM) and may also include non-volatile memory, such as at least one disk storage device.

[0149] The processor 61 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application. Specifically, when implementing the method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring as described in the foregoing method embodiments, the electronic device may be, for example, an electronic device with processing capabilities such as a server.

[0150] Optionally, the electronic device may also include a communication interface 63. In specific implementations, if the communication interface 63, memory 62, and processor 61 are implemented independently, they can be interconnected via a bus to complete communication. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc., but this does not imply that there is only one bus or one type of bus.

[0151] Optionally, in a specific implementation, if the communication interface 63, memory 62, and processor 61 are integrated on a single chip, then the communication interface 63, memory 62, and processor 61 can communicate through an internal interface.

[0152] The implementation principle and technical effects of the electronic device provided in this embodiment can be found in the foregoing embodiments, and will not be repeated here.

[0153] This application also provides a computer-readable storage medium storing computer-executable instructions. When the computer-executable instructions are executed, they are used to implement the method steps as described in the above method embodiments. The specific implementation methods and technical effects are similar and will not be repeated here.

[0154] The aforementioned computer-readable storage media can be implemented from any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read Only Memory (EEPROM), Erasable Programmable Read Only Memory (EPROM), Programmable Read Only Memory (PROM), Read Only Memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The readable storage medium can be any available medium accessible to a general-purpose or special-purpose computer.

[0155] An exemplary readable storage medium is coupled to a processor, enabling the processor to read information from and write information to the readable storage medium. Of course, the readable storage medium can also be a component of the processor. The processor and the readable storage medium can reside in an application-specific integrated circuit (ASIC). Alternatively, the processor and the readable storage medium can exist as discrete components in a leaf spring thickness determination device for a multibody dynamics model of a leaf spring.

[0156] This application also provides a computer program product, including a computer program, which, when executed, implements the method steps as described in the above method embodiments. The specific implementation and technical effects are similar and will not be repeated here.

[0157] In the above embodiments, the descriptions of each embodiment have their own emphasis. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions of other embodiments. The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as the combination of these technical features does not contradict each other, it should be considered within the scope of this specification.

[0158] Other embodiments of this application will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of this application that follow the general principles of this application and include common knowledge or customary techniques in the art not disclosed herein. The specification and examples are to be considered exemplary only, and the true scope and spirit of this application are indicated by the following claims.

[0159] It should be understood that this application is not limited to the precise structure described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of this application is limited only by the appended claims.

Claims

1. A method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, characterized in that, include: Obtain the geometric parameters and number of discrete elements of the leaf spring; Based on the geometric parameters and the number of discrete elements, the thickness values ​​of each discrete element that match the target mechanical properties are determined by a publicly verifiable deterministic algorithm to form a discrete element thickness distribution. The thickness distribution of the discrete unit is output as the input file for the multibody dynamics model of the leaf spring.

2. The method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring according to claim 1, characterized in that, The geometric parameters include the contour point coordinates and initial thickness of the leaf spring. Based on these geometric parameters and the number of discrete elements, a publicly verifiable deterministic algorithm is used to determine the thickness values ​​of each discrete element that match the target mechanical properties. This includes: Curve fitting is performed based on the coordinates of the contour points to obtain the working surface geometry of the leaf spring; Based on the working surface geometry and the initial thickness, the thickness value of each discrete unit in the normal direction at its corresponding position is determined by coordinate transformation.

3. The method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring according to claim 2, characterized in that, The step of curve fitting based on the coordinates of the contour points to obtain the working surface geometry of the leaf spring includes: The contour curve of the working surface of the leaf spring is generated by performing a fourth-order polynomial fitting on the coordinates of the contour points using the least squares method.

4. The method for determining the leaf spring thickness for a multibody dynamics model of a leaf spring according to claim 2 or 3, characterized in that, The process of determining the thickness value of each discrete element in the normal direction at its corresponding position based on the geometry of the working surface and the initial thickness, through coordinate transformation, includes: Based on the working surface geometry and the initial thickness, determine the matching surface geometry corresponding to the working surface geometry; Based on the working surface geometry and the paired surface geometry, the thickness value of each discrete element in the normal direction is determined.

5. The method for determining the thickness of a leaf spring for a multibody dynamics model of a leaf spring according to claim 4, characterized in that, The step of determining the paired surface geometry corresponding to the working surface geometry based on the working surface geometry and the initial thickness includes: Based on the contour curve of the working surface geometry, determine the tangent angle at each discrete point; Based on the initial thickness and the tangent angle, determine the coordinates of the contour points of the paired surfaces; The coordinates of the contour points of the paired surfaces are fitted with a fourth-order polynomial to generate a contour curve that characterizes the geometry of the paired surfaces.

6. The method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring according to claim 5, characterized in that, The coordinates of the contour points of the paired surfaces satisfy the following formula: in, Let z be the z-coordinate of the i-th discrete point on the contour curve of the working surface. Let z be the z-coordinate of the corresponding point on the contour curve of the paired surface. The initial thickness is... Let be the tangent angle at the i-th discrete point; The vertical distance between the working surface geometry and the paired surface geometry at the i-th discrete point.

7. The method for determining the thickness of a leaf spring for a multibody dynamics model of a leaf spring according to claim 5, characterized in that, The determination of the thickness value of each discrete element in the normal direction based on the working surface geometry and the paired surface geometry includes: Based on the contour curve of the working surface geometry and the contour curve of the paired surface geometry, the fitted vertical distances corresponding to each discrete point are determined. Based on the fitted vertical distances corresponding to each discrete point and the tangent angles of the working surface geometry at each discrete point, the thickness value of each discrete unit in the normal direction is determined through coordinate rotation transformation.

8. The method for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring according to claim 7, characterized in that, The determination of the thickness value of each discrete element in the normal direction through coordinate rotation transformation includes: An initial vector is constructed based on the fitted vertical distances corresponding to each of the discrete points; The initial vector is rotated clockwise around the origin of the coordinate system according to the tangent angle of the working surface geometry at each discrete point to obtain the rotated vector. Extract the component of the rotated vector in the normal direction, and use the absolute value of the component as the thickness value of the corresponding discrete element in the normal direction.

9. A device for determining the thickness of a leaf spring in a multibody dynamics model of a leaf spring, characterized in that, include: The acquisition module is used to acquire the geometric parameters and the number of discrete elements of the leaf spring; The determination module is used to determine the thickness value of each discrete element that matches the target mechanical properties based on the geometric parameters and the number of discrete elements, using a publicly verifiable deterministic algorithm, so as to form a discrete element thickness distribution. The output module is used to output the thickness distribution of the discrete unit as the input file of the multibody dynamics model of the leaf spring.

10. An electronic device, characterized in that, include: A processor, and a memory communicatively connected to the processor; The memory is used to store computer-executed instructions; The processor is configured to execute the computer execution instructions to implement the method as described in any one of claims 1 to 8.

11. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions, which, when executed, are used to implement the method as described in any one of claims 1 to 8.

Images