## Efficient and Accurate Simulation of Antenna Arrays in Ansys HFSS

### Unit-cell in HFSS

HFSS offers different method of creating and simulating a large array. The explicit method, shown in Figure 1(a) might be the first method that comes to our mind. This is where you create the exact CAD of the array and solve it. While this is the most accurate method of simulating an array, it is computationally extensive. This method may be non-feasible for the initial design of a large array. The use of unit cell (Figure 1(b)) and array theory helps us to start with an estimate of the array performance by a few assumptions. Finite Array Domain Decomposition (or FADDM) takes advantage of unit cell simplicity and creates a full model using the meshing information generated in a unit cell. In this blog we will review the creation of unit cell. In the next blog we will explain how a unit cell can be used to simulate a large array and FADDM.

In a unit cell, the following assumptions are made:

• The pattern of each element is identical.
• The array is uniformly excited in amplitude, but not necessarily in phase.
• Edge affects and mutual coupling are ignored

A unit cell works based on Master/Slave (or Primary/Secondary) boundary around the cell. Master/Slave boundaries are always paired. In a rectangular cell you may use the new Lattice Pair boundary that is introduced in Ansys HFSS 2020R1. These boundaries are means of simulating an infinite array and estimating the performance of a relatively large arrays. The use of unit cell reduces the required RAM and solve time.

Primary/Secondary (Master/Slave) (or P/S) boundaries can be combined with Floquet port, radiation or PML boundary to be used in an infinite array or large array setting, as shown in Figure 3.

To create a unit cell with P/S boundary, first start with a single element with the exact dimensions of the cell. The next step is creating a vacuum or airbox around the cell. For this step, set the padding in the location of P/S boundary to zero. For example, Figure 4 shows a microstrip patch antenna that we intend to create a 2D array based on this model. The array is placed on the XY plane. An air box is created around the unit cell with zero padding in X and Y directions.

You notice that in this example the vacuum box is larger than usual size of quarter wavelength that is usually used in creating a vacuum region around the antenna. We will get to calculation of this size in a bit, for now let’s just assign a value or parameter to it, as it will be determined later. The next step is to define P/S to generate the lattice. In AEDT 2020R1 this boundary is under “Coupled” boundary. There are two methods to create P/S: (1) Lattice Pair, (2) Primary/Secondary boundary.

Lattice Pair

The Lattice Pair works best for square lattices. It automatically assigns the primary and secondary boundaries. To assign a lattice pair boundary select the two sides that are supposed to create infinite periodic cells, right-click->Assign Boundary->Coupled->Lattice Pair, choose a name and enter the scan angles. Note that scan angles can be assigned as parameters. This feature that is introduced in 2020R1 does not require the user to define the UV directions, they are automatically assigned.

### Primary/Secondary

Primary/Secondary boundary is the same as what used to be called Master/Slave boundary. In this case, each Secondary (Slave) boundary should be assigned following a Primary (Master) boundary UV directions. First choose the side of the cell that Primary boundary. Right-click->Assign Boundary->Coupled->Primary. In Primary Boundary window define U vector. Next select the secondary wall, right-click->Assign Boundary->Couple->Secondary, choose the Primary Boundary and define U vector exactly in the same direction as the Primary, add the scan angles (the same as Primary scan angles)

### Floquet Port and Modes Calculator

Floquet port excites and terminates waves propagating down the unit cell. They are similar to waveguide modes. Floquet port is always linked to P/S boundaries. Set of TE and TM modes travel inside the cell. However, keep in mind that the number of modes that are absorbed by the Floquet port are determined by the user. All the other modes are short-circuited back into the model. To assign a Floquet port two major steps should be taken:

Defining Floquet Port

Select the face of the cell that you like to assign the Floquet port. This is determined by the location of P/S boundary. The lattice vectors A and B directions are defined by the direction of lattice (Figure 7).

The number of modes to be included are defined with the help of Modes Calculator. In the Mode Setup tab of the Floquet Port window, choose a high number of modes (e.g. 20) and click on Modes Calculator. The Mode Table Calculator will request your input of Frequency and Scan Angles. After selecting those, a table of modes and their attenuation using dB/length units are created. This is your guide in selecting the height of the unit cell and vaccume box. The attenation multiplied by the height of the unit cell (in the project units, defined in Modeler->Units) should be large enough to make sure the modes are attenuated enough so removing them from the calcuatlion does not cause errors. If the unit cell is too short, then you will see many modes are not attenuated enough. The product of the attenuatin and height of the airbox should be at least 50 dB. After the correct size for the airbox is calcualted and entered, the model with high attenuation can be removed from the Floquet port definition.

The 3D Refinement tab is used to control the inclusion of the modes in the 3D refinement of the mesh. It is recommended not to select them for the antenna arrays.

In our example, Figure 8 shows that the 5th mode has an attenuation of 2.59dB/length. The height of the airbox is around 19.5mm, providing 19.5mm*2.59dB/mm=50.505dB attenuation for the 5th mode. Therefore, only the first 4 modes are kept for the calculations. If the height of the airbox was less than 19.5mm, we would need to increase the height so accordingly for an attenuation of at least 50dB.

A simpler alternative for Floquet port is radiation boundary. It is important to note that the size of the airbox should still be kept around the same size that was calculated for the Floquet port, therefore, higher order modes sufficiently attenuated. In this case the traditional quarter wavelength padding might not be adequate.

Perfectly Matched Layer

Although using radiation boundary is much simpler than Floquet port, it is not accurate for large scan angles. It can be a good alternative to Floquet port only if the beam scanning is limited to small angles. Another alternative to Floquet port is to cover the cell by a layer of PML. This is a good compromise and provides very similar results to Floquet port models. However, the P/S boundary need to surround the PML layer as well, which means a few additional steps are required. Here is how you can do it:

1. Reduce the size of the airbox* slightly, so after adding the PML layer, the unit cell height is the same as the one that was generated using the Modes Calculation. (For example, in our model airbox height was 19mm+substrte thickness, the PML height was 3mm, so we reduced the airbox height to 16mm).
2. Choose the top face and add PML boundary.
3. Select each side of the airbox and create an object from that face (Figure 10).
4. Select each side of the PML and create objects from those faces (Figure 10).
5. Select the two faces that are on the same plane from the faces created from airbox and PML and unite them to create a side wall (Figure 10).
6. Then assign P/S boundary to each pair of walls (Figure 10).

*Please note for this method, an auto-size “region” cannot be used, instead draw a box for air/vacuum box. The region does not let you create the faces you need to combine with PML faces.

The advantage of PML termination over Floquet port is that it is simpler and sometimes faster calculation. The advantage over Radiation Boundary termination is that it provides accurate results for large scan angles. For better accuracy the mesh for the PML region can be defined as length based.

Seed the Mesh

To improve the accuracy of the PML model further, an option is to use length-based mesh. To do this select the PML box, from the project tree in Project Manager window right-click on Mesh->Assign Mesh Operation->On Selection->Length Based. Select a length smaller than lambda/10.

Scanning the Angle

In phased array simulation, we are mostly interested in the performance of the unit cell and array at different scan angles. To add the scanning option, the phase of P/S boundary should be defined by project or design parameters. The parameters can be used to run a parametric sweep, like the one shown in Figure 12. In this example the theta angle is scanned from 0 to 60 degrees.

Comparing PML and Floquet Port with Radiation Boundary

To see the accuracy of the radiation boundary vs. PML and Floquet Port, I ran the simulations for scan angles up to 60 degrees for a single element patch antenna. Figure 13 shows that the accuracy of the Radiation boundary drops after around 15 degrees scanning. However, PML and Floquet port show similar performance.

S Parameters

To compare the accuracy, we can also check the S parameters. Figure 14 shows the comparison of active S at port 1 for PML and Floquet port models. Active S parameters were used since the unit cell antenna has two ports. Figure 15 shows how S parameters compare for the model with the radiation boundary and the one with the Floquet port.

### Conclusion

The unit cell definition and options on terminating the cell were discussed here. Stay tuned. In the next blog we discuss how the unit cell is utilized in modeling antenna arrays.

## Using Material Designer To Perform Homogenization Studies

#### Editor’s Note:

3D Printing and other advanced manufacturing methods are driving the increased use of lattice-type structures in structural designs. This is great for reducing mass and increasing the stiffness of components but can be a real pain for those of us doing simulation. Modeling all of those tiny features across a part is difficult to mesh and takes forever to solve.

PADT has been doing a bit of R&D in this area recently, including a recent PHASE II NASA STTR with ASU and KSU. We see a lot of potential in combining generative design and 3D Printing to drive better structures. The key to this effort is efficient and accurate simulation.

The good news is that we do not have to model every unit cell. Instead, we can do some simulation on a single representative chunk and use the ANSYS Material Designer feature to create an approximate material property that we can use to represent the lattice volume as a homogeneous material.

In the post below, PADT’s Alex Grishin explains it all with theory, examples, and a clear step-by-step process that you can use for your lattice filled geometry.

## How to Use Lattice Optimization in ANSYS Mechanical and ANSYS SpaceClaim 19.2

One of the great new features in ANSYS Mechanical 19.2 is the ability to perform a lattice optimization.  Accomplished as an option within Topology Optimization, lattice optimization allows us to generate a lattice structure within our region of interest.  It includes varying thickness of the lattice members as part of the optimization.

Lattice structures can be very beneficial because weight can be substantially reduced compared to solid parts made using traditional manufacturing methods.  Further, recent advances in additive manufacturing enable the creation of lattice structures in ways that weren’t possible with traditional manufacturing.

Here I’ll explain how to perform a lattice optimization in ANSYS 19.2 step by step.

The procedure starts the same as a normal topology optimization in ANSYS Mechanical, with an initial static structural analysis on our original part or assembly.  If you’re not familiar with the process, this earlier PADT Focus blog should be helpful:  http://www.padtinc.com/blog/the-focus/topological-optimization-in-ansys-18-1-motorcycle-component-example

For the lattice optimization, I’m starting with a part I created that acts as a corner brace:

At this early point in the simulation, the Project Schematic looks like this:

I used the Multizone mesh method to get a hex mesh on the part:

Simple loads and constraints are recommended especially if you’ll be doing a downstream validation study.  That is because the downstream simulation on the resulting lattice geometry will most likely need to operate on the FE entities rather than geometric entities for load and constraint application. The boundary conditions in this simple model consisted of a fixed support on one side of the brace and a force load on the other side:

After solving, I reviewed the displacement as well as the stress results:

Satisfied with the results, the next step is to add a Topology Optimization block in the Project Schematic. The easiest way to do this is to right click on the Solution cell, then select Transfer Data to New > Topology Optimization:

You may need to re-solve the static structural simulation at this point.  You’ll know if you have yellow thunderbolts in the Project Schematic instead of green checkmarks for the Static Structural analysis.

At this point, the Project Schematic now looks like this:

The Mechanical window now has the Topology Optimization branch added:

The change to make to enable a lattice optimization is accomplished in the details view of the Optimization Region branch:

We then need to specify some settings for the lattice.  The first of these is the Lattice Type.  The various types are documented in the ANSYS 19.2 Help.  In my example I selected the Crossed option.

The other properties to define are:

• Minimum Density (to avoid lattice structures that are toothin.  Allowed bounds are 0 and 1)
• Maximum Density (elements are considered full/solid fordensities higher than this value, allowed bounds are 0 and 1)
• Lattice Cell Size (used in downstream geometry steps andadditive manufacturing)

Values I used in my example are shown here:

Assuming no other options need to be set, we solve the lattice optimization and review the results.  The results are displayed as a contour plot with values between zero and one, with values corresponding to the density settings as specified above.

Note that at this stage we don’t actually visualize the lattice structure – just a contour plot of where the lattice can be in the structure.  Where density values are higher than the maximum density specified, the geometry will end up being solid.  The lattice structure can exist where the results are between the minimum and maximum density values specified, with a varying thickness of lattice members corresponding to higher and lower densities.

The next step is to bring the lattice density information into SpaceClaim and generate actual lattice geometry.  This is done by adding a free standing Geometry block in the Workbench Project Schematic.

The next step is to drag and drop the Results cell from the Topology Optimization block onto the Geometry cell of the new free standing Geometry block:

The Project Schematic will now look like this:

Notice the Results cell in the Topology Optimization branch now has a yellow lightning bolt.  The next step is to right click on that Results cell and Update.  The Project Schematic will now look like this:

Before we can open SpaceClaim, we next need to right click on the Geometry cell in the downstream Geometry block and Update that as well:

After both Updates, the Project Schematic will now look like this:

The next step is to double click or right click on the now-updated Geometry cell to open SpaceClaim.  Note that both the original geometry and a faceted version of the geometry will exist in SpaceClaim:

It may seem counter intuitive, but we actually suppress the faceted geometry and only work with the original, solid geometry for the faceted process.  The faceted geometry should be automatically suppressed, as shown by the null symbol, ø, in the SpaceClaim tree.  At this point it will be helpful to hide the faceted geometry by unchecking its box in the tree:

Next we’ll utilize some capability in the Facets menu in SpaceClaim to create the lattice geometry, using the lattice distribution calculated by the lattice optimization.  Click on the Facets tab, then click on the Shell button:

Set the Infill option to be Basic:

At this point there should be a check box for “Use Density Attributes” below the word Shape.  This check box doesn’t always appear.  If it’s not there, first try clicking on the actual geometry object in the tree:

In one instance I had to go to %appdata%\Ansys and rename the v192 folder to v192.old to reset Workbench preferences and launch Workbench again.  That may have been ‘pilot error’ on my part as I was learning the process.

The next step is to check the Use density attributes box.  The Shape dropdown should be set to Lattices.  Once the Use density attributes box is checked, we can then one of the predefined lattice shapes, which will be used for downstream simulation and 3D printing.  The shape picked needs to match the lattice shape previously picked in the topology optimization.

In my case I selected the Cube Lattice with Side Diagonal Supports, which corresponds to the Crossed selection I made in the upsteam lattice optimization.  Note that a planar preview of this is displayed inside the geometry:

The next step is to click the green checkmark to have SpaceClaim create the lattice geometry based on the lattice distribution calculated by the lattice optimization:

When SpaceClaim is done with the lattice geometry generation, you should be able to see a ghosted image showing the lattice structure in the part’s interior:

Note that if you change views, etc., in SpaceClaim, you may then see the exterior surfaces of the part, but rest assured the lattice structure remains in the interior.

Your next step may need to be a validation.  To do this, we create a standalone Static Structural analysis block on the Project Schematic:

Next we drag and drop the Geometry cell from the faceted geometry block we just created onto the Geometry cell of the newly created Static Structural block:

We can now open Mechanical for the new Static Structural analysis.  Note that the geometry that comes into Mechanical in this manner will have a single face for the exterior, and a single face for the exterior. To verify that the lattice structure is actually in the geometry, I recommend creating a section plane so we can view the interior of the geometry:

To mesh the lattice structure, I’ve found that inserting a Mesh Method and setting it to the Tetrahedrons/Patch Independent option has worked for getting a reasonable mesh.  Care must be taken with element sizes or a very large mesh will be created.  My example mesh has about 500,000 nodes.  This is a section view, showing the mesh of the interior lattice structure (relatively coarse for the example).

For boundary condition application, I used Direct FE loads.  I used a lasso pick after aligned the view properly to select the nodes needed for the displacement and then the force loads, and created Named Selections for each of those nodal selections for easy load application.

Here are a couple of results plots showing a section view with the lattice in the interior (deflection followed by max principal stress):

Here is a variant on the lattice specifications, in which the variance in the thickness of the lattice members (a result of the optimization) is more evident:

Clearly, a lot more could be done with the geometry in SpaceClaim before a validation step or 3D printing.  However, hopefully this step by step guide is helpful with the basic process for performing a lattice optimization in ANSYS Mechanical and SpaceClaim 19.2.

## On the Functions of Cellular Structures in Nature

WHY did nature evolve cellular structures?

In a previous post, I laid out a structural classification of cellular structures in nature, proposing that they fall into 6 categories. I argued that it is not always apparent to a designer what the best unit cell choice for a given application is. While most mechanical engineers have a feel for what structure to use for high stiffness or energy absorption, we cannot easily address multi-objective problems or apply these to complex geometries with spatially varying requirements (and therefore locally optimum cellular designs). However, nature is full of examples where cellular structures possess multi-objective functionality: bone is one such well-known example. To be able to assign structure to a specific function requires us to connect the two, and to do that, we must identify all the functions in play. In this post, I attempt to do just that and develop a classification of the functions of cellular structures.

Any discussion of structure in nature has to contend with a range of drivers and constraints that are typically not part of an engineer’s concern. In my discussions with biologists (including my biochemist wife), I quickly run into justified skepticism about whether generalized models associating structure and function can address the diversity and nuance in nature – and I (tend to) agree. However, my attempt here is not to be biologically accurate – it is merely to construct something that is useful and relevant enough for an engineer to use in design. But we must begin with a few caveats to ensure our assessments consider the correct biological context.

## 1. Uniquely Biological Considerations

Before I attempt to propose a structure-function model, there are some legitimate concerns many have made in the literature that I wish to recap in the context of cellular structures. Three of these in particular are relevant to this discussion and I list them below.

### 1.1 Design for Growth

Engineers are familiar with “design for manufacturing” where design considers not just the final product but also aspects of its manufacturing, which often place constraints on said design. Nature’s “manufacturing” method involves (at the global level of structure), highly complex growth – these natural growth mechanisms have no parallel in most manufacturing processes. Take for example the flower stalk in Fig 1, which is from a Yucca tree that I found in a parking lot in Arizona.

At first glance, this looks like a good example of overlapping surfaces, one of the 6 categories of cellular structures I covered before. But when you pause for a moment and query the function of this packing of cells (WHY this shape, size, packing?), you realize there is a powerful growth motive for this design. A few weeks later when I returned to the parking lot, I found many of the Yucca stems simultaneously in various stages of bloom – and captured them in a collage shown in Fig 2. This is a staggering level of structural complexity, including integration with the environment (sunlight, temperature, pollinators) that is both wondrous and for an engineer, very humbling.

The lesson here is to recognize growth as a strong driver in every natural structure – the tricky part is determining when the design is constrained by growth as the primary force and when can growth be treated as incidental to achieving an optimum functional objective.

### 1.2 Multi-functionality

Even setting aside the growth driver mentioned previously, structure in nature is often serving multiple functions at once – and this is true of cellular structures as well. Consider the tessellation of “scutes” on the alligator. If you were tasked with designing armor for a structure, you may be tempted to mimic the alligator skin as shown in Fig. 3.

As you begin to study the skin, you see it is comprised of multiple scutes that have varying shape, size and cross-sections – see Fig 4 for a close-up.

The pattern varies spatially, but you notice some trends: there exists a pattern on the top but it is different from the sides and the bottom (not pictured here). The only way to make sense of this variation is to ask what functions do these scutes serve? Luckily for us, biologists have given this a great deal of thought and it turns out there are several: bio-protection, thermoregulation, fluid loss mitigation and unrestricted mobility are some of the functions discussed in the literature [1, 2]. So whereas you were initially concerned only with protection (armor), the alligator seeks to accomplish much more – this means the designer either needs to de-confound the various functional aspects spatially and/or expand the search to other examples of natural armor to develop a common principle that emerges independent of multi-functionality specific to each species.

### 1.3 Sub-Optimal Design

This is an aspect for which I have not found an example in the field of cellular structures (yet), so I will borrow a well-known (and somewhat controversial) example [3] to make this point, and that has to do with the giraffe’s Recurrent Laryngeal Nerve (RLN), which connects the Vagus Nerve to the larynx as shown in Figure 5, which it is argued, takes an unnecessarily long circuitous route to connect these two points.

We know that from a design standpoint, this is sub-optimal because we have an axiom that states the shortest distance between two points is a straight line. And therefore, the long detour the RLN makes in the giraffe’s neck must have some other evolutionary and/or developmental basis (fish do not have this detour) [3]. However, in the case of other entities such as the cellular structures we are focusing on, the complexity of the underlying design principles makes it hard to identify cases where nature has found a sub-optimal design space for the function of interest to us, in favor of other pressing needs determined by selection. What is sufficient for the present moment is to appreciate that such cases may exist and to bear them in mind when studying structure in nature.

## 2. Classifying Functions

Given the above challenges, the engineer may well ask: why even consider natural form in making determinations involving the design of engineering structures? The biomimic responds by reminding us that nature has had 3.8 billion years to develop a “design guide” and we would be wise to learn from it. Importantly, natural and engineering structures both exist in the same environment and are subject to identical physics and further, are both often tasked with performing similar functions. In the context of cellular structures, we may thus ask: what are the functions of interest to engineers and designers that nature has addressed through cellular design? Through my reading [1-4], I have compiled the classification of functions in Figure 6, though this is likely to grow over time.

This broad classification into structural and transport may seem a little contrived, but it emerges from an analyst’s view of the world. There are two reasons why I propose this separation:

1. Structural functions involve the spatial allocation of materials in the construction of the cellular structures, while transport functions involve the structure AND some other entity and their interactions (fluid or light for example) – thus additional physics needs to be comprehended for transport functions
2. Secondly, structural performance needs to be comprehended independent of any transport function: a cellular structure must retain its integrity over the intended lifetime in addition to performing any additional function

Each of these functions is a fascinating case study in its own right and I highly recommend the site AskNature.org [1] as a way to learn more on a specific application, but this is beyond the scope of the current post. More relevant to our high-level discussion is that having listed the various reasons WHY cellular structures are found in nature, the next question is can we connect the structures described in the previous post to the functions tabulated above? This will be the attempt of my next post. Until then, as always, I welcome all inputs and comments, which you can send by messaging me on LinkedIn.

## References

2. Foy (1983), The grand design: Form and colour in animals, Prentice-Hall, 1st edition
3. Dawkins (2010), The greatest show on earth: the evidence for evolution, Free Press, Reprint of 1st edition
4. Gibson, Ashby, Harley (2010), Cellular Materials in Nature and Medicine, Cambridge University Press; 1st edition
5. Ashby, Evans, Fleck, Gibson, Hutchinson, Wadley (2000), Metal Foams: A Design Guide, Butterworth-Heinemann, 1st edition

## Cellular Design Strategies in Nature: A Classification

What types of cellular designs do we find in nature?

Cellular structures are an important area of research in Additive Manufacturing (AM), including work we are doing here at PADT. As I described in a previous blog post, the research landscape can be broadly classified into four categories: application, design, modeling and manufacturing. In the context of design, most of the work today is primarily driven by software that represent complex cellular structures efficiently as well as analysis tools that enable optimization of these structures in response to environmental conditions and some desired objective. In most of these software, the designer is given a choice of selecting a specific unit cell to construct the entity being designed. However, it is not always apparent what the best unit cell choice is, and this is where I think a biomimetic approach can add much value. As with most biomimetic approaches, the first step is to frame a question and observe nature as a student. And the first question I asked is the one described at the start of this post: what types of cellular designs do we find in the natural world around us? In this post, I summarize my findings.

## Design Strategies

In a previous post, I classified cellular structures into 4 categories. However, this only addressed “volumetric” structures where the objective of the cellular structure is to fill three-dimensional space. Since then, I have decided to frame things a bit differently based on my studies of cellular structures in nature and the mechanics around these structures. First is the need to allow for the discretization of surfaces as well: nature does this often (animal armor or the wings of a dragonfly, for example). Secondly, a simple but important distinction from a modeling standpoint is whether the cellular structure in question uses beam- or shell-type elements in its construction (or a combination of the two). This has led me to expand my 4 categories into 6, which I now present in Figure 1 below.

Setting aside the “why” of these structures for a future post, here I wish to only present these 6 strategies from a structural design standpoint.

1. Volumetric – Beam: These are cellular structures that fill space predominantly with beam-like elements. Two sub-categories may be further defined:
• Honeycomb: Honeycombs are prismatic, 2-dimensional cellular designs extruded in the 3rd dimension, like the well-known hexagonal honeycomb shown in Fig 1. All cross-sections through the 3rd dimension are thus identical. Though the hexagonal honeycomb is most well known, the term applies to all designs that have this prismatic property, including square and triangular honeycombs.
• Lattice and Open Cell Foam: Freeing up the prismatic requirement on the honeycomb brings us to a fully 3-dimensional lattice or open-cell foam. Lattice designs tend to embody higher stiffness levels while open cell foams enable energy absorption, which is why these may be further separated, as I have argued before. Nature tends to employ both strategies at different levels. One example of a predominantly lattice based strategy is the Venus flower basket sea sponge shown in Fig 1, trabecular bone is another example.
2. Volumetric – Shell:
• Closed Cell Foam: Closed cell foams are open-cell foams with enclosed cells. This typically involves a membrane like structure that may be of varying thickness from the strut-like structures. Plant sections often reveal a closed cell foam, such as the douglas fir wood structure shown in Fig 1.
• Periodic Surface: Periodic surfaces are fascinating mathematical structures that often have multiple orders of symmetry similar to crystalline groups (but on a macro-scale) that make them strong candidates for design of stiff engineering structures and for packing high surface areas in a given volume while promoting flow or exchange. In nature, these are less commonly observed, but seen for example in sea urchin skeletal plates.
3. Surface:
• Tessellation: Tessellation describes covering a surface with non-overlapping cells (as we do with tiles on a floor). Examples of tessellation in nature include the armored shells of several animals including the extinct glyptodon shown in Fig 1 and the pineapple and turtle shell shown in Fig 2 below.
• Overlapping Surface: Overlapping surfaces are a variation on tessellation where the cells are allowed to overlap (as we do with tiles on a roof). The most obvious example of this in nature is scales – including those of the pangolin shown in Fig 1.

This separation into 6 categories is driven from a designer’s and an analyst’s perspective – designers tend to think in volumes and surfaces and the analyst investigates how these are modeled (beam- and shell-elements are at the first level of classification used here). However, this is not sufficient since it ignores the function of the cellular design, which both designer and analyst need to also consider. In the case of tessellation on the skin of an alligator for example as shown in Fig 3, was it selected for protection, easy of motion or for controlling temperature and fluid loss?

In a future post, I will attempt to develop an approach to classifying cellular structures that derives not from its structure or mechanics as I have here, but from its function, with the ultimate goal of attempting to reconcile the two approaches. This is not a trivial undertaking since it involves de-confounding multiple functional requirements, accounting for growth (nature’s “design for manufacturing”) and unwrapping what is often termed as “evolutionary baggage,” where the optimum solution may have been sidestepped by natural selection in favor of other, more pressing needs. Despite these challenges, I believe some first-order themes can be discerned that can in turn be of use to the designer in selecting a particular design strategy for a specific application.

## References

This is by no means the first attempt at a classification of cellular structures in nature and while the specific 6 part separation proposed in this post was developed by me, it combines ideas from a lot of previous work, and three of the best that I strongly recommend as further reading on this subject are listed below.

1. Gibson, Ashby, Harley (2010), Cellular Materials in Nature and Medicine, Cambridge University Press; 1st edition
2. Naleway, Porter, McKittrick, Meyers (2015), Structural Design Elements in Biological Materials: Application to Bioinspiration. Advanced Materials, 27(37), 5455-5476
3. Pearce (1980), Structure in Nature is a Strategy for Design, The MIT Press; Reprint edition

As always, I welcome all inputs and comments – if you have an example that does not fit into any of the 6 categories mentioned above, please let me know by messaging me on LinkedIn and I shall include it in the discussion with due credit. Thanks!

## Modeling 3D Printed Cellular Structures: Approaches

How can the mechanical behavior of cellular structures (honeycombs, foams and lattices) be modeled?

This is the second in a two-part post on the modeling aspects of 3D printed cellular structures. If you haven’t already, please read the first part here, where I detail the challenges associated with modeling 3D printed cellular structures.

The literature on the 3D printing of cellular structures is vast, and growing. While the majority of the focus in this field is on the design and process aspects, there is a significant body of work on characterizing behavior for the purposes of developing analytical material models. I have found that these approaches fall into 3 different categories depending on the level of discretization at which the property is modeled: at the level of each material point, or at the level of the connecting member or finally, at the level of the cell. At the end of this article I have compiled some of the best references I could find for each of the 3 broad approaches.

## 1. Continuum Modeling

The most straightforward approach is to use bulk material properties to represent what is happening to the material at the cellular level [1-4]. This approach does away with the need for any cellular level characterization and in so doing, we do not have to worry about size or contact effects described in the previous post that are artifacts of having to characterize behavior at the cellular level. However, the assumption that the connecting struts/walls in a cellular structure behave the same way the bulk material does can particularly be erroneous for AM processes that can introduce significant size specific behavior and large anisotropy. It is important to keep in mind that factors that may not be significant at a bulk level (such as surface roughness, local microstructure or dimensional tolerances) can be very significant when the connecting member is under 1 mm thick, as is often the case.

The level of error introduced by a continuum assumption is likely to vary by process: processes like Fused Deposition Modeling (FDM) are already strongly anisotropic with highly geometry-specific meso-structures and an assumption like this will generate large errors as shown in Figure 1. On the other hand, it is possible that better results may be had for powder based fusion processes used for metal alloys, especially when the connecting members are large enough and the key property being solved for is mechanical stiffness (as opposed to fracture toughness or fatigue life).

## 2. Cell Level Homogenization

The most common approach in the literature is the use of homogenization – representing the effective property of the cellular structure without regard to the cellular geometry itself. This approach has significantly lower computational expense associated with its implementation. Additionally, it is relatively straightforward to develop a model by fitting a power law to experimental data [5-8] as shown in the equation below, relating the effective modulus E* to the bulk material property Es and their respective densities (ρ and ρs), by solving for the constants C and n.

While a homogenization approach is useful in generating comparative, qualitative data, it has some difficulties in being used as a reliable material model in analysis & simulation. This is first and foremost since the majority of the experiments do not consider size and contact effects. Secondly, even if these were considered, the homogenization of the cells only works for the specific cell in question (e.g. octet truss or hexagonal honeycomb) – so every new cell type needs to be re-characterized. Finally, the homogenization of these cells can lose insight into how structures behave in the transition region between different volume fractions, even if each cell type is calibrated at a range of volume fractions – this is likely to be exacerbated for failure modeling.

## 3. Member Modeling

The third approach involves describing behavior not at each material point or at the level of the cell, but at a level in-between: the connecting member (also referred to as strut or beam). This approach has been used by researchers [9-11] including us at PADT [12] by invoking beam theory to first describe what is happening at the level of the member and then use that information to build up to the level of the cells.

This approach, while promising, is beset with some challenges as well: it requires experimental characterization at the cellular level, which brings in the previously mentioned challenges. Additionally, from a computational standpoint, the validation of these models typically requires a modeling of the full cellular geometry, which can be prohibitively expensive. Finally, the theory involved in representing member level detail is more complex, makes assumptions of its own (e.g. modeling the “fixed” ends) and it is not proven adequately at this point if this is justified by a significant improvement in the model’s predictability compared to the above two approaches. This approach does have one significant promise: if we are able to accurately describe behavior at the level of a member, it is a first step towards a truly shape and size independent model that can bridge with ease between say, an octet truss and an auxetic structure, or different sizes of cells, as well as the transitions between them – thus enabling true freedom to the designer and analyst. It is for this reason that we are focusing on this approach.

## Conclusion

Continuum models are easy to implement and for relatively isotropic processes and materials such as metal fusion, may be a good approximation of stiffness and deformation behavior. We know through our own experience that these models perform very poorly when the process is anisotropic (such as FDM), even when the bulk constitutive model incorporates the anisotropy.

Homogenization at the level of the cell is an intuitive improvement and the experimental insights gained are invaluable – comparison between cell type performances, or dependencies on member thickness & cell size etc. are worthy data points. However, caution needs to be exercised when developing models from them for use in analysis (simulation), though the relative ease of their computational implementation is a very powerful argument for pursuing this line of work.

Finally, the member level approach, while beset with challenges of its own, is a promising direction forward since it attempts to address behavior at a level that incorporates process and geometric detail. The approach we have taken at PADT is in line with this approach, but specifically seeks to bridge the continuum and cell level models by using cellular structure response to extract a point-wise material property. Our preliminary work has shown promise for cells of similar sizes and ongoing work, funded by America Makes, is looking to expand this into a larger, non-empirical model that can span cell types. If this is an area of interest to you, please connect with me on LinkedIn for updates. If you have questions or comments, please email us at info@padtinc.com or drop me a message on LinkedIn.

## References (by Approach)

### Bulk Property Models

[1] C. Neff, N. Hopkinson, N.B. Crane, “Selective Laser Sintering of Diamond Lattice Structures: Experimental Results and FEA Model Comparison,” 2015 Solid Freeform Fabrication Symposium

[2] M. Jamshidinia, L. Wang, W. Tong, and R. Kovacevic. “The bio-compatible dental implant designed by using non-stochastic porosity produced by Electron Beam Melting®(EBM),” Journal of Materials Processing Technology214, no. 8 (2014): 1728-1739

[3] S. Park, D.W. Rosen, C.E. Duty, “Comparing Mechanical and Geometrical Properties of Lattice Structure Fabricated using Electron Beam Melting“, 2014 Solid Freeform Fabrication Symposium

[4] D.M. Correa, T. Klatt, S. Cortes, M. Haberman, D. Kovar, C. Seepersad, “Negative stiffness honeycombs for recoverable shock isolation,” Rapid Prototyping Journal, 2015, 21(2), pp.193-200.

### Cell Homogenization Models

[5] C. Yan, L. Hao, A. Hussein, P. Young, and D. Raymont. “Advanced lightweight 316L stainless steel cellular lattice structures fabricated via selective laser melting,” Materials & Design 55 (2014): 533-541.

[6] S. Didam, B. Eidel, A. Ohrndorf, H.‐J. Christ. “Mechanical Analysis of Metallic SLM‐Lattices on Small Scales: Finite Element Simulations versus Experiments,” PAMM 15.1 (2015): 189-190.

[7] P. Zhang, J. Toman, Y. Yu, E. Biyikli, M. Kirca, M. Chmielus, and A.C. To. “Efficient design-optimization of variable-density hexagonal cellular structure by additive manufacturing: theory and validation,” Journal of Manufacturing Science and Engineering 137, no. 2 (2015): 021004.

[8] M. Mazur, M. Leary, S. Sun, M. Vcelka, D. Shidid, M. Brandt. “Deformation and failure behaviour of Ti-6Al-4V lattice structures manufactured by selective laser melting (SLM),” The International Journal of Advanced Manufacturing Technology 84.5 (2016): 1391-1411.

### Beam Theory Models

[9] R. Gümrük, R.A.W. Mines, “Compressive behaviour of stainless steel micro-lattice structures,” International Journal of Mechanical Sciences 68 (2013): 125-139

[10] S. Ahmadi, G. Campoli, S. Amin Yavari, B. Sajadi, R. Wauthle, J. Schrooten, H. Weinans, A. Zadpoor, A. (2014), “Mechanical behavior of regular open-cell porous biomaterials made of diamond lattice unit cells,” Journal of the Mechanical Behavior of Biomedical Materials, 34, 106-115.

[11] S. Zhang, S. Dilip, L. Yang, H. Miyanji, B. Stucker, “Property Evaluation of Metal Cellular Strut Structures via Powder Bed Fusion AM,” 2015 Solid Freeform Fabrication Symposium

[12] D. Bhate, J. Van Soest, J. Reeher, D. Patel, D. Gibson, J. Gerbasi, and M. Finfrock, “A Validated Methodology for Predicting the Mechanical Behavior of ULTEM-9085 Honeycomb Structures Manufactured by Fused Deposition Modeling,” Proceedings of the 26th Annual International Solid Freeform Fabrication, 2016, pp. 2095-2106

## SFF Symposium 2016 Paper: Predicting the Mechanical Behavior of ULTEM-9085 Honeycomb Structures

Our work on  3D printed honeycomb modeling that started as a Capstone project with students from ASU in September 2015 (described in a previous blog post), was published in a peer-reviewed paper released last week in the proceedings of the SFF Symposium 2016. The full title of the paper is “A Validated Methodology for Predicting the Mechanical Behavior of ULTEM-9085 Honeycomb Structures Manufactured by Fused Deposition Modeling“. This was the precursor work that led to a us winning an 18-month award to pursue this work further with America Makes.

Abstract
ULTEM-9085 has established itself as the Additive Manufacturing (AM) polymer of choice for end-use applications such as ducts, housings, brackets and shrouds. The design freedom enabled by AM processes has allowed us to build structures with complex internal lattice structures to enhance part performance. While solutions exist for designing and manufacturing cellular structures, there are no reliable ways to predict their behavior that account for both the geometric and process complexity of these structures. In this work, we first show how the use of published values of elastic modulus for ULTEM-9085 honeycomb structures in FE simulation results in 40- 60% error in the predicted elastic response. We then develop a methodology that combines experimental, analytical and numerical techniques to predict elastic response within a 5% error. We believe our methodology is extendable to other processes, materials and geometries and discuss future work in this regard.

## Modeling 3D Printed Cellular Structures: Challenges

In this post, I discuss six challenges that make the modeling of 3D printed cellular structures (such as honeycombs and lattices) a non-trivial matter. In a following post, I will present how some of these problems have been addressed with different approaches.

At the outset, I need to clarify that by modeling I mean the analytical representation of material behavior, primarily for use in predictive analysis (simulation). Here are some reasons why this is a challenging endeavor for 3D printed cellular solids – some of these reasons are unique to 3D printing, others are a result of aspects that are specific to cellular solids, independent of how they are manufactured. I show examples with honeycombs since that is the majority of the work we have data for, but I expect that these ideas apply to foams and lattices as well, just with varying degrees of sensitivity.

## 1. Complex Geometry with Non-Uniform Local Conditions

I state the most well-appreciated challenge with cellular structures first: they are NOT fully-dense solid materials that have relatively predictable responses governed by straightforward analytical expressions. Consider a dogbone-shaped specimen of solid material under tension: it’s stress-strain response can be described fairly well using continuum expressions that do not account for geometrical features beyond the size of the dogbone (area and length for stress and strain computations respectively). However, as shown in Figure 1, such is not the case for cellular structures, where local stress and strain distributions are non-uniform. Further, they may have variable distributions of bending, stretching and shear in the connecting members that constitute the structure. So the first question becomes: how does one represent such complex geometry – both analytically and numerically?

## 2. Size Effects

A size effect is said to be significant when an observed behavior varies as a function of the size of the sample whose response is being characterized even after normalization (dividing force by area to get stress, for example). Here I limit myself to size effects that are purely a mathematical artifact of the cellular geometry itself, independent of the manufacturing process used to make them – in other words this effect would persist even if the material in the cellular structure was a mathematically precise, homogeneous and isotropic material.

It is common in the field of cellular structure modeling to extract an “effective” property – a property that represents a homogenized behavior without explicitly modeling the cellular detail. This is an elegant concept but introduces some practical challenges in implementation – inherent in the assumption is that this property, modulus for example, is equivalent to a continuum property valid at every material point. The reality is the extraction of this property is strongly dependent on the number of cells involved in the experimental characterization process. Consider experimental work done by us at PADT, and shown in Figure 2 below, where we varied both the number of axial and longitudinal cells (see inset for definition) when testing hexagonal honeycomb samples made of ULTEM-9085 with FDM. The predicted effective modulus increases with increasing number of cells in the axial direction, but reduces (at a lower rate) for increasing number of cells in the longitudinal direction.

This is a significant challenge and deserves a full form post to do justice (and is forthcoming), but the key to remember is that testing a particular cellular structure does not suffice in the extraction of effective properties. So the second question here becomes: what is the correct specimen design for characterizing cellular properties?

## 3. Contact Effects

In the compression test shown in the inset in Figure 2, there is physical contact between the platen and the specimen that creates a local effect at the top and bottom that is different from the experience of the cells closer the center. This is tied to the size effect discussed above – if you have large enough cells in the axial direction, the contribution of this effect should reduce – but I have called it out as a separate effect here for two reasons: Firstly, it raises the question of how best to design the interface for the specimen: should the top and bottom cells terminate in a flat plate, or should the cells extend to the surface of contact (the latter is the case in the above image). Secondly, it raises the question of how best to model the interface, especially if one is seeking to match simulation results to experimentally observed behavior. Both these ideas are shown in Figure 3 below. This also has implications for product design – how do we characterize and model the lattice-skin interface? As such, independent of addressing size effects, there is a need to account for contact behavior in characterization, modeling and analysis.

## 4. Macrostructure Effects

Another consideration related to specimen design is demonstrated in an exaggerated manner in the slowed down video below, showing a specimen flying off the platens under compression – the point being that for certain dimensions of the specimen being characterized (typically very tall aspect ratios), deformation in the macrostructure can influence what is perceived as cellular behavior. In the video below, there is some induced bending on a macro-level.

## 5. Dimensional Errors

While all manufacturing processes introduce some error in dimensional tolerances, the error can have a very significant effect for cellular structures – a typical industrial 3D printing process has tolerances within 75 microns (0.003″) – cellular structures (micro-lattices in particular) very often are 250-750 microns in thickness, meaning the tolerances on dimensional error can be in the 10% and higher error range for thickness of these members. This was our finding when working with Fused Deposition Modeling (FDM), where on a 0.006″ thick wall we saw about a 10% larger true measurement when we scanned the samples optically, as shown in Figure 4. Such large errors in thickness can yield a significant error in measured behavior such as elastic modulus, which often goes by some power to the thickness, amplifying the error. This drives the need for some independent measurement of the manufactured cellular structure – made challenging itself by the need to penetrate the structure for internal measurements. X-ray scanning is a popular, if expensive approach. But the modeler than has the challenge of devising an average thickness for analytical calculations and furthermore, the challenge of representation of geometry in simulation software for efficient analysis.

## 6. Mesostructural Effects

The layerwise nature of Additive Manufacturing introduces a set of challenges that are somewhat unique to 3D Printed parts. Chief among these is the resulting sensitivity to orientation, as shown for the laser-based powder bed fusion process in Figure 5 with standard materials and parameter sets. Overhang surfaces (unsupported) tend to have down-facing surfaces with different morphology compared to up-facing ones. In the context of cellular structures, this is likely to result in different thickness effects depending on direction measured.

For the FDM process, in addition to orientation, the toolpaths that effectively determine the internal meso-structure of the part (discussed in a previous blog post in greater detail) have a very strong influence on observed stiffness behavior, as shown in Figure 6. Thus orientation and process parameters are variables that need to be comprehended in the modeling of cellular structures – or set as constants for the range of applicability of the model parameters that are derived from a certain set of process conditions.

## Summary

Modeling cellular structures has the above mentioned challenges – most have practical implications in determining what is the correct specimen design – it is our mission over the next 18 months to address some of these challenges to a satisfactory level through an America Makes grant we have been awarded. While these ideas have been explored in other manufacturing contexts,  much remains to be done for the AM community, where cellular structures have a singular potential in application.

In future posts, I will discuss some of these challenges in detail and also discuss different approaches to modeling 3D printed cellular structures – they do not always address all the challenges here satisfactorily but each has its pros and cons. Until then, feel free to send us an email at info@padtinc.com citing this blog post, or connect with me on LinkedIn so you get notified whenever I write a post on this, or similar subjects in Additive Manufacturing (1-2 times/month).

## Classification of Cellular Solids (and why it matters)

Updated (8/30/2016): Two corrections made following suggestions by Gilbert Peters: the first corrects the use of honeycomb structures in radiator grille applications as being for flow conditioning, the second corrects the use of the Maxwell stability criterion, replacing the space frame example with an octet truss.

~

This is my first detailed post in a series on cellular structures for additive manufacturing, following an introductory post I wrote where I classified the research landscape in this area into four elements: design, analysis, manufacturing and implementation.

Within the design element, the first step in implementing cellular structures in Additive Manufacturing (AM) is selecting the appropriate unit cell(s). The unit cell is selected based on the performance desired of it as well as the manufacturability of the cells. In this post, I wish to delve deeper into the different types of cellular structures and why the classification is important. This will set the stage for defining criteria for why certain unit cell designs are preferable over others, which I will attempt in future posts. This post will also explain in greater detail what a “lattice” structure, a term that is often erroneously used to describe all cellular solids, truly is.

# 1. Honeycomb

1.1 Definition
Honeycombs are prismatic, 2-dimensional cellular designs extruded in the 3rd dimension, like the well-known hexagonal honeycomb shown in Figure 1. All cross-sections through the 3rd dimension are thus identical, making honeycombs somewhat easy to model. Though the hexagonal honeycomb is most well known, the term applies to all designs that have this prismatic property, including square and triangular honeycombs. Honeycombs have a strong anisotropy in the 3rd dimension – in fact, the modulus of regular hexagonal and triangular honeycombs is transversely isotropic – equal in all directions in the plane but very different out-of-plane.

1.2 Design Implications
The 2D nature of honeycomb structures means that their use is beneficial when the environmental conditions are predictable and the honeycomb design can be oriented in such a way to extract maximum benefit. One such example is the crash structure in Figure 2 as well as a range of sandwich panels. Several automotive radiator grilles are also of a honeycomb design to condition the flow of air. In both cases, the direction of the environmental stimulus is known – in the former, the impact load, in the latter, airflow.

# 2. Open-Cell Foam

2.1 Definition
Freeing up the prismatic requirement on the honeycomb brings us to a fully 3-dimensional open-cell foam design as shown in one representation of a unit cell in Figure 3. Typically, open-cell foams are bending-dominated, distinguishing them from stretch-dominated lattices, which are discussed in more detail in a following section on lattices.

2.2 Design Implications
Unlike the honeycomb, open cell foam designs are more useful when the environmental stimulus (stress, flow, heat) is not as predictable and unidirectional. The bending dominated mechanism of deformation make open-cell foams ideal for energy absorption – stretch dominated structures tend to be stiffer. As a result of this, applications that require energy absorption such as mattresses and crumple zones in complex structures. The interconnectivity of open-cell foams also makes them a candidate for applications requiring fluid flow through the structure.

# 3. Closed-Cell Foam

3.1 Definition
As the name suggests, closed cell foams are open-cell foams with enclosed cells, such as the representation shown in Figure 6. This typically involves a membrane like structure that may be of varying thickness from the strut-like structures, though this is not necessary. Closed-cell foams arise from a lot of natural processes and are commonly found in nature. In man-made entities, they are commonly found in the food industry (bread, chocolate) and in engineering applications where the enclosed cell is filled with some fluid (like air in bubble wrap, foam for bicycle helmets and fragile packaging).

3.2 Design Implications
The primary benefit of closed cell foams is the ability to encapsulate a fluid of different properties for compressive resilience. From a structural standpoint, while the membrane is a load-bearing part of the structure under certain loads, the additional material and manufacturing burden can be hard to justify. Within the AM context, this is a key area of interest for those exploring 3D printing food products in particular but may also have value for biomimetic applications.

# 4. Lattice

4.1 Definition
Lattices are in appearance very similar to open cell foams but differ in that lattice member deformation is stretch-dominated, as opposed to bending*. This is important since for the same material allocation, structures tend to be stiffer in tension and/or compression compared to bending – by contrast, bending dominated structures typically absorb more energy and are more compliant.

So the question is – when does an open cell foam become stretch dominated and therefore, a lattice? Fortunately, there is an app equation for that.

Maxwell’s Stability Criterion
Maxwell’s stability criterion involves the computation of a metric M for a lattice-like structure with b struts and j joints as follows:

In 2D structures: M = b – 2j + 3
In 3D structures:
M = b – 3j + 6

Per Maxwell’s criterion, for our purposes here where the joints are locked (and not pinned), if M < 0, we get a structure that is bending dominated. If M >= 0, the structure is stretch dominated. The former constitutes an open-cell foam, the latter a lattice.

There are several approaches to establishing the appropriateness of a lattice design for a structural applications (connectivity, static and kinematic determinism etc.) and how they are applied to periodic structures and space frames. It is easy for one (including for this author) to confuse these ideas and their applicability. For the purposes of AM, Maxwell’s Stability Criterion for 3D structures is a sufficient condition for static determinancy. Further, for a periodic structure to be truly space-filling (as we need for AM applications), there is no simple rigid polyhedron that fits the bill – we need a combination of polyhedra (such as an octahedron and tetrahedron in the octet truss shown in the video below) to generate true space filling, and rigid structures. The 2001 papers by Deshpande, Ashby and Fleck illustrate these ideas in greater detail and are referenced at the end of this post.

Video: The octet truss is a classic stretch-dominated structure, with b = 36 struts, j = 14 joints and M = 0 [Attr. Lawrence Livermore National Labs]

4.2 Design Implications
Lattices are the most common cellular solid studied in AM – this is primarily on account of their strong structural performance in applications where high stiffness-to-weight ratio is desired (such as aerospace), or where stiffness modulation is important (such as in medical implants). However, it is important to realize that there are other cellular representations that have a range of other benefits that lattice designs cannot provide.

# Conclusion: Why this matters

It is a fair question to ask why this matters – is this all just semantics? I would like to argue that the above classification is vital since it represents the first stage of selecting a unit cell for a particular function. Generally speaking, the following guidelines apply:

• Open cell foams where energy absorption and compliance is important
• Closed cell foams for fluid-filled and hydrostatic applications
• Lattice structures where stiffness and resistance to bending is critical

Finally, another reason it is important to retain the bigger picture on all cellular solids is it ensures that the discussion of what we can do with AM and cellular solids includes all the possibilities and is not limited to only stiffness driven lattice designs.

Note: This blog post is part of a series on “Additive Manufacturing of Cellular Solids” that I am writing over the coming year, diving deep into the fundamentals of this exciting and fast evolving topic. To ensure you get each post (~2 a month) or to give me feedback for improvement, please connect with me on LinkedIn.

## References

[1] Ashby, “Materials Selection in Mechanical Design,” Fourth Edition, 2011
[2] Gibson & Ashby, “Cellular Solids: Structure & Properties,” Second Edition, 1997
[3] Gibson, Ashby & Harley, “Cellular Materials in Nature & Medicine,” First Edition, 2010
[4] Ashby, Evans, Fleck, Gibson, Hutchinson, Wadley, “Metal Foams: A Design Guide,” First Edition, 2000
[5] Deshpande, Ashby, Fleck, “Foam Topology Bending versus Stretching Dominated Architectures,” Acta Materialia 49, 2001
[6] Deshpande, Fleck, Ashby, “Effective properties of the octet-truss lattice material,”  Journal of the Mechanics and Physics of Solids, 49, 2001

## Notes

* We defer to reference [1] in distinguishing lattice structures as separate from foams – this is NOT the approach used in [2] and [3] where lattices are treated implicitly as a subset of open-cell foams. The distinction is useful from a structural perspective and as such is retained here.

## The Additive Manufacturing Cellular Solids Research Landscape

I am writing this post after visiting the 27th SFF Symposium, a 3-day Additive Manufacturing (AM) conference held annually at the University of Texas at Austin. The SFF Symposium stands apart from other 3D printing conferences held in the US (such as AMUG, RAPID and Inside3D) in the fact that about 90% of the attendees and presenters are from academia. This year had 339 talks in 8 concurrent tracks and 54 posters, with an estimated 470 attendees from 20 countries – an overall 50% increase over the past year.

As one would expect from a predominantly academic conference, the talks were deeper in their content and tracks were more specialized. The track I presented in (Lattice Structures) had a total of 15 talks – 300 minutes of lattice talk, which pretty much made the conference for me!

In this post, I wish to summarize the research landscape in AM cellular solids at a high level: this classification dawned on me as I was listening to the talks over two days and taking in all the different work going on across several universities. My attempt in this post is to wrap my arms around the big picture and show how all these elements are needed to make cellular solids a routine design feature in production AM parts.

## Classification of Cellular Solids

First, I feel the need to clarify a technicality that bothered me a wee bit at the conference: I prefer the term “cellular solids” to “lattices” since it is more inclusive of honeycomb and all foam-like structures, following Gibson and Ashby’s 1997 seminal text of the same name. Lattices are generally associated with “open-cell foam” type structures only – but there is a lot of room for honeycomb structures and close-cell foams, each having different advantages and behaviors, which get excluded when we use the term “lattice”.

## The AM Cellular Solids Research Landscape

The 15 papers at the symposium, and indeed all my prior literature reviews and conference visits, suggested to me that all of the work in this space falls into one or more of four categories shown in Figure 2. For each of the four categories (design, analysis, manufacturing & implementation), I have listed below the current list of capabilities (not comprehensive), many of which were discussed in the talks at SFF. Further down I list the current challenges from my point of view, based on what I have learned studying this area over the past year.

Over the coming weeks I plan to publish a post with more detail on each of the four areas above, summarizing the commercial and academic research that is ongoing (to the best of my knowledge) in each area. For now, I provide below a brief elaboration of each area and highlight some important research questions.

### 1. Representation (Design)

This deals with how we incorporate cellular structures into our designs for all downstream activities. This involves two aspects: the selection of the specific cellular design (honeycomb or octet truss, for example) and its implementation in the CAD framework. For the former, a key question is: what is the optimum unit cell to select relative to performance requirements, manufacturability and other constraints? The second set of challenges arises from the CAD implementation: how does one allow for rapid iteration with minimal computational expense, how do cellular structures cover the space and merge with the external skin geometry seamlessly?

### 2. Optimization (Analysis)

Having tools to incorporate cellular designs is not enough – the next question is how to arrange these structures for optimum performance relative to specified requirements? The two most significant challenges in this area are performing the analysis at reasonable computational expense and the development of material models that accurately represent behavior at the cellular structure level, which may be significantly different from the bulk.

### 3. Realization (Manufacturing)

Manufacturing cellular structures is non-trivial, primarily due to the small size of the connecting members (struts, walls). The dimensions required are often in the order of a few hundred microns and lower, which tends to push the capabilities of the AM equipment under consideration. Additionally, in most cases, the cellular structure needs to be self-supporting and specifically for powder bed fusion, must allow for removal of trapped powder after completion of the build. One way to address this is to develop a map that identifies acceptable sizes of both the connecting members and the pores they enclose. For this, we need robust ways of monitoring quality of AM cellular solids by using in-situ and Non-Destructive techniques to guard against voids and other defects.

### 4. Application (Implementation)

Cellular solids have a range of potential applications. The well established ones include increasing stiffness-to-weight ratios, energy absorption and thermal performance. More recent applications include improving bone integration for implants and modulating stiffness to match biological distributions of material (biomimicry), as well as a host of ideas involving meta-materials. The key questions here include how do we ensure long term reliability of cellular structures in their use condition? How do we accurately identify and validate these conditions? How do we monitor quality in the field? And how do we ensure the entire life cycle of the product is cost-effective?

## So What?

I wrote this post for two reasons: I love to classify information and couldn’t help myself after 5 hours of hearing and thinking about this area. But secondly, I hope it helps give all of us working in this space context to engage and communicate more seamlessly and see how our own work fits in the bigger picture.

A lot of us have a singular passion for the overlapping zone of AM and cellular solids and I can imagine in a few years we may well have a conference, an online journal or a forum of some sort just dedicated to this field – in fact, I’d love to assess interest in such an effort or an equivalent collaborative exercise. If this idea resonates with you, please connect with me on LinkedIn and drop me a note, or send us an email (info@padtinc.com) and cite this blog post so it finds its way to me.

## Additive Manufacturing – Back to the Future!

Most histories of Additive Manufacturing (3D printing) trace the origins of the technology back to Charles Hull’s 1984 patent, the same year production began on the first of the Back to the Future movies. Which is something of a shock when you see 3D printing dotting the Gartner Hype Cycle like it was invented in the post-Seinfeld era. But that is not what this post is about.

When I started working on Additive Manufacturing (AM), I was amazed at the number of times I was returning to text books and class notes I had used in graduate school a decade ago. This led me to reflect on how AM is helping bring back to the forefront disciplines that had somehow lost their cool factor – either by becoming part of the old normal, or because they contained ideas that were ahead of their time. I present three such areas of research that I state, with only some exaggeration, were waiting for AM to come along.

• Topology Optimization: I remember many a design class where we would discuss topology optimization, look at fancy designs and end with a conversation that involved one of the more cynical students asking “All that’s fine, but how are you going to make that?”. Cue the elegant idea of building up a structure layer-by layer. AM is making it possible to manufacture parts with geometries that look like they came right out of a stress contour plot. And firms such as ANSYS, Autodesk and Altair, as well as universities and labs are all working to improve their capabilities at the intersection of topology optimization and additive manufacturing.

• Lattice Structures: One of the first books I came across when I joined PADT was a copy of Cellular Solids by Lorna Gibson and M.F. Ashby. Prof. Gibson’s examples of these structures as they occur in nature demonstrate how they provide an economy of material usage for the task at hand. Traditionally, in engineering structures, cellular designs are limited to foams or consistent shapes like sandwich panels where the variation in cell geometry is limited – this is because manufacturing techniques do not normally lend themselves well to building complex, three dimensional structures like those found in nature. With AM technologies however, cell sizes and structures can be varied and densities modified depending on the design of the structure and the imposed loading conditions, making this an exciting area of research.

• Metallurgy: As I read the preface to my “Metallurgy for the Non-Metallurgist” text book, I was surprised to note the author openly bemoan the decline of interest in metallurgy, and subsequently, fewer metallurgists in the field. And I guess it makes sense: materials science is today mostly concerned with much smaller scales than the classical metallurgist trained in. Well, lovers of columnar grain growth and precipitation hardening can now rejoice – metallurgy is at the very heart of AM technology today – most of the projected growth in AM is in metals. The science of powder metallurgy and the microstructure-property-process relationships of the metal AM technologies are vital building blocks to our understanding of metal 3D printing. Luckily for me, I happen to possess a book on powder metallurgy. And it too, is from 1984.