# A Kind of Filling Lattice Structure Based on DFAM for Mechanical Parts: The Diamond Lattice Structure

Authors: Xuan Yin; Wenjun Meng; Jinzhao Cheng; Hailong Wang; Xiaoxia Zhao
Journal: Advances in Materials Science and Engineering
 (2021)

---

## Abstract

Thanks to the geometric and material complexity of additive manufacturing, the design space of mechanical parts has been developed, in which lattice filling structure customization can be applied to the solid filling of mechanical parts to achieve the goal of mechanical structure lightweight. A kind of diamond lattice structure unit is designed by imitating the natural method based on Design for Additive Manufacturing of mechanical parts. The mathematical model of the relative density and mechanical properties of the unit are established, and the relationship between the two is obtained, which is verified by simulations; then the relatively uniform results are obtained. The variable density hypothesis of diamond lattice structure is proposed, the methods of simulations and compression tests are used to verify the hypothesis, and the results show that the variable density structure with the density of the filling element decreasing gradually with the stress point as the center has better compression performance and concurrently verify the correctness and applicability of the equivalent modulus of elasticity mathematical model. The results of this study can be applied to the solid sandwich filling of pressure mechanical parts, and the stress density matching relationship can be carried out to further specific design.


---

## Body

## 1. Introduction

The characteristics of additive manufacturing, such as geometric complexity and material complexity, provide new opportunities for customization of products, which also develop new design space for designers. The concept of Design for Additive Manufacturing (DFAM) emerged.Design for Additive Manufacturing (DFAM) is defined as follows [1, 2]: The process of achieving expected performance and other life cycle objectives within the capability of additive manufacturing through the synthesis of shape, size, hierarchy, and material composition. DFAM, which gives full play to the customizable and multifunctional characteristics of additive manufacturing, focuses on maximizing the functional value of products. The design freedom of additive manufacturing [3] is shown in Figure 1.Figure 1
Design freedom of additive manufacturing.Different from the reduction process, additive manufacturing allows complex structural adjustment from the inside without additional manufacturing costs. Therefore, the application of sandwich structure as an important embodiment of additive manufacturing structure design is conducive to the lightweight of parts. Lattice structure is a kind of structure which is easy to be parameterized and controlled. Taking it as the basic unit of sandwich structure is the basis of this paper.Lattice structure is a kind of periodically connected space filling element, which is usually created by truss structure and minimalist surface [4], which not only has the advantages of saving materials, reducing energy use, and saving time, but also has high energy absorption rate, high-strength weight ratio, and thermal management ability. The lattice structure has not been widely used in the past due to the limitation of processing methods, and the rise of additive manufacturing makes mesoscopic lattice structure get a lot of attention and research.Rosen [2] and others [1, 5] of Georgia Institute of Technology in the United States have carried out systematic research on DFAM and proposed systematic design method of cellular structure. Li Yang and others of the University of Louisville in the United States have carried out theoretical and experimental research on metal cellular unit made of AM and summarized the size effect of cellular structure [6, 7]. Sustainable Manufacturing and Life Cycle Engineering Research Group of the University of New South Wales in Australia is also committed to the research of lattice structure; the design, analysis and manufacture of lattice structure are reviewed [4]; in addition, it cooperates with Supachai Vongbunyong in Thailand to develop the integrated software tool for cellular structure design [8]. The research team of Universidade de Lisboa in Portugal has done in-depth research on the mechanical properties of the lattice structure [9, 10]. Ahmed Yussuf Hussein of the University of Exeter in UK has studied the development of light porous structure of metal additive manufacturing [11]; in addition, many scholars have done a lot of research on the additive manufacturing and mechanical properties of porous metal biomaterials [12, 13]. Panda of the National Institute of Technology in India designed and developed the cellular structure of additive manufacturing [14]. In China, Beijing University of Science and Technology is dedicated to the research of spatial porous structure; Song et al. proposed an effective optimization method for the design of irregular cellular structure [15]. Zhao et al. studied the design method of variable-density porous structure based on local relative density mapping [16]. Liu et al. from Beijing University of Aeronautics and Astronautics studied the design method of honeycomb structure generation for additive manufacturing [17]. Liao et al. of South China University of Technology proposed a topological optimization-based variable density lattice structure optimization design method [18].It is not difficult to find that the research on the lattice structure of additive manufacturing is relatively limited at this stage. Most of the lattice structure units are conventional which are generated manually; few of them are optimized design. Most of the lattice structures are in periodic arrangement, lacking the variable scale mapping of the functional requirements. In addition, the experimental research on the lattice structure is very rare, so the simulation analysis cannot be verified, and the parameter optimization of the lattice structure has no experience to follow.This paper focuses on the theoretical and experimental study of the lattice structure with variable density based on compressive stress. In this paper, a kind of optimized lattice structure unit is proposed. A fine mathematical model of lattice structure density and elastic modulus is established on the premise of considering deformation cross section. On this basis, a series of simulation and experiment methods are used to verify the correctness of the above theory.


## 2. Lattice Structure Unit Design

The spatial lattice structure mentioned in the existing literature includes octahedral structure, body centered cubic structure, and orthotropic cubic structure. There are many rods connected at the joints of these structures, which is easy to cause stress concentration and rod waste. Furthermore, in the context of additive manufacturing, the lattice structure unit should not be limited to the traditional regular structure, but should be combined with the load-bearing situation of parts to find the most suitable internal lattice structure unit under specific load.There are some natural structures in nature, which have unique advantages under some stress conditions. For example, the spider web structure has great tension, the laminated structure of plant stem makes its bending and torsion resistance remarkable, the vein structure of lotus leaf makes its pressure bearing capacity extremely strong, and the honeycomb structure has a certain strength while playing its accommodation function.Diamond is one of the substances with the highest hardness in the world; because each C atom in its crystal is connected with four C atoms around it by covalent bond with the same angle, the inter atomic force is quite firm, which makes the material very stable. This kind of natural and superior molecular structure is worth learning from in the design of lattice structure. In this paper, a preliminary study on the lattice elements of pressure parts is carried out by using the simulated natural design.


### 2.1. Three-Dimensional Modeling of Diamond Lattice Structure

The diamond lattice structure unit (DLSU) designed in this paper is shown in Figure2. C atom is replaced by node and covalent bond is replaced by rod. Each node and four surrounding nodes are connected by rods to form space tetrahedron structure. The lengths of rods are equal, and the angle between rods is 109° 28 ′(109.47°). Fillet transition connection is designed between nodes and rods in order to increase the stability of slender members. This kind of structure with space center symmetry has relatively uniform internal force, and only four rods are connected to each node, resulting in a smaller stress concentration at the node, which is more robust than the ordinary structure.Figure 2
The diamond lattice structure unit.In order to divide the solid into space grids in the design and parametric modeling of the later structure, the concept of structural cell is introduced. The space of each structural unit is Hexagonal Prism (HP) cell, as shown in Figure3; the structural cells are connected end to end to fill the solid structure without gaps.where ∅ is the diameter of the inner tangent circle of the hexagonal prism; h is the height of HP cell; b is the bottom edge length of HP cell; and l is the rod length of the DLSU. After calculation, the parameter relationships are as follows:Figure 3
HP cell and its superposition.(1)h=63∅,l=64∅,b=33∅,


### 2.2. Calculation of Cell’s Relative Density of Diamond Lattice Structure

According to the Gibson–Ashby micromechanical model of porous structure (1997), the equivalent mechanical properties of lattice structure, such as elastic modulus, shear modulus, yield strength, etc., are all represented by polynomials related to the relative density of lattice structure. The relative density of lattice structure can be expressed as(2)ρx=VsVc,where ρx represents the relative density of structural cell; Vs is the volume of solid part in structural cell; and Vc is the total volume of structural cell, for HP cell Vc=2/2∅3.


#### 2.2.1. The Case of Low Density

The rods constituting the unit structure are slender rods with poor stability when the relative density is low. Therefore, the lattice unit structure with spherical nodes and fillet transition is adopted. The schematic diagram is shown in Figure4, in which the volume of equivalent rod is accurately calculated by equation (3).Figure 4
Model of fillet connection between node and rod.(3)Vg=πr2⋅l−2R−rR+r+2Ry+2π∫0mg2y−∫R−rR+r+2Ry−Rmf2ydy.In the design, the rod is the main force object, and its volume and mechanical properties are changed by changing its section radiusr. The node, as the auxiliary structure of stress, the determination of its size mainly follows the principle of noninterference of intersecting line. Therefore, the node radius R should be determined by the rod radius r and the fillet radius Ry; that is, R=6/2r+6−2/2Ry.The calculation process of relative density of structure cell is shown in Table1. The fillet radius Ry=r is selected to simplify the calculation; the final results are as follows:(4)ρx=10.8833r∅2−20.1912r∅3.Table 1
Calculation process of relative density.
ParametersHP cellMesh parameterInscribed circle diameter∅Solid volume parametersRod section radiusr, joint fillet radius Ry, mesh parameter ∅Number of equivalent nodes (nj)2Number of equivalent rods (ng)4Length of equivalent rod (l)6/4∅Length of rod transition area (m)Ry6−2r6r+2Ry/6−1r+RyVolume of equivalent node (Vj)43π62r+6−22Ry3Volume of equivalent rod (Vg)πr2⋅l−2R−rR+3r+2π∫0mg2ydy−∫R−rR+3r−Rmf2ydySolid volume (Vs)Vs = njVj + ngVgStructure cell volume (Vc)2/2∅3Relative density (ρx)2njVj+ngVg/∅3It can be seen that the value of relative density is a cubic function of the ratio of the section radius of the rodr to the cell parameter∅. When the ratio increases gradually, the space interference of the nodes geometric structure will be formed, which should be avoided by controlling the ratio in the range of 0≤r/∅≤3/8.


#### 2.2.2. The Case of High Density

The spherical node is changed into an ordinary node as the stability increases with the increasing relative density of the unit; that is, the rods are directly connected. The structure unit and the analyzed equivalent rods are shown in Figure5, in which (a) shows the basic unit of diamond lattice structure after omitting the spherical node and fillet; (b), (c), and (d) show the equivalent rod and its partial enlarged drawing after unit superposition. It can be seen from Figure 5(d) that the two ends of the equivalent rods are all three equal sector sections with variable radius. As shown in Figure 5(e), the equivalent rod is made up into a cylindrical rod in order to calculate the volume. 5(f) is shown to establish a Cartesian coordinate system for the complete cylindrical rod, by which calculating the volume of the equivalent rod shown in Figure 5(c), and thus establishes a mathematical model for the relative density of lattice DLSU with ordinary nodes.Figure 5
DLSU with ordinary nodes and their equivalent rods.In Figure5(f), the sector plane A is expressed asz=x/tan109.47°/2, the Cartesian coordinate system is transformed into polar coordinate system for the convenience of calculation, where x=ρcosθ, y=ρsinθ, and dxdy=ρdρdθ, and the complement part is divided, calculated, and then removed to obtain the volume of the equivalent rod:(5)Vg′=πr2l−12∫0π/3∫0rρcosθtan109.47°/2ρdρdθ,where r is the section radius of the rod, l is the length of the rod (related to the cell size), and the integral region D is the sector region with radius r and included angle of 120°.The relative density of HP cell isρx=42Vg′/∅3. The calculation result is as follows:(6)ρx=23πr∅2−86tan109.47°2−1r∅3.From the properties of the function in formulae (5) and (8), it can be seen that the ratio of the section radius r of the element member to the cell size parameter (a or ∅) is an important value, which is recorded in this paper as parameter ratio. The larger the parameter ratio is, the closer the structure is to the solid. When the relative density ρx equals 1, the structure appears as a solid in the macroscopic aspect, while the parameter ratio should be within the range of r/∅≤0.5239, which makes the relative density ρx meaningful.


#### 2.2.3. Cell’s Relative Density of Diamond Lattice Structure

The conclusions in sections 2.2.1 and 2.2.2 are combined to form the mathematical models of the relative density of the HP cell of the diamond lattice structure, as shown in the following equation:(7)ρx=10.8833r∅2−20.1912r∅3,0<r∅≤38,23πr∅2−86tan109.47°2−1r∅3,38<r∅≤0.5239.The function valueρx obtained by different parameter combinations of the combined mathematical model varies from 0 to 1. Plot the function to get the result shown in Figure 6.Figure 6
Final results of structure cell relative density.The contour lines of function value shown in horizontal plane 0 in Figure6 are a series of straight lines passing through the origin, which indicates that when the parameter ratio r/∅ is fixed, the relative density is the same, and with the increase of parameter ratio, the relative density increases nonlinearly.


## 2.1. Three-Dimensional Modeling of Diamond Lattice Structure

The diamond lattice structure unit (DLSU) designed in this paper is shown in Figure2. C atom is replaced by node and covalent bond is replaced by rod. Each node and four surrounding nodes are connected by rods to form space tetrahedron structure. The lengths of rods are equal, and the angle between rods is 109° 28 ′(109.47°). Fillet transition connection is designed between nodes and rods in order to increase the stability of slender members. This kind of structure with space center symmetry has relatively uniform internal force, and only four rods are connected to each node, resulting in a smaller stress concentration at the node, which is more robust than the ordinary structure.Figure 2
The diamond lattice structure unit.In order to divide the solid into space grids in the design and parametric modeling of the later structure, the concept of structural cell is introduced. The space of each structural unit is Hexagonal Prism (HP) cell, as shown in Figure3; the structural cells are connected end to end to fill the solid structure without gaps.where ∅ is the diameter of the inner tangent circle of the hexagonal prism; h is the height of HP cell; b is the bottom edge length of HP cell; and l is the rod length of the DLSU. After calculation, the parameter relationships are as follows:Figure 3
HP cell and its superposition.(1)h=63∅,l=64∅,b=33∅,


## 2.2. Calculation of Cell’s Relative Density of Diamond Lattice Structure

According to the Gibson–Ashby micromechanical model of porous structure (1997), the equivalent mechanical properties of lattice structure, such as elastic modulus, shear modulus, yield strength, etc., are all represented by polynomials related to the relative density of lattice structure. The relative density of lattice structure can be expressed as(2)ρx=VsVc,where ρx represents the relative density of structural cell; Vs is the volume of solid part in structural cell; and Vc is the total volume of structural cell, for HP cell Vc=2/2∅3.


### 2.2.1. The Case of Low Density

The rods constituting the unit structure are slender rods with poor stability when the relative density is low. Therefore, the lattice unit structure with spherical nodes and fillet transition is adopted. The schematic diagram is shown in Figure4, in which the volume of equivalent rod is accurately calculated by equation (3).Figure 4
Model of fillet connection between node and rod.(3)Vg=πr2⋅l−2R−rR+r+2Ry+2π∫0mg2y−∫R−rR+r+2Ry−Rmf2ydy.In the design, the rod is the main force object, and its volume and mechanical properties are changed by changing its section radiusr. The node, as the auxiliary structure of stress, the determination of its size mainly follows the principle of noninterference of intersecting line. Therefore, the node radius R should be determined by the rod radius r and the fillet radius Ry; that is, R=6/2r+6−2/2Ry.The calculation process of relative density of structure cell is shown in Table1. The fillet radius Ry=r is selected to simplify the calculation; the final results are as follows:(4)ρx=10.8833r∅2−20.1912r∅3.Table 1
Calculation process of relative density.
ParametersHP cellMesh parameterInscribed circle diameter∅Solid volume parametersRod section radiusr, joint fillet radius Ry, mesh parameter ∅Number of equivalent nodes (nj)2Number of equivalent rods (ng)4Length of equivalent rod (l)6/4∅Length of rod transition area (m)Ry6−2r6r+2Ry/6−1r+RyVolume of equivalent node (Vj)43π62r+6−22Ry3Volume of equivalent rod (Vg)πr2⋅l−2R−rR+3r+2π∫0mg2ydy−∫R−rR+3r−Rmf2ydySolid volume (Vs)Vs = njVj + ngVgStructure cell volume (Vc)2/2∅3Relative density (ρx)2njVj+ngVg/∅3It can be seen that the value of relative density is a cubic function of the ratio of the section radius of the rodr to the cell parameter∅. When the ratio increases gradually, the space interference of the nodes geometric structure will be formed, which should be avoided by controlling the ratio in the range of 0≤r/∅≤3/8.


### 2.2.2. The Case of High Density

The spherical node is changed into an ordinary node as the stability increases with the increasing relative density of the unit; that is, the rods are directly connected. The structure unit and the analyzed equivalent rods are shown in Figure5, in which (a) shows the basic unit of diamond lattice structure after omitting the spherical node and fillet; (b), (c), and (d) show the equivalent rod and its partial enlarged drawing after unit superposition. It can be seen from Figure 5(d) that the two ends of the equivalent rods are all three equal sector sections with variable radius. As shown in Figure 5(e), the equivalent rod is made up into a cylindrical rod in order to calculate the volume. 5(f) is shown to establish a Cartesian coordinate system for the complete cylindrical rod, by which calculating the volume of the equivalent rod shown in Figure 5(c), and thus establishes a mathematical model for the relative density of lattice DLSU with ordinary nodes.Figure 5
DLSU with ordinary nodes and their equivalent rods.In Figure5(f), the sector plane A is expressed asz=x/tan109.47°/2, the Cartesian coordinate system is transformed into polar coordinate system for the convenience of calculation, where x=ρcosθ, y=ρsinθ, and dxdy=ρdρdθ, and the complement part is divided, calculated, and then removed to obtain the volume of the equivalent rod:(5)Vg′=πr2l−12∫0π/3∫0rρcosθtan109.47°/2ρdρdθ,where r is the section radius of the rod, l is the length of the rod (related to the cell size), and the integral region D is the sector region with radius r and included angle of 120°.The relative density of HP cell isρx=42Vg′/∅3. The calculation result is as follows:(6)ρx=23πr∅2−86tan109.47°2−1r∅3.From the properties of the function in formulae (5) and (8), it can be seen that the ratio of the section radius r of the element member to the cell size parameter (a or ∅) is an important value, which is recorded in this paper as parameter ratio. The larger the parameter ratio is, the closer the structure is to the solid. When the relative density ρx equals 1, the structure appears as a solid in the macroscopic aspect, while the parameter ratio should be within the range of r/∅≤0.5239, which makes the relative density ρx meaningful.


### 2.2.3. Cell’s Relative Density of Diamond Lattice Structure

The conclusions in sections 2.2.1 and 2.2.2 are combined to form the mathematical models of the relative density of the HP cell of the diamond lattice structure, as shown in the following equation:(7)ρx=10.8833r∅2−20.1912r∅3,0<r∅≤38,23πr∅2−86tan109.47°2−1r∅3,38<r∅≤0.5239.The function valueρx obtained by different parameter combinations of the combined mathematical model varies from 0 to 1. Plot the function to get the result shown in Figure 6.Figure 6
Final results of structure cell relative density.The contour lines of function value shown in horizontal plane 0 in Figure6 are a series of straight lines passing through the origin, which indicates that when the parameter ratio r/∅ is fixed, the relative density is the same, and with the increase of parameter ratio, the relative density increases nonlinearly.


## 2.2.1. The Case of Low Density

The rods constituting the unit structure are slender rods with poor stability when the relative density is low. Therefore, the lattice unit structure with spherical nodes and fillet transition is adopted. The schematic diagram is shown in Figure4, in which the volume of equivalent rod is accurately calculated by equation (3).Figure 4
Model of fillet connection between node and rod.(3)Vg=πr2⋅l−2R−rR+r+2Ry+2π∫0mg2y−∫R−rR+r+2Ry−Rmf2ydy.In the design, the rod is the main force object, and its volume and mechanical properties are changed by changing its section radiusr. The node, as the auxiliary structure of stress, the determination of its size mainly follows the principle of noninterference of intersecting line. Therefore, the node radius R should be determined by the rod radius r and the fillet radius Ry; that is, R=6/2r+6−2/2Ry.The calculation process of relative density of structure cell is shown in Table1. The fillet radius Ry=r is selected to simplify the calculation; the final results are as follows:(4)ρx=10.8833r∅2−20.1912r∅3.Table 1
Calculation process of relative density.
ParametersHP cellMesh parameterInscribed circle diameter∅Solid volume parametersRod section radiusr, joint fillet radius Ry, mesh parameter ∅Number of equivalent nodes (nj)2Number of equivalent rods (ng)4Length of equivalent rod (l)6/4∅Length of rod transition area (m)Ry6−2r6r+2Ry/6−1r+RyVolume of equivalent node (Vj)43π62r+6−22Ry3Volume of equivalent rod (Vg)πr2⋅l−2R−rR+3r+2π∫0mg2ydy−∫R−rR+3r−Rmf2ydySolid volume (Vs)Vs = njVj + ngVgStructure cell volume (Vc)2/2∅3Relative density (ρx)2njVj+ngVg/∅3It can be seen that the value of relative density is a cubic function of the ratio of the section radius of the rodr to the cell parameter∅. When the ratio increases gradually, the space interference of the nodes geometric structure will be formed, which should be avoided by controlling the ratio in the range of 0≤r/∅≤3/8.


## 2.2.2. The Case of High Density

The spherical node is changed into an ordinary node as the stability increases with the increasing relative density of the unit; that is, the rods are directly connected. The structure unit and the analyzed equivalent rods are shown in Figure5, in which (a) shows the basic unit of diamond lattice structure after omitting the spherical node and fillet; (b), (c), and (d) show the equivalent rod and its partial enlarged drawing after unit superposition. It can be seen from Figure 5(d) that the two ends of the equivalent rods are all three equal sector sections with variable radius. As shown in Figure 5(e), the equivalent rod is made up into a cylindrical rod in order to calculate the volume. 5(f) is shown to establish a Cartesian coordinate system for the complete cylindrical rod, by which calculating the volume of the equivalent rod shown in Figure 5(c), and thus establishes a mathematical model for the relative density of lattice DLSU with ordinary nodes.Figure 5
DLSU with ordinary nodes and their equivalent rods.In Figure5(f), the sector plane A is expressed asz=x/tan109.47°/2, the Cartesian coordinate system is transformed into polar coordinate system for the convenience of calculation, where x=ρcosθ, y=ρsinθ, and dxdy=ρdρdθ, and the complement part is divided, calculated, and then removed to obtain the volume of the equivalent rod:(5)Vg′=πr2l−12∫0π/3∫0rρcosθtan109.47°/2ρdρdθ,where r is the section radius of the rod, l is the length of the rod (related to the cell size), and the integral region D is the sector region with radius r and included angle of 120°.The relative density of HP cell isρx=42Vg′/∅3. The calculation result is as follows:(6)ρx=23πr∅2−86tan109.47°2−1r∅3.From the properties of the function in formulae (5) and (8), it can be seen that the ratio of the section radius r of the element member to the cell size parameter (a or ∅) is an important value, which is recorded in this paper as parameter ratio. The larger the parameter ratio is, the closer the structure is to the solid. When the relative density ρx equals 1, the structure appears as a solid in the macroscopic aspect, while the parameter ratio should be within the range of r/∅≤0.5239, which makes the relative density ρx meaningful.


## 2.2.3. Cell’s Relative Density of Diamond Lattice Structure

The conclusions in sections 2.2.1 and 2.2.2 are combined to form the mathematical models of the relative density of the HP cell of the diamond lattice structure, as shown in the following equation:(7)ρx=10.8833r∅2−20.1912r∅3,0<r∅≤38,23πr∅2−86tan109.47°2−1r∅3,38<r∅≤0.5239.The function valueρx obtained by different parameter combinations of the combined mathematical model varies from 0 to 1. Plot the function to get the result shown in Figure 6.Figure 6
Final results of structure cell relative density.The contour lines of function value shown in horizontal plane 0 in Figure6 are a series of straight lines passing through the origin, which indicates that when the parameter ratio r/∅ is fixed, the relative density is the same, and with the increase of parameter ratio, the relative density increases nonlinearly.


## 3. Equivalent Mechanical Properties of Diamond Lattice Structure

### 3.1. Equivalent Elastic Modulus of Diamond Lattice Structure

As shown in Figure7(a), stress σ is applied on the cell of diamond lattice structure, which results in deformation of lattice structure. The dotted line in Figure 7(b) represents the original structural axis, and the solid line represents the deformed structural axis. Since rods are rigidly connected with each other, the included angle between the tangent lines of each rod end at the rigid node O remains unchanged after deformation. In order to obtain the equivalent elastic modulus, the displacement at point A is required, that is, the total deformation of the unitΔl.Figure 7
Stress and deformation diagram of structure cell.
(a)(b)According to the functional principle, the work done by the pressure is equal to the stored energy of rods [19]; that is, the work done by the external force is equal to the sum of the strain energy of the rod system in numerical value, which is expressed as follows:(8)12FΔl=U1+U2,where U1 is the sum of strain energy in tension and compression and U2 is the sum of strain energy in bending. It shows that the total deformation of the structure unit can be obtained by calculating the strain energy of each rod. The specific analysis is as follows.Figure8(a) is the force analysis of the rod. There is no slip in the horizontal direction on node O for its resultant force in the horizontal direction is 0, so that node O displaces only in the vertical direction. In order to learn the strain energy conveniently, the node O is set as the fixed point, and the rods OA and OB are taken as the research objects for the force and deformation analysis.Figure 8
Force analysis and profile simplification of rods.
(a)(b)The resultant force of OA rod is an axial force ofF, which leads to tension and compression deformation. For OB rod, the force FB is decomposed into axial force and vertical force to obtain the forces N and P, respectively, and the rod is subject to axial force and bending moment M, resulting in tension and compression deformation and bending deformation. According to the symmetry of the structure and the uniformity of the force, the force condition of rods OC and OD is the same as that of OB. Therefore, tension and compression strain energy is produced in all the rods of the deformed unit, and bending strain energy is produced in the rods OB, OC, and OD.The equations of strain energy in tension and compression and in bending are, respectively, as follows:(9)U1=∫0lF2x2EAxdx,(10)U2=∫0lM2x2EIdx,where E represents the inherent elastic modulus of the raw material.F (x) in equation (9) is a function of the force on the rod, which is a constant function in this example. Ax is a function of the cross-sectional area of the rod. Due to the complexity of the change of the cross-sectional dimension of the designed member, as shown in Figure 8(b), the profile of the rod is simplified from the change of curved surface to the change of wedge shape for the convenience of calculation. In this research, (11)Ax=π6−1r−3−2x2,0≤x≤2r,πr2,2r<x<l−2r,π6−1r−3−2l−x2,l−2r≤x≤l.In equation (10), Mx is a function of bending moment, which of OB rod is(12)Mx=Px−l2.The value ofP is FBcos 19.47°. I is the section moment of inertia of the rod, that of cylinder section is I0=πr4/4; for the wedge section, the reduced moment of inertia of the variable section rod Iz is calculated as follows [20]:(13)Iz=161/I1+2/I2+9/I3+4/I4=πr41.3159,where I1, I2, I3, and I4 are the section moments of inertia when x=32/4r, x=2/2r, x=2/4r, and x=0, respectively, as shown in Figure 8(b).In conclusion, the total tensile and compressive strain energy of the structure unit is as follows:(14)U1=U1OA+3U1OB=∫0lF22EAxdx+3∫0lN22EAxdx.The value ofN is FBsin 19.47°. The total bending strain energy of the rods in the structure unit is expressed as follows:(15)U2=3U2OB=3∫0lM2x2EIdx=3∫02rM2x2EIzdx+∫2rl−2rM2x2EI0dx+∫l−2rlM2x2EIzdx.From the above calculation process, it can be observed that rod OA mainly bears load by tension and compression deformation, while rods OB, OC, and OD have very small tension and compression deformation and mainly bear load by bending deformation.From the above equations (8) ∼ (15), the total deformation of the structure can be obtained. That is,(16)Δl=σπE0.0196∅5r4−0.1826∅4r3+2.0220∅3r2−2.1427∅2r.Therefore, the equivalent elastic modulusE∗ of diamond lattice structure is obtained by equation (17)(17)E∗=σε=σΔhΔl,where ε represents the strain produced by the stress of the structure cell and h is the original height of the cell, whose value is 6/3∅. After sorting out, the following result is obtained:(18)E∗=E0.0077∅/r4−0.0712∅/r3+0.7883∅/r2−0.8353∅/r.Obviously, the equivalent elastic modulus of diamond lattice structure is related to the elastic modulusE of the raw material itself and the parameter ratio of structure cell ∅/r; the parameter ratio here is the inverse of that described in Section 2. The value of E/E∗ is a quartic function of the parameter ratio ∅/r, while from the above study, it is concluded that the relative density is a cubic function of the parameter ratio ∅/r; it is shown that the equivalent elastic modulus E∗ has a nonlinear relationship with the relative density ρx. Therefore, when the parameter ratio is fixed, the relative density and the equivalent elastic modulus of diamond lattice structure are fixed, respectively. After calculation, the following relationship is obtained:(19)E∗=E0.4239ρx3/2−0.0008.


### 3.2. Simulation of Mechanical Properties of Diamond Lattice Structure Unit

For lattice structure, the stress-strain curve changes with the increase of relative density, and the corresponding elastic modulus, plastic yield stress, and ultimate stress of brittle fracture increase [21]. With the increase of the parameter ratio r/∅, the density of diamond structure increases, and the mechanical properties change.A structure unit is intercepted to simulate the upper unit transferring three equal forcesF along the direction of the rods, and the stress state of the unit under different relative density is observed. When the maximum stress of the structure reaches the yield stress, the value of F reaches the upper limit, which is recorded as Fmax. Aluminum alloy 1061 was selected as the simulation material, with yield stress of 27.57 MPa and elastic modulus of 69 GPa. Figure 9 shows the mechanical simulation analysis of diamond lattice structure units under different rod section radius (0.5 mm, 1.0 mm, 1.5 mm, and 2 mm) when the inner tangent diameter of HP cell section ∅ is 10 mm. It can be seen from the von Mises stress that the most dangerous point of the structure occurs at the joint between the node and the rod.Figure 9
Simulation analysis of diamond lattice structure unit.As seen in Table2 which shows the comparison of simulation result data, the larger the section radius r of the rod, the greater the relative density of the structure, the greater the load-bearing capacity, and the smaller the deformation of the structure; all these changes are nonlinear as shown in Figure 10.Table 2
Comparison of simulation data of diamond lattice structure unit.
Serial numberRod section radiusr (mm)Relative densityρxForceFmax (N)Maximum stressσs(MPa)Maximum displacements (mm) (×10−3)Simulation equivalent modulus of elasticityE1∗ (GPa)Theoretical equivalent modulus of elasticityE2∗ (GPa)10.50.02653.64427.575.6020.06170.071020.60.03756.57327.574.9130.12610.157230.70.050710.70027.574.4070.22890.278740.80.065715.90227.573.9710.37760.437450.90.082524.31927.573.8110.60160.637961.00.101132.34627.573.6100.84480.885071.10.121442.69727.573.4401.17021.182081.20.143356.20127.573.3711.57181.531591.30.166970.79727.573.2642.04501.9391101.40.192188.59227.573.2182.59562.4075111.50.2188110.53927.573.2313.22562.9383121.60.2469130.45027.573.1123.95213.5331131.70.2765160.65027.573.1734.77354.1974141.80.3075190.28027.573.1475.70064.9323151.90.3398219.68127.573.0636.76195.73841620.3734250.63627.572.9667.96716.6186Figure 10
Variation trend of parameters of structure unit with section radius of rod.Comparing the equivalent elastic modulusE1∗ calculated by the simulation data with the equivalent elastic modulus E2∗ calculated by the above theoretical research, as shown in Figure 11, the overall trend of the curve is consistent.Figure 11
Comparison of theoretical and simulation results of equivalent modulus of elasticity.However, with the increase of the parameter ratio, the bending moment condition of OB, OC, and OD rods in the model changes into the shear condition. Therefore, the actual equivalent elastic modulus of the lattice structure is higher than the calculated value of the theoretical model in the region of high parameter ratio, as shown in Figure11. The existence of this gap makes it necessary to adjust the theoretical model for the high parameter ratio in future work.


## 3.1. Equivalent Elastic Modulus of Diamond Lattice Structure

As shown in Figure7(a), stress σ is applied on the cell of diamond lattice structure, which results in deformation of lattice structure. The dotted line in Figure 7(b) represents the original structural axis, and the solid line represents the deformed structural axis. Since rods are rigidly connected with each other, the included angle between the tangent lines of each rod end at the rigid node O remains unchanged after deformation. In order to obtain the equivalent elastic modulus, the displacement at point A is required, that is, the total deformation of the unitΔl.Figure 7
Stress and deformation diagram of structure cell.
(a)(b)According to the functional principle, the work done by the pressure is equal to the stored energy of rods [19]; that is, the work done by the external force is equal to the sum of the strain energy of the rod system in numerical value, which is expressed as follows:(8)12FΔl=U1+U2,where U1 is the sum of strain energy in tension and compression and U2 is the sum of strain energy in bending. It shows that the total deformation of the structure unit can be obtained by calculating the strain energy of each rod. The specific analysis is as follows.Figure8(a) is the force analysis of the rod. There is no slip in the horizontal direction on node O for its resultant force in the horizontal direction is 0, so that node O displaces only in the vertical direction. In order to learn the strain energy conveniently, the node O is set as the fixed point, and the rods OA and OB are taken as the research objects for the force and deformation analysis.Figure 8
Force analysis and profile simplification of rods.
(a)(b)The resultant force of OA rod is an axial force ofF, which leads to tension and compression deformation. For OB rod, the force FB is decomposed into axial force and vertical force to obtain the forces N and P, respectively, and the rod is subject to axial force and bending moment M, resulting in tension and compression deformation and bending deformation. According to the symmetry of the structure and the uniformity of the force, the force condition of rods OC and OD is the same as that of OB. Therefore, tension and compression strain energy is produced in all the rods of the deformed unit, and bending strain energy is produced in the rods OB, OC, and OD.The equations of strain energy in tension and compression and in bending are, respectively, as follows:(9)U1=∫0lF2x2EAxdx,(10)U2=∫0lM2x2EIdx,where E represents the inherent elastic modulus of the raw material.F (x) in equation (9) is a function of the force on the rod, which is a constant function in this example. Ax is a function of the cross-sectional area of the rod. Due to the complexity of the change of the cross-sectional dimension of the designed member, as shown in Figure 8(b), the profile of the rod is simplified from the change of curved surface to the change of wedge shape for the convenience of calculation. In this research, (11)Ax=π6−1r−3−2x2,0≤x≤2r,πr2,2r<x<l−2r,π6−1r−3−2l−x2,l−2r≤x≤l.In equation (10), Mx is a function of bending moment, which of OB rod is(12)Mx=Px−l2.The value ofP is FBcos 19.47°. I is the section moment of inertia of the rod, that of cylinder section is I0=πr4/4; for the wedge section, the reduced moment of inertia of the variable section rod Iz is calculated as follows [20]:(13)Iz=161/I1+2/I2+9/I3+4/I4=πr41.3159,where I1, I2, I3, and I4 are the section moments of inertia when x=32/4r, x=2/2r, x=2/4r, and x=0, respectively, as shown in Figure 8(b).In conclusion, the total tensile and compressive strain energy of the structure unit is as follows:(14)U1=U1OA+3U1OB=∫0lF22EAxdx+3∫0lN22EAxdx.The value ofN is FBsin 19.47°. The total bending strain energy of the rods in the structure unit is expressed as follows:(15)U2=3U2OB=3∫0lM2x2EIdx=3∫02rM2x2EIzdx+∫2rl−2rM2x2EI0dx+∫l−2rlM2x2EIzdx.From the above calculation process, it can be observed that rod OA mainly bears load by tension and compression deformation, while rods OB, OC, and OD have very small tension and compression deformation and mainly bear load by bending deformation.From the above equations (8) ∼ (15), the total deformation of the structure can be obtained. That is,(16)Δl=σπE0.0196∅5r4−0.1826∅4r3+2.0220∅3r2−2.1427∅2r.Therefore, the equivalent elastic modulusE∗ of diamond lattice structure is obtained by equation (17)(17)E∗=σε=σΔhΔl,where ε represents the strain produced by the stress of the structure cell and h is the original height of the cell, whose value is 6/3∅. After sorting out, the following result is obtained:(18)E∗=E0.0077∅/r4−0.0712∅/r3+0.7883∅/r2−0.8353∅/r.Obviously, the equivalent elastic modulus of diamond lattice structure is related to the elastic modulusE of the raw material itself and the parameter ratio of structure cell ∅/r; the parameter ratio here is the inverse of that described in Section 2. The value of E/E∗ is a quartic function of the parameter ratio ∅/r, while from the above study, it is concluded that the relative density is a cubic function of the parameter ratio ∅/r; it is shown that the equivalent elastic modulus E∗ has a nonlinear relationship with the relative density ρx. Therefore, when the parameter ratio is fixed, the relative density and the equivalent elastic modulus of diamond lattice structure are fixed, respectively. After calculation, the following relationship is obtained:(19)E∗=E0.4239ρx3/2−0.0008.


## 3.2. Simulation of Mechanical Properties of Diamond Lattice Structure Unit

For lattice structure, the stress-strain curve changes with the increase of relative density, and the corresponding elastic modulus, plastic yield stress, and ultimate stress of brittle fracture increase [21]. With the increase of the parameter ratio r/∅, the density of diamond structure increases, and the mechanical properties change.A structure unit is intercepted to simulate the upper unit transferring three equal forcesF along the direction of the rods, and the stress state of the unit under different relative density is observed. When the maximum stress of the structure reaches the yield stress, the value of F reaches the upper limit, which is recorded as Fmax. Aluminum alloy 1061 was selected as the simulation material, with yield stress of 27.57 MPa and elastic modulus of 69 GPa. Figure 9 shows the mechanical simulation analysis of diamond lattice structure units under different rod section radius (0.5 mm, 1.0 mm, 1.5 mm, and 2 mm) when the inner tangent diameter of HP cell section ∅ is 10 mm. It can be seen from the von Mises stress that the most dangerous point of the structure occurs at the joint between the node and the rod.Figure 9
Simulation analysis of diamond lattice structure unit.As seen in Table2 which shows the comparison of simulation result data, the larger the section radius r of the rod, the greater the relative density of the structure, the greater the load-bearing capacity, and the smaller the deformation of the structure; all these changes are nonlinear as shown in Figure 10.Table 2
Comparison of simulation data of diamond lattice structure unit.
Serial numberRod section radiusr (mm)Relative densityρxForceFmax (N)Maximum stressσs(MPa)Maximum displacements (mm) (×10−3)Simulation equivalent modulus of elasticityE1∗ (GPa)Theoretical equivalent modulus of elasticityE2∗ (GPa)10.50.02653.64427.575.6020.06170.071020.60.03756.57327.574.9130.12610.157230.70.050710.70027.574.4070.22890.278740.80.065715.90227.573.9710.37760.437450.90.082524.31927.573.8110.60160.637961.00.101132.34627.573.6100.84480.885071.10.121442.69727.573.4401.17021.182081.20.143356.20127.573.3711.57181.531591.30.166970.79727.573.2642.04501.9391101.40.192188.59227.573.2182.59562.4075111.50.2188110.53927.573.2313.22562.9383121.60.2469130.45027.573.1123.95213.5331131.70.2765160.65027.573.1734.77354.1974141.80.3075190.28027.573.1475.70064.9323151.90.3398219.68127.573.0636.76195.73841620.3734250.63627.572.9667.96716.6186Figure 10
Variation trend of parameters of structure unit with section radius of rod.Comparing the equivalent elastic modulusE1∗ calculated by the simulation data with the equivalent elastic modulus E2∗ calculated by the above theoretical research, as shown in Figure 11, the overall trend of the curve is consistent.Figure 11
Comparison of theoretical and simulation results of equivalent modulus of elasticity.However, with the increase of the parameter ratio, the bending moment condition of OB, OC, and OD rods in the model changes into the shear condition. Therefore, the actual equivalent elastic modulus of the lattice structure is higher than the calculated value of the theoretical model in the region of high parameter ratio, as shown in Figure11. The existence of this gap makes it necessary to adjust the theoretical model for the high parameter ratio in future work.


## 4. Variable Density Design of Diamond Lattice Structure

### 4.1. Variable Density Hypothesis and Simulation Verification

When the structure is stressed, the internal stress decreases with the increase of the distance from the stressed point. According to the above research, the load-bearing capacity of the structure increases with the increase of the relative density of the structure.In the observation of organisms in nature, their tissues are often distributed in areas with large stress and deformation, in order to keep the original shape and physiological function of organism under external load without damage. Some organisms can dynamically adjust the flow direction of internal tissue to the area according to the specific position of external load, so as to reduce the influence of stress and deformation on themselves.Therefore, according to the tissue distribution law of natural organisms and the above research conclusions, it is assumed that the variable density structure has better bearing capacity; that is, with the distance from the stress point, the structure density is smaller and smaller, and such structure combination is better than the equal density structure.According to the hypothesis, the diamond structure is designed with variable density by changing the section radius distribution of the rods filled with lattice cell on the premise that the total material volume and total relative density of the structure remain unchanged. The constant density structure cell is adjusted by variable density, and the stress distribution and displacement are simulated and compared under the same stress conditions; the calculation results shown in Figures12 and 13 are obtained, respectively.Figure 12
Comparison of stress state. (a) The superposition of constant density units (r = 1 mm). (b) The superposition of variable density units (r = 1.2 mm and 0.8 mm).
(a)(b)Figure 13
Comparison of displacement state. (a) The superposition of constant density units (r = 1 mm). (b) The superposition of variable density units (r = 1.2 mm and 0.8 mm).
(a)(b)From Figure12, it can be seen intuitively that for the equal volume structure of materials under the same stress condition, the superposition of internal variable density units makes the structure stress distribution more uniform than that of constant density units. The maximum stress is also reduced from 23.65 MPa to 18.34 MPa, a decrease of 22.45%. As seen from Figure 13, the maximum displacement, which occurs in the stress center of the structure, is reduced from 3.320 × 10−3 to 2.587 × 10−3 mm, a decrease of 22.08%. It shows that the structure adjusted by variable density according to stress conditions can bear more force and the stress distribution of the whole structure is more uniform and reasonable. The above hypothesis is verified.


### 4.2. Additive Manufacturing and Compression Contrast Test of Diamond Lattice Structure

AlSi10 Mg aluminum alloy powder is used to make diamond lattice structure by SLS technology. The mechanical properties of the material are as follows: modulus of elasticity 70 ± 5 GPa, yield strength 245 ± 10 MPa, and forming accuracy of 0.05 mm; the formed structure is shown in Figure14. From left to right, there are traditional octahedral lattice structure, constant density diamond lattice structure, and variable density diamond lattice structure, and they are named specimen 1, specimen 2, and specimen 3, respectively. Various variables of the three kinds of specimens, including total volume, size, relative density of lattice part, and thickness of upper and lower bottom surfaces, are designed to be equal. A stress boss with a thickness of 2 mm is designed at the center of the top of the structure to simulate the local stress of the structure.Figure 14
Specimens with AM and their compression test.The compression tests are carried out on the universal mechanical testing machine, and the force-displacement curves shown in Figure15 are finally obtained.Figure 15
Comparison of force-displacement curves of three kinds of specimens.The test curves are analyzed as follows:First of all, the three kinds of specimens have gone through the elastic stage, yield stage, and strengthening stage in the compression process. Their curves basically coincide in the elastic stage and yield stage because they are of the same relative density, and they all enter into a short yield stage when the pressure is about 13 kN. Secondly, the three kinds of specimens show different trends in the strengthening stage: the displacement of specimen 1 changes faster and faster with the increase of force until the bottom edge of the specimen breaks; specimen 2 and specimen 3 show higher compressive capacity in the strengthening stage, and the displacement basically increases at a constant speed with the increase of force, but small cracks appear due to the problem of the forming quality of the specimen, resulting in the trend of the force-displacement curve being step like in the later stage until reaching the failure condition. Specimens 2 and 3 do not form obvious fracture on the structure after failure.Finally, the failure loads of the three kinds of specimens are 21.1 kN, 28.1 kN, and 38.8 kN, respectively. It can be seen through comparison that the compressive strength of specimen 3 is higher than that of specimen 2.Table3 is the analysis and comparison of the compression test results. The average equivalent modulus of elasticity of the diamond lattice structure specimens obtained from the test results is 4.165, and the equivalent modulus of elasticity calculated from equation (19) is 4.084, the two are relatively consistent, which preliminarily proves the correctness of the theoretical research.Table 3
Analysis of test results of three kinds of specimens.
NameSpecimen 1Specimen 2Specimen 3Quality (g)313131Overall size (mm)∅30 h40.66∅30 h40.66∅30 h40.66Relative density0.2690.2690.269Yield displacement (mm)0.33650.32900.3330Yield load (kN)13.0412.9113.13Failure load (kN)21.128.138.8Failure time (s)356.634478.112684.657Equivalent modulus of elasticity by formula (GPa)—4.0844.084Equivalent modulus of elasticity by test (GPa)4.1034.1554.175According to the comparison results of Figure15 and Table 3, the following test conclusions can be obtained:(i)
Compared with the octahedral lattice structure; the diamond lattice structure shows better compressive properties in the strengthening stage.(ii)
The compressive strength of the variable density lattice structure adapted to local stress is higher than that of the uniform one.(iii)
The formula of equivalent modulus of elasticity studied in Section2 has certain applicability.


## 4.1. Variable Density Hypothesis and Simulation Verification

When the structure is stressed, the internal stress decreases with the increase of the distance from the stressed point. According to the above research, the load-bearing capacity of the structure increases with the increase of the relative density of the structure.In the observation of organisms in nature, their tissues are often distributed in areas with large stress and deformation, in order to keep the original shape and physiological function of organism under external load without damage. Some organisms can dynamically adjust the flow direction of internal tissue to the area according to the specific position of external load, so as to reduce the influence of stress and deformation on themselves.Therefore, according to the tissue distribution law of natural organisms and the above research conclusions, it is assumed that the variable density structure has better bearing capacity; that is, with the distance from the stress point, the structure density is smaller and smaller, and such structure combination is better than the equal density structure.According to the hypothesis, the diamond structure is designed with variable density by changing the section radius distribution of the rods filled with lattice cell on the premise that the total material volume and total relative density of the structure remain unchanged. The constant density structure cell is adjusted by variable density, and the stress distribution and displacement are simulated and compared under the same stress conditions; the calculation results shown in Figures12 and 13 are obtained, respectively.Figure 12
Comparison of stress state. (a) The superposition of constant density units (r = 1 mm). (b) The superposition of variable density units (r = 1.2 mm and 0.8 mm).
(a)(b)Figure 13
Comparison of displacement state. (a) The superposition of constant density units (r = 1 mm). (b) The superposition of variable density units (r = 1.2 mm and 0.8 mm).
(a)(b)From Figure12, it can be seen intuitively that for the equal volume structure of materials under the same stress condition, the superposition of internal variable density units makes the structure stress distribution more uniform than that of constant density units. The maximum stress is also reduced from 23.65 MPa to 18.34 MPa, a decrease of 22.45%. As seen from Figure 13, the maximum displacement, which occurs in the stress center of the structure, is reduced from 3.320 × 10−3 to 2.587 × 10−3 mm, a decrease of 22.08%. It shows that the structure adjusted by variable density according to stress conditions can bear more force and the stress distribution of the whole structure is more uniform and reasonable. The above hypothesis is verified.


## 4.2. Additive Manufacturing and Compression Contrast Test of Diamond Lattice Structure

AlSi10 Mg aluminum alloy powder is used to make diamond lattice structure by SLS technology. The mechanical properties of the material are as follows: modulus of elasticity 70 ± 5 GPa, yield strength 245 ± 10 MPa, and forming accuracy of 0.05 mm; the formed structure is shown in Figure14. From left to right, there are traditional octahedral lattice structure, constant density diamond lattice structure, and variable density diamond lattice structure, and they are named specimen 1, specimen 2, and specimen 3, respectively. Various variables of the three kinds of specimens, including total volume, size, relative density of lattice part, and thickness of upper and lower bottom surfaces, are designed to be equal. A stress boss with a thickness of 2 mm is designed at the center of the top of the structure to simulate the local stress of the structure.Figure 14
Specimens with AM and their compression test.The compression tests are carried out on the universal mechanical testing machine, and the force-displacement curves shown in Figure15 are finally obtained.Figure 15
Comparison of force-displacement curves of three kinds of specimens.The test curves are analyzed as follows:First of all, the three kinds of specimens have gone through the elastic stage, yield stage, and strengthening stage in the compression process. Their curves basically coincide in the elastic stage and yield stage because they are of the same relative density, and they all enter into a short yield stage when the pressure is about 13 kN. Secondly, the three kinds of specimens show different trends in the strengthening stage: the displacement of specimen 1 changes faster and faster with the increase of force until the bottom edge of the specimen breaks; specimen 2 and specimen 3 show higher compressive capacity in the strengthening stage, and the displacement basically increases at a constant speed with the increase of force, but small cracks appear due to the problem of the forming quality of the specimen, resulting in the trend of the force-displacement curve being step like in the later stage until reaching the failure condition. Specimens 2 and 3 do not form obvious fracture on the structure after failure.Finally, the failure loads of the three kinds of specimens are 21.1 kN, 28.1 kN, and 38.8 kN, respectively. It can be seen through comparison that the compressive strength of specimen 3 is higher than that of specimen 2.Table3 is the analysis and comparison of the compression test results. The average equivalent modulus of elasticity of the diamond lattice structure specimens obtained from the test results is 4.165, and the equivalent modulus of elasticity calculated from equation (19) is 4.084, the two are relatively consistent, which preliminarily proves the correctness of the theoretical research.Table 3
Analysis of test results of three kinds of specimens.
NameSpecimen 1Specimen 2Specimen 3Quality (g)313131Overall size (mm)∅30 h40.66∅30 h40.66∅30 h40.66Relative density0.2690.2690.269Yield displacement (mm)0.33650.32900.3330Yield load (kN)13.0412.9113.13Failure load (kN)21.128.138.8Failure time (s)356.634478.112684.657Equivalent modulus of elasticity by formula (GPa)—4.0844.084Equivalent modulus of elasticity by test (GPa)4.1034.1554.175According to the comparison results of Figure15 and Table 3, the following test conclusions can be obtained:(i)
Compared with the octahedral lattice structure; the diamond lattice structure shows better compressive properties in the strengthening stage.(ii)
The compressive strength of the variable density lattice structure adapted to local stress is higher than that of the uniform one.(iii)
The formula of equivalent modulus of elasticity studied in Section2 has certain applicability.


## 5. Conclusion

The aim of this paper is to study the theoretical effect of relative density on the mechanical properties of the diamond lattice structure and to use this study to design a varying density lattice structure optimized for equivalent elastic modulus in compression based on additive manufacturing.A highly detailed geometric model is developed to determine the exact volume of the diamond lattice within a greater cell superposition, thereby allowing an accurate determination of the relative density of the lattice structure given the rod diameter and diameter of inscribed circle representing the size of the unit fill space. The mathematical model of the relative density and mechanical properties of the unit is established, and the relationship between them is obtained. A series of simulations are carried out to verify them, and the results are relatively consistent. The data of compression test in Section4 also confirms the correctness of the mathematical model accuracy.The hypothesis of variable density is put forward and verified by simulation and experiment. The results show that the diamond lattice structure has more advantages than the traditional octahedral structure in compression, and the variable density diamond lattice structure whose density of the filling unit gradually decreasing with the stress point as the center has better compression performance.



---

*Source: 8816593-2021-01-16.xml*