Mechanical Properties of Optimized Diamond Lattice Structure for Bone Scaffolds Fabricated via Selective Laser Melting

Developments in selective laser melting (SLM) have enabled the fabrication of periodic cellular lattice structures characterized by suitable properties matching the bone tissue well and by fluid permeability from interconnected structures. These multifunctional performances are significantly affected by cell topology and constitutive properties of applied materials. In this respect, a diamond unit cell was designed in particular volume fractions corresponding to the host bone tissue and optimized with a smooth surface at nodes leading to fewer stress concentrations. There were 33 porous titanium samples with different volume fractions, from 1.28 to 18.6%, manufactured using SLM. All of them were performed under compressive load to determine the deformation and failure mechanisms, accompanied by an in-situ approach using digital image correlation (DIC) to reveal stress–strain evolution. The results showed that lattice structures manufactured by SLM exhibited comparable properties to those of trabecular bone, avoiding the effects of stress-shielding and increasing longevity of implants. The curvature of optimized surface can play a role in regulating the relationship between density and mechanical properties. Owing to the release of stress concentration from optimized surface, the failure mechanism of porous titanium has been changed from the pattern of bottom-up collapse by layer (or cell row) to that of the diagonal (45°) shear band, resulting in the significant enhancement of the structural strength.


Introduction
Recently, cellular lattice structures featuring multifunctional performances [1] including high strength, lightweight and good energy absorption have been extensively studied as a suitable candidate for biomedical applications such as osseointegration or bone grafting [2,3] due to their inner topological complexity and comparable properties to the host bone. In such applications, a biomaterial should be biocompatible, not only possessing similar structural and mechanical properties to that of the bone it replaces, especially Young's modulus, but also retaining biological activities for tissue ingrowth and optimal osseointegration [4,5].
Lattice structures characterized by non-stochastic open unit cells have better mechanical properties in comparison to stochastic foams that exhibit localized deformations from internal imperfections [5][6][7]. Their controllable morphological parameters, such as pore architecture, pore size and volume fraction, can tailor the biomechanical properties matching the host tissue well, and the permeability

Design of Cellular Lattice Structures
The diamond lattice structures studied in this paper were generated by means of a CAD-based approach. Figure 1 shows the design method and process of the diamond lattice structure. As Figure 1d shows, a diamond unit cell that has fourteen nodes and sixteen equal edges is an isotropic geometry [4]. For a unit cell, the length of each strut between two nodes (L, and size of the unit cell (C) are related to each other by using the formula L = √ 3 4 c. The angle (θ) between strut and the horizontal plane equals 35 • 16 , which is a safe degree for avoiding deformation during the SLM process, resulting in a feature of being self-supporting. Table 1 shows the geometrical parameters of the diamond lattice structures studied in this article. designed by smooth optimized surfaces at unit nodes with several parameters accurately controlling porosity and pore size. Secondly, Ti6Al4V alloy samples with volume fractions ranging from 1.28 to 18.6% were fabricated using SLM system, followed by uniaxial compression experiments to measure the constitutive properties of SLM Ti6Al4V lattice structures, including compressive strength and elastic modulus to match the levels of human bones. Finally, the performance of diamond structures were demonstrated to be modified with surface optimization, resulting from the stress concentration reduction and the difference in the failure mechanism.

Design of Cellular Lattice Structures
The diamond lattice structures studied in this paper were generated by means of a CAD-based approach. Figure 1 shows the design method and process of the diamond lattice structure. As Figure  1d shows, a diamond unit cell that has fourteen nodes and sixteen equal edges is an isotropic geometry [4]. For a unit cell, the length of each strut between two nodes (L, and size of the unit cell (C) are related to each other by using the formula L = √3 4 . The angle (θ) between strut and the horizontal plane equals 35°16′, which is a safe degree for avoiding deformation during the SLM process, resulting in a feature of being self-supporting. Table 1 shows the geometrical parameters of the diamond lattice structures studied in this article. In the unit cell architecture, each optimized node is defined by rod diameter (D) and optimized radius (R). Combined with unit cell size, parameters including C, D, R are introduced as the basic elements that can accurately control pore size and volume fraction of this interconnected network. Compared with the straight beam-like nodes (Figure 1c), optimized nodes have smooth surfaces whose curvature is regulated by the value of optimized radius (Figure 1e), which is considered to In the unit cell architecture, each optimized node is defined by rod diameter (D) and optimized radius (R). Combined with unit cell size, parameters including C, D, R are introduced as the basic elements that can accurately control pore size and volume fraction of this interconnected network. Compared with the straight beam-like nodes (Figure 1c), optimized nodes have smooth surfaces whose curvature is regulated by the value of optimized radius (Figure 1e), which is considered to own better manufacturability for SLM. Then, they are assembled into the diamond lattice (Figure 1d), forming the finished unit cell architecture.
A final architecture can be carried out by a periodic array of the unit cell along with three mutually perpendicular directions. Figure 1a,b depict diamond-type 3 × 3 × 3 cellular structures, sample 0600 with straight beam-like struts and sample 0616 with optimized surface at its unit nodes, respectively, where '06' denotes D = 0.6 mm, '16' denotes R = 1.6 mm, and other numbers (for example '1008') are similar. In this paper, to better study the performance of diamond structures, the value of the unit cell size is set to be 5.5 mm to provide the ability for SLM manufacturing of cellular lattice structures with a wide range of volume fractions [19]. Dimensions of the test samples are fixed at 16.5 × 16.5 × 16.5 mm 3 , because 3 × 3 × 3 cell types are a compromise, combining economy with accuracy, as Smith et al. [33] studied cell size effect and stated that many properties could be well predicted with lower amounts of unit cells.  Based on this design approach, different combinations of pore sizes, volume fraction values and mechanical properties can be obtained for specific application requirements by modifying the parameters above. The volume fraction, as one of the key factors controlling the mechanical properties of porous parts, is determined by rod diameter, optimized radius and unit cell size. As shown in Figure 2, when the unit size is set as a constant, increasing the rod diameter is an effective way to achieve the rapid growth of volume fraction, in order to get the desired mechanical properties, while an optimized radius that contributes little to the increase of volume fraction may play a role to improve and adjust the stiffness and strength of diamond structures. own better manufacturability for SLM. Then, they are assembled into the diamond lattice (Figure 1d), forming the finished unit cell architecture. A final architecture can be carried out by a periodic array of the unit cell along with three mutually perpendicular directions. Figure 1a,b depict diamond-type 3 × 3 × 3 cellular structures, sample 0600 with straight beam-like struts and sample 0616 with optimized surface at its unit nodes, respectively, where '06' denotes D = 0.6 mm, '16' denotes R = 1.6 mm, and other numbers (for example '1008') are similar. In this paper, to better study the performance of diamond structures, the value of the unit cell size is set to be 5.5 mm to provide the ability for SLM manufacturing of cellular lattice structures with a wide range of volume fractions [19]. Dimensions of the test samples are fixed at 16.5 × 16.5 × 16.5 mm 3 , because 3 × 3 × 3 cell types are a compromise, combining economy with accuracy, as Smith et al. [33] studied cell size effect and stated that many properties could be well predicted with lower amounts of unit cells. Table 1. Geometrical parameters of the diamond lattice structures.
Based on this design approach, different combinations of pore sizes, volume fraction values and mechanical properties can be obtained for specific application requirements by modifying the parameters above. The volume fraction, as one of the key factors controlling the mechanical properties of porous parts, is determined by rod diameter, optimized radius and unit cell size. As shown in Figure 2, when the unit size is set as a constant, increasing the rod diameter is an effective way to achieve the rapid growth of volume fraction, in order to get the desired mechanical properties, while an optimized radius that contributes little to the increase of volume fraction may play a role to improve and adjust the stiffness and strength of diamond structures.

Materials
Titanium alloy is an important metal for medical devices because of its excellent osseointegration, superior corrosion resistance and favorable mechanical properties. Thus, cellular lattice structures were made from Ti6Al4V powder that was produced by a gas atomization process. As shown in Figure 3, SEM images of Ti6Al4V powder exhibit a feature of a nearly spherical shape and smooth surfaces, indicating a good flowability. The particle size of the powder is in the range of 3-50 μm.

Materials
Titanium alloy is an important metal for medical devices because of its excellent osseointegration, superior corrosion resistance and favorable mechanical properties. Thus, cellular lattice structures were made from Ti6Al4V powder that was produced by a gas atomization process. As shown in Figure 3, SEM images of Ti6Al4V powder exhibit a feature of a nearly spherical shape and smooth surfaces, indicating a good flowability. The particle size of the powder is in the range of 3-50 µm. The chemical composition is listed in Table 2. The powder has a low content of impurities (oxygen, hydrogen and nitrogen), which might avoid negative effect of embrittlement [34]. The chemical composition is listed in Table 2. The powder has a low content of impurities (oxygen, hydrogen and nitrogen), which might avoid negative effect of embrittlement [34].

Fabrication by SLM
Cellular lattice structures with rod diameter from 0.2 to 1.0 mm and optimized radius from 0 to 1.6 mm were manufactured by a SLM machine (EOSINT-M280, EOS GmbH, Munich, Germany) equipped with a 200 W CW ytterbium fiber laser. The diameter of the laser beam focused on the powder bed is 0.1 mm. The laser melting process occurs in a protective argon atmosphere with O 2 content less than 0.1 vol.%, and the processing parameters were set up as follows: the laser power was 170 W; the track length of the laser beam was 5 mm; the hatch spacing was 75 μm; the scanning speed was 1250 mm/s; the beam offset was 0.015 mm; and the layer thickness was 30 μm. Using these parameters, the density of bulk Ti6Al4V was measured as 4.42 g/cm 3 (≥ 99.5%) by the Archimedes method, tensile strength was 930 MPa and modulus of elasticity was 116 GPa after a post processing of heat treatment. There were 33 samples in a uniform dimension of 16.5 × 16.5 × 16.5 mm 3 built in one batch to reduce deviation during fabrication. When the SLM process was finished, the proposed samples were removed from the base plate via wire electrical discharge machining (wire-EDM); Figure 4 shows the samples for compressive test. Figure 5 exhibits the comparison between CAD models and as-built samples that was observed under an optical microscope, revealing a good consistency. Even the sample 0200 with the volume fraction of 1.28% (theoretical value) and rod diameter of 0.2 mm, which is considered to be the most difficult to form, has been manufactured successfully with no defects or broken cells, confirming the capability of SLM to fabricate these cellular lattice structures. For the compressive test, the as-built samples were subjected to a post processing of heat treatment (annealing for 3 h at 650 °C) in an argon gas atmosphere with subsequent furnace cooling for thermal stress releasing.

Fabrication by SLM
Cellular lattice structures with rod diameter from 0.2 to 1.0 mm and optimized radius from 0 to 1.6 mm were manufactured by a SLM machine (EOSINT-M280, EOS GmbH, Munich, Germany) equipped with a 200 W CW ytterbium fiber laser. The diameter of the laser beam focused on the powder bed is 0.1 mm. The laser melting process occurs in a protective argon atmosphere with O 2 content less than 0.1 vol.%, and the processing parameters were set up as follows: the laser power was 170 W; the track length of the laser beam was 5 mm; the hatch spacing was 75 µm; the scanning speed was 1250 mm/s; the beam offset was 0.015 mm; and the layer thickness was 30 µm. Using these parameters, the density of bulk Ti6Al4V was measured as 4.42 g/cm 3 (≥ 99.5%) by the Archimedes method, tensile strength was 930 MPa and modulus of elasticity was 116 GPa after a post processing of heat treatment. There were 33 samples in a uniform dimension of 16.5 × 16.5 × 16.5 mm 3 built in one batch to reduce deviation during fabrication. When the SLM process was finished, the proposed samples were removed from the base plate via wire electrical discharge machining (wire-EDM); Figure 4 shows the samples for compressive test. Figure 5 exhibits the comparison between CAD models and as-built samples that was observed under an optical microscope, revealing a good consistency. Even the sample 0200 with the volume fraction of 1.28% (theoretical value) and rod diameter of 0.2 mm, which is considered to be the most difficult to form, has been manufactured successfully with no defects or broken cells, confirming the capability of SLM to fabricate these cellular lattice structures. For the compressive test, the as-built samples were subjected to a post processing of heat treatment (annealing for 3 h at 650 • C) in an argon gas atmosphere with subsequent furnace cooling for thermal stress releasing.

Test Details
Uniaxial compressive tests of the lattice structures were conducted by using a universal testing machine equipped with a 100 KN load cell. Most of the samples subjected to load were parallel to the building direction, except for those samples signed "XY" whose direction of loading was parallel to the layer. The compressive tests were carried out at a displacement rate of 2 mm/min in a continuous process without intermittence, accompanied by a video camera with a resolution of 10 megapixels and a frame rate of 29 fps for continuously recording the deformation of each sample. The camera was equipped with a macro lens and a spotlight, ensuring the proper illumination necessary for digital image correlation (DIC) analyses. Then, frames correlated with strain were extracted from these videos, providing information about the failure behavior of the lattice structures. In order to compare the uniaxial test results, the strain was calculated by dividing the displacement by sample height of 16.5 mm, while the stress was computed from the load and apparent cross-sectional area of 16.5 × 16.5 mm 2 . Elastic modulus (E) was considered as the slope of linear fits of the stress-strain curve. And stress at the first peak of the experimental data was regarded as ultimate stress (σ).

Test Details
Uniaxial compressive tests of the lattice structures were conducted by using a universal testing machine equipped with a 100 KN load cell. Most of the samples subjected to load were parallel to the building direction, except for those samples signed "XY" whose direction of loading was parallel to the layer. The compressive tests were carried out at a displacement rate of 2 mm/min in a continuous

Test Details
Uniaxial compressive tests of the lattice structures were conducted by using a universal testing machine equipped with a 100 KN load cell. Most of the samples subjected to load were parallel to the building direction, except for those samples signed "XY" whose direction of loading was parallel to the layer. The compressive tests were carried out at a displacement rate of 2 mm/min in a continuous

Analysis of Mechanical Properties
In general, comparisons for the same topology of micro-lattice are made according to the volume fraction or relative density which is calculated by dividing the density of micro-lattice block (ρ) by that of parent material (ρ 0 ). So the volume fraction (Vf ), porosity (P), were calculated as: In addition, in order to evaluate the effect of optimized radius on mechanical properties of diamond lattice structure, relative modulus (E * ) and relative ultimate stress (σ * ) are adopted to assess its material efficiency, namely: Stress-strain curves of samples with rod diameter of 0.4, 0.6, 0.8, 1 (mm) are shown in Figure 6a-d, respectively. These stress-strain responses are characterized by a nearly linear region followed by yield at a plastic region that maintained constant stress in a longer strain range due to the influence of plasticity enhancement with the increase of the rod radius. The corresponding σ * and E * values with optimized radiuses from 0.6 mm to 1.6 mm are depicted in Figure 6e-h. It is observed that the relative ultimate stress and relative modulus of the samples are distinctly linearly increasing with the increase of optimized radius, whatever the sizes of the rod diameter are, indicating that the surface optimization at nodes can improve the mechanical properties of diamond lattice structure and enhance the performance efficiency of this cellular materials. The introduction of surface optimization was considered to produce more uniform transmission of external force, leading to the reduction of stress concentration and the change of stress distribution on each strut cell. Therefore, it can be speculated that the mechanical properties of other lattice structures can be further improved by taking measure such as, as a simple way, the surface optimization at nodes.
On the other hand, an important concept to note is the anisotropy. By comparing the data obtained from different loading direction shown in Table 2, it is obvious that both the elastic modulus and ultimate strength of samples perpendicular to the building direction (signed "XY") are higher than the parallel ones. The anisotropy of mechanical properties that are commonly found in SLM parts can be attributed to the oriented growth process along the direction of substrate layer-by-layer [35].
It is well known that the higher a volume fraction is generated, the better mechanical properties of porous materials. In order to derive a predictable relationship between mechanical properties and volume fraction, a significant theoretical model of Gibson-Ashby was used as following [10,36]: where E 0 , ρ 0 and σ 0 are the elastic modulus, density and strength of fully dense solid materials, and E, ρ and σ are the apparent modulus, density and strength of studied lattice samples, respectively. The exponents n 1 and n 2 are the density factors. According to the previous studies, n 1 and n 2 were calculated to be 2 and 1.5, respectively, by theoretical computational models [16] and by experimental approach [14,36]. The Gibson-Ashby prefactor C 1 and C 2 are influenced by the pore shape and size, the value of these could be changed by the introduction of optimized radius. Using experimental data can yield the corresponding value, expressed as Materials 2018, 11, x FOR PEER REVIEW 8 of 17 Figure 6. (a-d) The stress-strain curves; (e-h) the corresponding relative ultimate stress * and relative modulus * .
It is well known that the higher a volume fraction is generated, the better mechanical properties of porous materials. In order to derive a predictable relationship between mechanical properties and volume fraction, a significant theoretical model of Gibson-Ashby was used as following [10,36]: According to the experimental values in Table 3, the corresponding constant C 1 and C 2 values can be calculated by the Equations (6) and (7), as shown in Figure 7. It is observed that, with the increase of the optimized radius, C 1 keeps almost constant or has increased by a small margin, while C 2 has increased significantly, which indicates that the surface optimization at the nodes contributes more to the improvement of strength than elastic modulus. In comparison, although relative elastic modulus and the relative strength in previous section have been improved simultaneously. While in the view of the Gibson-Ashby model, only strength is enhanced in fact. That is to say, there is a difference in influence of optimized surface at lattice nodes on their elastic modulus and the compressive strength, which may provide a way to adjusting the relationship between them by varying the curvature of optimized surface. pore shape and size, the value of these could be changed by the introduction of optimized radius. Using experimental data can yield the corresponding value, expressed as According to the experimental values in Table 3, the corresponding constant 1 and 2 values can be calculated by the Equations (6) and (7), as shown in Figure 7. It is observed that, with the increase of the optimized radius, 1 keeps almost constant or has increased by a small margin, while      Table 3 lists the performance data of all tested samples. The diamond cellular lattice structures studied in this paper have an interconnected high porosity of 81-97% that is considered suitable for tissue ingrowth and vascularization [13]. These three-dimensional structures with high interconnected porosity provide biological anchorage for the ingrowth of mineralized tissue into the pore space, enhancing the fixation of the implant in the surrounding bone as well as improving its long-term stability [10]. As Figure 2 shows a wide range of the porosity with controllability, the studied lattice structures provide great promise for creating biological alternatives for implants, in which the interconnected porosity plays a vital role in cell seeding, nutrients transfer to stimulate bone formation.

Bone Modulus Matching
Some of the literatures cataloged the properties of cortical bone and trabecular bone from multiple locations within the human skeleton, as summarized in Table 4. One of the most striking features of bone's properties is the huge variation in modulus and strength [41]. Comparison between Tables 3 and 4 shows that the elastic modulus of diamond structures and human bone are in a similar range, which is a significant conclusion for the application of the lattice titanium structures built by SLM as biomaterials for bone replacement. As a consequence, the diamond structure reveals properties comparable with trabecular bone, which can minimize the stress-shielding effect and increase the longevity of implants. What needs to be pointed out is that the listed values in Table 4 should only represent the dimensions of bone's properties, because the properties of bone depend on numerous factors, such as anatomical site, the age and sex of the bone tissue, methods of storage and testing conditions [13]. Thus, more requirements should be given to scaffold for replacing natural tissues, for example, by constructing gradient or heterogeneous structures for satisfying multi-scale structures in human bone. The performance of diamond lattice shown in Table 3 is expected to provide a reference in the future study. Figure 8a,b show plots of ultimate stress versus volume fraction values and the elastic modulus versus volume fraction values, respectively. It is clear that strength and modulus increase linearly with the increase of density, and the linear coefficients of optimized surface structures are higher than the ones of non-optimization. It was indicated that the optimized radius played a role in regulating the relationship between density and mechanical properties including the strength and modulus. This regulation controlled by the value of optimized radius may extend the range of mechanical properties of the diamond structure without changing volume fraction and permeability. In addition, the smooth surface inner complex shape contributes to removing the unmelted or semi-melted powders and enhancing surface finish by post-processing such as magnetic fluid-assisted polishing and electrochemical polishing, as these micro powders on the implant surface could trigger inflammation.  Figure 8a,b show plots of ultimate stress versus volume fraction values and the elastic modulus versus volume fraction values, respectively. It is clear that strength and modulus increase linearly with the increase of density, and the linear coefficients of optimized surface structures are higher than the ones of non-optimization. It was indicated that the optimized radius played a role in regulating the relationship between density and mechanical properties including the strength and modulus. This regulation controlled by the value of optimized radius may extend the range of mechanical properties of the diamond structure without changing volume fraction and permeability. In addition, the smooth surface inner complex shape contributes to removing the unmelted or semimelted powders and enhancing surface finish by post-processing such as magnetic fluid-assisted polishing and electrochemical polishing, as these micro powders on the implant surface could trigger inflammation. A linear relationship was found between the elastic modulus and strength in Figure 8c, which is similar to the properties of bone stated in the previous article [41]. Figure 8c also shows that the A linear relationship was found between the elastic modulus and strength in Figure 8c, which is similar to the properties of bone stated in the previous article [41]. Figure 8c also shows that the compression strength of the optimized structures is superior to that of non-optimized structures under the same condition of elastic modulus, increasing fatigue endurance for long-term application in the human body. Moreover, it is possible to accurately predict and even tailor the mechanical properties of the diamond lattices by governing the volume fraction and adjusting the curvature of node surface, because the relationship between the mechanical properties and volume fraction has been exactly established.

Deformation Behavior
The compressive force-deformation behaviors of 0600, 0610 and 0616 samples are plotted in Figure 9. All of them experienced four distinct regions: a region of elastic deformation up to the first peak value, followed by a collapse in which most of the strength is lost, then a region of fluctuation in stress, and a densification region characterized by a rapid increase in stress, as shown in Figure 9. These regions are in good agreement with that observed in the previous investigations [14,42].
properties of the diamond lattices by governing the volume fraction and adjusting the curvature of node surface, because the relationship between the mechanical properties and volume fraction has been exactly established.

Deformation Behavior
The compressive force-deformation behaviors of 0600, 0610 and 0616 samples are plotted in Figure 9. All of them experienced four distinct regions: a region of elastic deformation up to the first peak value, followed by a collapse in which most of the strength is lost, then a region of fluctuation in stress, and a densification region characterized by a rapid increase in stress, as shown in Figure 9. These regions are in good agreement with that observed in the previous investigations [14,42]. The stress-strain curve of the 0600 sample shows a significant plateau of stress in the fluctuant region, indicating a typical characteristic of brittle deformation behavior [5]. In comparison, with the increase of the optimized radius, the stiffness and strength have been greatly enhanced, whereas the amplitude in the fluctuant region becomes higher gradually. It should be noted that, in the region of collapse, all of these curves bore high maximum stress before a steep drop that occurred at strains of 8.74%, 8.78%, and 7.71% respectively. The bigger optimized radius was used, the steeper and deeper stress drop occurred, resulting in a long-term low stress in 0616. These differences, on one hand, can be attributed to the effect of the increase of volume fraction from optimized surface at nodes. On the other hand, the more important reason is that the distribution of stress in the unit structure has been changed by the smooth surface, resulting in the release of stress concentration of unit cell and the increase of stability of the lattice structure. Moreover, strain-stress curves show that higher load bearing capacity can be expected when stress concentration has been released, leading to the mechanism of the shearing crack for sample 0616, as shown in Table 5. Table 5 shows the video frames of samples in the complete stress-strain process. Observing their deformation behavior, each fracture occurred as a random event in a non-uniform way on the local scale of each sample, but the overall deformation behavior was determined by the structure itself and The stress-strain curve of the 0600 sample shows a significant plateau of stress in the fluctuant region, indicating a typical characteristic of brittle deformation behavior [5]. In comparison, with the increase of the optimized radius, the stiffness and strength have been greatly enhanced, whereas the amplitude in the fluctuant region becomes higher gradually. It should be noted that, in the region of collapse, all of these curves bore high maximum stress before a steep drop that occurred at strains of 8.74%, 8.78%, and 7.71% respectively. The bigger optimized radius was used, the steeper and deeper stress drop occurred, resulting in a long-term low stress in 0616. These differences, on one hand, can be attributed to the effect of the increase of volume fraction from optimized surface at nodes. On the other hand, the more important reason is that the distribution of stress in the unit structure has been changed by the smooth surface, resulting in the release of stress concentration of unit cell and the increase of stability of the lattice structure. Moreover, strain-stress curves show that higher load bearing capacity can be expected when stress concentration has been released, leading to the mechanism of the shearing crack for sample 0616, as shown in Table 5. Table 5 shows the video frames of samples in the complete stress-strain process. Observing their deformation behavior, each fracture occurred as a random event in a non-uniform way on the local scale of each sample, but the overall deformation behavior was determined by the structure itself and parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns.
Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45 • ) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45 • ) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. parent material. In addition, the change of stress distribution has a significant influence on the crushing behavior of these structures. The impact of the control parameter of optimized radius on the mechanical properties can be revealed for different deformation patterns. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve. Failure mechanism of sample 0600 without optimized surface at nodes was accompanied by a layer-by-layer deformation from the bottom up. Due to the effect of the successive layers of collapse, 0600 has given a periodic stress fluctuation curve with a small amplitude until the densification region, indicating that the stretch dominated deformation behavior played a key role in the deformation mechanism of such structures. Furthermore, fracture occurred almost exclusively at the lattice nodes, and then rod-shaped struts separated from junction nodes were observed during the compression process owing to the buckling at the junction. SEM morphologies of this fracture surfaces shown in Figure 10 provide a fracture mechanism of tensile yield with ductile dimpling.
The in-situ images of 0610 show that, when the strain went to 10-15%, a diagonal (45°) shear band was observed; at the same time, the lowest unit cells were also destroyed. Deformation was accompanied by a layer-by-layer buckling and a diagonal (45°) shear band simultaneously, leading to a slight fluctuation in the lowest stress region, which is similar to the secondary peak in the 0600 stress-strain curve.  Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616. Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45 • ) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45 • ) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a long-term process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern. Table 6. Failure mechanism comparison with previous article.

Structures Optimization Volume Fraction, (%) Deformation Modes Deformation Pattern
Materials 2018, 11, x FOR PEER REVIEW 14 of 17 Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.  Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.  Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.

Non 22
Materials 2018, 11, x FOR PEER REVIEW 14 of 17 Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.  Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.  Figure 10. SEM morphologies of fracture surface: (a-c) tensile yield at node with ductile dimpling for sample 0600, (d-g) shear rupture at rod for sample 0616.
In contrast, testing sample 0616, destruction was mainly presented in the pattern of a diagonal (45°) shear band, whereas upward-layer failure pattern was not observed. The SEM micrographs of sample 0616 in Figure 10 provide evidence that the fracture micro-mechanisms are quite different from that of 0600. They show a typical shear rupture along the slip plane, ensuring the assumption mentioned above that the loading stress has been transferred from the node to the rod center under the influence of optimized surface. This transference gives rise in the reduction of stress concentration, resulting in the improvement of mechanical properties.
Nevertheless, effected by the collective failure on the diagonal (45°) shear band, the curve of 0616 showed a more severe collapse, in which 95% of the structure's strength was lost, following a longterm process of low-stress and a secondary linear elasticity behavior. The occurrence of this secondary linear elasticity is principally because other structures outside the shear band are involved in the deformation after the contacts of the struts in the shear band. As shown in Table 6, failure mechanism of sample 0616 was in good agreement with diamond structures as in previous article [5], but their optimum conditions were different. Both sample 0616 and structure of [5] showed the same deformation pattern of bending dominated under different condition in which non-optimized structure had a volume fraction of 22%, while optimized one had 8.97%. Similarly, an entirely different deformation pattern was found in the same structure, the main difference between 0600 and [5] is the volume fraction. Comparison in Table 6 indicates that the deformation pattern is determined not only by the unit cell type consisting of porous structures [43] but also by the volume fraction of the structure and the stress distribution of the nodes. According to Deshpande et al. [44], the deformation behavior of the diamond structure is bending dominated, but the sample 0600 exhibits a tendency of stretch dominated in the pattern of a layer-by-layer failure, which demonstrated the effect of volume fraction on deformation pattern.

Conclusions
In this study, a diamond unit cell with an optimized surface at nodes was designed in order to improve and adjust its mechanical properties. The diamond cellular lattice structures with a wide range of volume fractions from 1.28% to 18.6%, having a unit cell size (5.5 mm), were manufactured by the SLM using the Ti6Al4V alloy, showing the great flexibility of this design approach and its suitability for matching the properties of porous human bone. The effects of the node optimization on the mechanical properties and deformation behavior under compressive test were investigated. The major findings of this research are summarized as follows.

1.
Mechanical properties including strength and modulus increase with the increase of optimized radius at a fixed volume fraction, indicating that the stress concentration at lattice nodes have been released by the introduction of surface optimization; 2.
A minitype database that correlates porosity levels to mechanical properties has been established. The compressive strength and elastic modulus of diamond lattice structures are in the range of 2-78 MPa and 34-1403 MPa, respectively, which are comparable with trabecular bone; 3.
By optimizing the nodes with different curvature in each unit cell, the relationship between compressive strength and elastic modulus can be regulated, enabling the cellular structures more suitable for bone implants implication; 4.
The change of stress distribution stemming from the optimized surface at nodes has a significant influence on the deformation behavior of these structures. With the analysis of in-situ images and SEM, it is revealed that the improvement of mechanical properties was mainly due to the transition of the deformation behavior from upward-layer collapse in non-optimization structures to diagonal (45 • ) shear band in optimized ones.
Although diamond structures fabricated by SLM show a good potential in terms of mechanical tailoring for bone scaffolds, there are some defects to be resolved at the present stage. It is a time-consuming and tedious manual process for CAD-based lattice structures, especially for more complex functionally graded scaffolds. In order to better mimic host tissue, scaffolds with hybridizing structures are needed to satisfy the biological requirements for the better regeneration of bone tissue. The database of mechanical properties studied in this article can be a reference for bone scaffold applications. Besides this, effected by the limitations in the diameter of the laser beam and powder particle size, the geometric feature size of the studied scaffolds is still too large to match that of human bone. With the development of SLM manufacturing in the future work, more accurate parts to be fabricated will be desirable. There still is a long distance to application due to some other biological factors which are not strictly addressed in this work. Other measures such as spraying hydroxyapatite to improve the biological activity and testing the permeability are supposed to be taken in the future work.