A discretization method for predicting the equivalent elastic parameters of the graded lattice structure

In recent years, cellular materials have been widely studied and applied in aerospace and other fields due to the advantages of lightweight and multi-function. However, it is difficult to predict the equivalent elastic properties of the graded lattice structure and other non-uniform cellular materials because of the complex configuration and non-uniformity. A new discretization method for predicting the equivalent elastic parameters of the graded lattice structure is proposed based on the strain energy equivalent method and the discretization method in this paper. The graded lattice structure is discretized into lattice cells, the equivalent elastic properties are predicted by calculating the global equivalent elastic parameters with the parameters of lattice cells, and the calculation formulas are derived. After that, taking edge cube, face-centered cubic and body-centered cubic lattice as examples, the effectiveness and accuracy of the method are verified by theoretical calculation, numerical analysis, and experiment. The results show that the calculation errors of equivalent elastic parameters are between 4.5%–9.7%, and the errors can be significantly improved by reducing the graded factor. It proves that the proposed discretization method can predict the equivalent elastic parameters of the graded lattice structure effectively, and is suitable for different lattice structures.


Introduction
With the development of additive manufacturing technology, the research and application of cellular materials are greatly promoted. 1 Cellular materials have attracted more and more attention because of the structural characteristics such as high porosity, low relative density, and periodic arrangement, 2 which have the performance advantages of lightweight, high specific stiffness, high specific strength, and good designability. 3,4 The cellular materials are used to replace the traditional continuous material to realize the geometric structure reconstruction so that the lightweight of the structure can be effectively realized with the elastic properties unchanged. 5 Cellular materials have become a recognized lightweight and high-performance materials, which have important practical value and broad application prospects in aerospace, rail transit, engineering machinery, and other fields. 6 At present, scholars at home and abroad have already studied a lot on cellular materials. The cellular materials are divided into 2D lattice (grid structure, honeycomb structure, etc.) 7 and 3D lattice (truss structure, sandwich structure, cubic lattice structure, etc.) 8 according to the configuration rules and are divided into uniform lattice and non-uniform lattice (the graded lattice structure, variable-density lattice structure, etc.) 9 according to the material distribution. With the deepening of the research, the configuration mode, 10,11 distribution characteristics, 12 equivalent elastic properties 13,14 of the cellular materials and related research methods are gradually studied. 15,16 In order to further explore the performance advantages of cellular materials, people's focus has been extended from uniform lattice structures to nonuniform lattice structures such as graded lattice and variable-density lattice, 17 and the performance of cellular materials has expanded from structural lightweight to multi-functional characteristics, such as vibration resistance, shock resistance, heat and heat transfer, sound absorption and noise reduction, and zero/ negative thermal expansion, ect. 18,19 Seharing et al. 20 analyzed the design, mechanical behaviors, manufacturability, and application of the graded lattice structure manufactured via metallic additive manufacturing technology. Wang et al. 21 proposed a design method of cellular lattice structures combined with topology optimization technology, which realized the optimization design of overall and micro geometric structures. Jin et al. 22 carried out topology optimization on continuum geometry and realized variable-density lattice structures modeling. It was pointed out that the static performance of variable-density lattice structures was better than that of traditional uniform lattice structures. Choy et al. 23 designed and manufactured cubic lattice and honeycomb lattice structures with different support diameters and densities, and explored the energy absorption performance of lattice structures in potential impact protection applications. Liu et al. 24 have studied the dynamic response of nonlinear density gradient foam under impact. It is pointed out that adopting appropriate nonlinear density gradient distribution can improve the impact resistance and energy absorption ability of foam rod. Mridula et al. 25 and Long et al. 26 have proved through experiments that the buffer effect and energy absorption performance of gradient cellular materials are better than that of uniform cellular materials under high-speed impact, and the higher the speed, the more obvious the advantages. Dawei et al. 27 used the asymptotic homogenization method to calculate the equivalent elastic properties of various the graded lattice structure. Lin et al. 28 proposed a new topology optimization algorithm to design variable-density lattice structures, which maximized the first natural frequency of lattice structures. A large number of studies show that, compared with uniform lattice structures, the graded lattice structure and other non-uniform cellular materials have obvious advantages in vibration resistance and impact resistance.
However, the present research on the graded lattice structure and other non-uniform cellular materials have some limitations. It is difficult to accurately calculate the equivalent elastic parameters and predict the equivalent elastic properties due to the complex configuration and non-uniformity. The strain energy equivalent method and progressive homogenization method are mainly used to characterize the performance by macro equivalent method, ignoring the refinement distribution of materials, and the calculation errors are usually relatively high. The research methods mainly rely on the finite element numerical analysis and experimental testing, which are based on a large number of numerical analysis and experiments. These reasons make it difficult to predict the equivalent elastic properties of gradient lattice structure and other non-uniform cellular materials during theoretical design phase, which hinders further research.
Therefore, the graded lattice structure is taken as the research object, and the prediction method of equivalent elastic properties of the graded lattice structure is mainly explored in this paper. A new method for predicting the equivalent elastic properties of the graded lattice structure is proposed based on the strain energy equivalent method and discretization method. According to the periodic characteristics of lattice structures, the graded lattice structure is discretized into lattice cells by the discretization method. Three basic equivalent elastic parameters (equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus) are calculated to predict the equivalent elastic properties with the equivalent elastic parameters of lattice cells, and the calculation formulas of equivalent elastic parameters are derived. The concept of graded factor is defined to characterize the distribution characteristics of materials. After that, the effectiveness and accuracy of the proposed method are verified by numerical analysis and experiments, taking the edge cube (EC), face-centered cubic (FCC) and bodycentered cubic (BCC) as examples. It is hoped that this work can provide theoretical help for the research of the graded lattice structure and other non-uniform cellular materials.

Mechanical model
Different from the uniform lattice structures, the material distribution of the graded lattice structure changes along the gradient direction, and the member size and relative density decrease or increase by layers, while the lattice at the same layer presents the uniform distribution, which is a simple non-uniform lattice structure and easy to describe. 29 In order to analyze the equivalent elastic properties of the graded lattice structure, uniform lattice structures and the graded lattice structure were used to fill the given geometric structure (l z 3 l x 3 l y ). The layer along the gradient direction (zaxis) is n, the column along the horizontal direction (xaxis) is m, and the rows along the thickness direction (y-axis) is s. Compression and shear forces are applied to the geometry, and the mechanical model is shown in Figure 1. Figure 1 The left and bottom of the geometric structure are constrained by sliding support, and the freedom along the x-axis and z-axis are limited, as shown in Figure  1(a) and (c). The force F z along the z-axis is applied on the upper surface, and the compression deformation u x and u z of the geometric structure along the x-axis and z-axis are produced under the compression F z . The bottom of the geometric structure is fixed, and the freedom in the z-axis is limited, as shown in Figure 1(b) and (d). The shear force F x along the x-axis is applied on the upper surface, and the sliding deformation u xz of the geometric structure along the x-axis direction is generated under the shear force F x . It should be noted that the dotted lines in each Figure (especially in Figure(c)) indicate the deformation of geometric structure under applied loading, which does not necessarily conform to the actual geometric deformation but is only for the convenience of description.

Calculation of equivalent parameters of uniform lattice structures
The dotted line indicates that the displacement along the x-axis of the geometric structure under the compression F z is recorded as u x , and the displacement along the z-axis is recorded as u z , as shown in Figure  1(a). The dotted line indicates that the displacement along the x-axis of the geometric structure under the force F x is recorded as u xz , as shown in Figure 1(b). According to the stress-strain formulas, the elastic modulus, Poisson's ratio and shear modulus are calculated as follows: where E H , v H , and G H are elastic modulus, Poisson's ratio and shear modulus, s z and s xz are compression stress and shear stress respectively, e z , e z , and e xz are compression strain and shear strain respectively, A z is the surface area of the geometric structure, l x and l z are geometric dimensions along x-axis and z-axis. Xin et al. 15 derived the mathematical formulas of equivalent elastic parameters of uniform cellular materials based on the strain energy equivalent method and homogenization method, the formulas are as follows: where E is the elastic modulus of the original material; V is the volume of the cell, a is the envelope size of the cell; A is the cross-sectional area of the member; I is the moment of inertia; a, b and g are undetermined coefficients.
Taking BCC lattice as an example, substituting the parameters into the formulas, the formulas (2) is rewritten into formulas (3), which is the formulas for calculating the equivalent parameters of uniform BCC lattice structure.
According to formulas (2) and formulas (3), the specific calculation formulas of equivalent elastic parameters can be obtained by fitting the undetermined coefficients with the finite element method and numerical analysis. Thus, the equivalent elastic parameters of uniform lattice structures can be calculated, which is prepared for the calculation of equivalent elastic parameters of the graded lattice structure.

Discrete calculation of equivalent parameters of graded lattice
In mathematics, discretization refers to the segmentation of continuous data into a segment of discrete intervals. The discretization is used to transform the continuous problem into a discrete problem that is easy to analyze and calculate. 30,31 The discretization method is used in the calculation of the equivalent elastic parameters of graded lattice in this paper. The graded cellular materials are discretized into lattice cells by discretization method according to the periodicity and modularity of cellular materials. The equivalent elastic properties of the graded lattice structure are characterized by mathematical calculation with the equivalent elastic parameters of lattice cells. The calculation formulas of equivalent elastic parameters (equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus) of the graded lattice structure are derived. 32,33 The cell size of the graded lattice structure changes by layers along the gradient direction. The deformation and stress distribution of the graded lattice structure under certain loading also show gradient distribution characteristics. 34 According to the deformation coordination law 35 under compression, the deformation of each layer lattice cell along the gradient direction is obviously different, but the deformation of the lattice cell in the same layer is approximate, as shown in Figure 1(c). The total displacement u z of the whole the graded lattice structure along the z-axis is regarded as the deformation superposition of each layer lattice cell along the z-axis, and the total displacement u x along the x-axis is the deformation superposition of each layer lattice cell along the x-axis. The formulas of total displacement u z and u x are as follows: The sliding deformation of each layer lattice along the gradient direction of the geometric structure under the shear stress condition is obviously different, and the sliding deformation of the lattice at the same layer is approximate, as shown in Figure 1(d). The total sliding deformation u xz of the whole the graded lattice structure along the x-axis is regarded as the superposition of the sliding deformation of each layer lattice cell along the x-axis. The formulas of total sliding deformation u xz is as follows: The lattice material is equivalent to continuous material and analyzed according to the linear strain relation. The formulas of equivalent elastic modulus, Poisson's ratio and shear modulus of the graded lattice structure are derived as follows: (1) Equivalent elastic modulus (2) Equivalent Poisson's ratio In the same way, it can be concluded that: (3) Equivalent shear modulus In the same way, it can be concluded that: Therefore, the formulas of equivalent elastic parameters (equivalent elastic modulus, Poisson's ratio and shear modulus) of the graded lattice structure are as follows: According to formulas (2) and formulas (9), the design factors such as layer number, column number, cell sizes and effective effect areas are also determined, when the lattice type, material distribution and cell arrangement are determined. The equivalent elastic parameters of any the graded lattice structure can be calculated from the equivalent elastic parameters of uniform lattice structures, and the equivalent elastic properties of the graded lattice structure can be predicted.

Graded factor
To characterize the variation of material distribution along the gradient direction, the graded factor f i is defined as the ratio between the size difference of the cellular member and the initial cellular member size. The formula is as follows: where f i can be expressed as a function of i. Taking the body-centered lattice structures as an example, the envelope dimension size of the lattice cell is recorded as a, the initial cellular member size is denoted as t 0 , and the dimension of the i-layer lattice cellular member is recorded as t i . Keeping a unchanged, the relative density is changed by changing the cellular member size. The materials gradient distribution rules are constrained by the graded factor f i . When f i . 0, it is positive the graded lattice structure; when f i = 0, it is uniform lattice structures; when f i \ 0, it is negative the graded lattice structure. According to formulas (2) and formulas (9), when the design parameters such as lattice type, number of layers, initial cellular member size and graded factor are determined, the overall material distribution and structure of the graded lattice structure are determined accordingly, and the equivalent elastic parameters can be obtained by calculation.

Simulation model
In order to verify the validity and accuracy of the formulas of the equivalent elastic parameters of the graded lattice structure, three classical lattice types (edge cube (EC), face-centered cubic (FCC) and body-centered cubic (BCC)) are selected for numerical analysis.
According to the mechanical model in Section 2, the dimension of the geometry structure is set as 40 3 40 3 120 mm, F z = 2000 N, F x = 2000 N. Uniform lattice and graded lattice are used to fill the geometry, and 3D numerical analysis models are established. Nylon PA12 was used with a density of 0.95 g/cm -3 , an elastic modulus of 1.8 GPa and Poisson's ratio of 0.33. The envelope size of the lattice cell is a, and the dimension of the cellular member is t. The cellular member size of uniform lattice remains unchanged, and the cellular member of graded lattice changes along the gradient direction. According to the size of the dimension, the appropriate envelope size and member size are selected, a = 10 mm, n = 4, m = 12, s = 4. For the EC and BCC lattice, the thickness of each edge is increased by 1 mm to avoid the weak area around the edge. The parameters of the graded lattice structure are shown in Table 1.
Three kinds of numerical analysis model of the graded lattice structure are shown in Figure 2. Figure  2(a) is the edge cubic (EC) the graded lattice structure, Figure 2(b) is face-centered cubic (FCC) the graded lattice structure, and Figure 2(c) is body-centered cubic (BCC) the graded lattice structure.

Simulation analysis
The finite element software ANSYS Workbench was used to simulate and analyze the elastic properties of the graded lattice structure. The deformation and displacement results of uniform lattice structures and the graded lattice structure under the compression and shear forces are obtained by numerical analysis. The 4 Graded factor f i (i-1)/6 (i-1)/7 (i-1)/5 displacement distribution nephogram of three kinds of the graded lattice structure is shown in Figure 3. The displacement oz of the EC the graded lattice structure under the compression is shown in Figure 3(a). The displacement ox of the EC the graded lattice structure under the shear force is shown in Figure 3(b). The displacement oz of the FCC the graded lattice structure under the compression is shown in Figure 3(c). The displacement ox of the FCC cubic the graded lattice structure under the shear force is shown in Figure 3(d). The displacement oz of the BCC the graded lattice structure under the compression is shown in Figure 3(e). The displacement ox of the BCC the graded lattice structure under the shear force is shown in Figure 3(f). Other simulation results are shown in the form of curves.
The effective displacement parameters are extracted to calculate and analyze the equivalent elastic properties according to the displacement nephogram obtained by simulation. Due to the obvious discontinuity of  material distribution in lattice structures, to improve the consistency of elastic parameters of each part of the material and reduce the numerical errors caused by local deformation, the average displacement of the upper surface and the right surface is selected as the effective parameter u z , u x , u xz . The equivalent elastic parameters of uniform lattice structures can be calculated from various displacement parameters according to formulas (1). The displacement curves and equivalent elastic parameter curves of uniform lattice structures are drawn, as shown in Figure 4. The displacement curves of three kinds of uniform lattice structures is shown in Figure 4(a). The size ratio t/a of the cell is taken as the horizontal axis variable, and the displacement parameter u is taken as the longitudinal axis variable. It can be seen from the displacement curves that under the same loading, the displacement of the uniform lattice structures along all directions decreases with the increase of the cellular member size. The displacement parameters of the graded lattice structure are identified by the scatter diagram T.
According to the formulas (2), the specific formulas of equivalent elastic parameters of uniform lattice structures can be obtained by fitting undetermined coefficient. By substituting three groups of simulation results (a, t, E H , v H , G H ) into the formulas, the problem is transformed into a system of linear equations of three variables, and the undetermined coefficients can be approximately calculated, the formulas and undetermined coefficients are shown in Table 2. The fitting curves of equivalent elastic parameters (including equivalent elastic modulus, Poisson's ratio and shear modulus) of the EC, FCC and BCC uniform lattice structures are presented, in Figure 4(b)-(d). The size ratio t/a of the cell is taken as the horizontal axis variable, the equivalent elastic modulus E H and equivalent variable, the equivalent elastic modulus E H and equivalent shear modulus G H are taken as the left longitudinal variables, marked as C, with the unit of GPa. The Poisson's ratio v H is the right longitudinal variable, which is a dimensionless parameter. It can be seen from the curves that under the same loading, the equivalent elastic modulus E H and equivalent shear modulus G H of the three kinds of uniform lattice structures present a trend of increasing with the increase of cellular member size. However, the trend of the equivalent Poisson's ratio v H varies with the lattice types. The equivalent Poisson's ratio v H of the EC and the BCC uniform lattice structures increases with the increase of the cellular member size, while the equivalent Poisson's ratio v H of the FCC uniform lattice structures decreases with the increase of the cellular member size.
The equivalent elastic parameters of the graded lattice structure are calculated theoretically according to the equivalent elastic parameters of uniform lattice structures. According to Table 1 and Figure 4, the equivalent elastic parameters of uniform lattice structures are extracted. And the equivalent elastic parameters of the graded lattice structure are calculated according to the formulas (8). The simulation and calculation results of equivalent elastic properties of the graded lattice structure are shown in Figure 5. The s ec , s fcc , and s bcc are the simulation results of the equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus of EC, FCC, and BCC graded cellular materials respectively, as shown in Figure 5(a). The c ec , c fcc , and c bcc are calculation results. It can be seen from the Figure 5 that the calculation results of equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus of three kinds of the graded lattice structure are very close to the simulation results.
Simulation-calculation errors are further analyzed, as shown in Figure 5(b). The e ec-sc , e fcc-sc and e bcc-sc are simulation-calculation errors of equivalent elastic parameters of EC, FCC, and BCC the graded lattice structure respectively. It can be seen from the Figure 5 that the simulation-calculation errors of equivalent elastic parameters of three kinds of the graded lattice structure   is kept in small values. The simulation-calculation errors e ec-sc of E H T , n H T , and G H T of EC the graded lattice structure are 3.5%, 7.6%, and 5.6%. The simulationcalculation errors e fcc-sc of E H T , n H T , and G H T of FCC the graded lattice structure is 2.8%, 6.8%, and 8.1%. The simulation-calculation errors e bcc-sc of E H T , n H T , and G H T of BCC the graded lattice structure is 3.6%, 9.8%, and 7.6%. Generally speaking, the error is between 2.6%-9.8%, and the prediction results are more accurate. It is proved by simulation results that the proposed discretization method can calculate the equivalent elastic parameters to predict the equivalent elastic properties of the graded lattice structure.

Influence of graded factor
Furthermore, the effects of the size variation and material distribution inhomogeneity on the prediction of equivalent elastic properties are analyzed. By increasing the number of lattice layers in the same size range to change the gradient factor, the influence of gradient factor on the accurate and qualitative prediction of equivalent parameters of the graded lattice structures is studied.
According to the geometric structure and mechanical model in Section 2, the gradient factor is reduced by increasing the number of lattice layers. Set the external dimension of the geometric structure as l 3 40 3 120 mm, take l = 40, 50, 60, 70, 80 mm, a = 10 mm, m = 12, s = 4, n = 4, 5, 6, 7, 8. The initial member size t 0 and maximum member size t max are selected according to Table 1. The gradient factor is calculated as follows: where Dt is the size range of cellular members, Dt = t max 2t 0 . The gradient factor f i decreases with the increase of n.
The cellular member size of each layer is calculated according to formulas (11). According to formulas (9) in Section 3.2, the equivalent elastic parameters of the graded lattice structure can be calculated from the equivalent elastic parameters of uniform lattice structures. Taking BCC lattice as an example, the simulation models of the graded lattice structure with different layers is established, and the equivalent elastic properties are simulated and analyzed. The simulation and calculation results of equivalent elastic parameters of the graded lattice structure with different layers are shown in Figure 6.
The E H T ÀS , n H T ÀS , and G H T ÀS are the simulation results of equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus of gradient bodycentered lattice structures, the E H TÀS , n H T ÀS , and G H T ÀS are calculation results respectively, as shown in Figure  6(a). It can be seen from the Figure that when the number of lattice layers changes, the calculation results of equivalent elastic parameters of the graded lattice structure gradually approach the simulation results, and the curves approximately coincide. The simulationcalculation errors are analyzed as shown in Figure 6(b).
The e E-sc , e v-sc , and e G-sc are the simulation-calculation errors of equivalent elastic modulus, equivalent Poisson's ratio and equivalent shear modulus of BCC lattice structures respectively. It can be seen from Figure 6 that when the number of layers increases and the graded factor decreases, the error curves of the equivalent parameters e E-sc , e v-sc , and e G-sc of the graded lattice structure decrease gradually, and the calculation error decreases gradually. When the number of lattice layers n increases from 4 to 8 and the graded factor decreases by more than 50%, the errors are between 1.6%-5.1%, with e E-sc 1.6%, e v-sc 5.1%, and e G-sc 2.1%. The results show that the proposed method and formulas are still available, when the lattice layers n changes, and the accuracy of equivalent parameters can be effectively improved by increasing the lattice layers and decreasing the graded factor, keeping the size range of the lattice member unchanged.

Experiment
The sample of the graded lattice structure was processed by additive manufacturing technology (SLS) according to the 3D models of the graded lattice structure. The material was nylon PA12, with density 1.05 g/ cm 23 , elastic modulus 1.80 GPa, and Poisson's ratio 0.33. WDW-100 electronic universal testing machine is used to carry out the experiments of compression and shear, and the loading F = 2 KN. The accuracy of the sensor is 0.003 mm. According to the GB/T 1453-2005 and GB/T 28,889-2012, the compression and shear experiment equipment is built, as shown in Figure 7. Sample 3 is fixed by plate 2 to keep the upper plane and lower plane parallel and uniformly stressed during the compression test, as shown in Figure 7(a). Then, the sample is placed on the working platform 5, and the compression force F z is applied along the z-axis through the loading shaft 1. Displacement sensor 4 is used to collect the displacement u z and u x of the sample along z-axis and x-axis under the compression force. Sample 3 is fixed by plate 2 to keep the left and right planes of the sample parallel and uniformly stressed during the shear test, as shown in Figure 7(b). The sample is placed on the working platform 5, and the shear force F x is applied along the x-axis through the loading shaft 1. Displacement sensor 4 is used to collect the displacement u xz along the x-axis of the sample under the shear force.
The experimental samples of the graded lattice structure are shown in Figure 8. Samples 1-3 are EC, FCC and BCC graded lattice samples respectively for compression experiment; Samples 4-6 are EC, FCC and BCC graded lattice samples respectively for shear experiment.

Results and discussion
The displacement parameters of three kinds of the graded lattice structure samples along each direction under the compression force and shear force were collected in experiments, and the experimental data were analyzed. The experimental results show that when the horizontal compression loading reaches 4270 N, 10390 N, and 25,010 N respectively, and the shear load reaches 2450 N, 2860 N, and 3760 N respectively, the deformation of the EC, FCC, and BCC lattice structures sample will appear obvious yield phenomenon, and the loading-displacement curve tends to be horizontal and begins to decline. Therefore, the load of 2000 N which is far less than the yield load is selected as the target load for further analysis. The loadingdisplacement curves near the load of 2000 N are shown and approximately regarded as an elastic segment for simplified analysis.
The experimental load-displacement curves of the sample are shown in Figure 9. The u ec-z , u fcc-z , and u bcc-z are the displacements along the x-axis of the EC, FCC, and BCC the graded lattice structure samples under the compression force, as shown in Figure 9(a). The u ec-x , u fcc-x , and u bcc-x are the displacements along the x-axis of three kinds of samples under the compression force, as shown in Figure 9(b). The u ec-xz , u fcc-xz , and u bcc-xz are the displacements along the x-axis of three kinds of samples under the shear force, as shown in Figure 9(c).
The displacement parameters of the sample along each direction when the loading was of 2000 N were selected to calculate the elastic parameters. According to formulas (1), the parameters of the sample were calculated according to the deformation displacement measured, and compared with the calculation and simulation results, as shown in Figure 10(a). The s ec , s fcc , and s bcc are the simulation results of equivalent elastic parameters of three the graded lattice structure; c ec , c fcc , and c bcc are theoretical calculation results; and e ec , e fcc , and e bcc are the experimental results respectively. It can  be seen from the Figure 10 that the simulation and calculation results of the equivalent elastic parameters of the three kinds of graded lattice structure samples are very close to the experimental results.
Furthermore, calculation errors are analyzed. The e ec-ce , e fcc-ce , and e bcc-ce are the calculation -experiment errors of the three kinds of the graded lattice structure samples, as shown in Figure 10(b). It can be seen from T , and G H T of FCC the graded lattice structure is 4.5%, 7.1%, and 8.5%. The calculation-experiment errors e bcc-ce of E H T , n H T , and G H T of BCC the graded lattice structure is 8.6%, 9.7%, and 6.4%. Generally speaking, the errors are between 4.5%-9.7%, and the prediction results are more accurate. On the one hand, the results of the experiment show that the simulation method is reliable for the study of the equivalent performance of the graded lattice structure. On the other hand, it is proved that the proposed discretization method can accurately predict the equivalent elastic properties of the graded lattice structure.
Of course, there are some shortcomings in this paper. The amount of work and research content is limited in this paper, with three types of lattice studied. The size of the design dimension and the number of experimental samples is small. It is worth noting that the accuracy of the prediction of equivalent elastic properties of the graded lattice structure is verified by taking three lattice types as examples, but the influence of lattice layers and  graded factor changes on the accuracy of equivalent elastic performance prediction is not verified in the experiment verification part. However, according to the results obtained, the simulation results of the equivalent elastic parameters of the three kinds of the graded lattice structure samples are very close to the experimental results, which proves the reliability of the simulation method. It is reasonable to believe that the results obtained by theory and numerical analysis in this paper are credible. In addition, the conclusion in this paper is inevitably limited, and the relevant influencing factors and laws need to be further explored and improved by other researchers.

Conclusions
The graded lattice structure is taken as the research object to explore the prediction method of equivalent elastic properties of non-uniform cellular materials in this paper. A discretization method for calculating the equivalent elastic parameters of the graded lattice structure is proposed. The theoretical calculation, numerical analysis and experimental verification of the equivalent elastic parameters of the graded lattice structure are carried out. According to the research, the conclusions are as follows: (1) The graded lattice structure is discretized into lattice cells by the discretization method, and the calculation formulas of equivalent elastic parameters of the graded lattice structure are derived. The formulas show that the equivalent elastic properties of the graded lattice structure can be predicted by mathematical calculation of the equivalent elastic parameters of uniform lattice structures. (2) Taking EC, FCC and BCC lattice as examples, the equivalent elastic properties of the graded lattice structure are simulated and calculated. The results show that the calculation results of equivalent elastic parameters of the graded lattice structure are consistent with the simulation results, and the errors are between 2.6%-9.8%. And it has been proved that the errors can be significantly improved by increasing the number of lattice layers and reducing the graded factor. When the number of lattice layers n increases from 4 to 8 and the graded factor decreases by more than 50%, the errors are between 1.6%-5.1%. (3) Three kinds of the graded lattice structure (EC, FCC, and BCC) are processed by additive manufacturing technology (SLS), and the static compression and shear tests are carried out. The experiment results show that the simulation results, theoretical calculation results and experimental results are closely approximate, and the calculation-experiment errors are kept between 4.5%-9.7%. The experimental results show that the proposed discretization method can calculate the equivalent elastic parameters and realize the prediction of the equivalent elastic properties of the graded lattice structure effectively.
The discretization method proposed in this paper can predict the equivalent elastic parameters of the graded lattice structure during the theoretical design phase, and reduce the workload of simulation and experiment. This work can provide theoretical support for further research and application of the graded lattice structure and other non-uniform cellular materials