3D topology optimization of continuous fiber-reinforced structures via natural evolution method

The method to optimize a topology of 3D continuous fiber-reinforced additively manufactured structures is discussed. The proposed method makes it possible to simultaneously search for density distribution and local reinforcement layup in 3D composite structures of transversely isotropic materials. The approach uses a dynamical systems method to find density distribution, combined with the method for rotation of reinforcement direction to align it in the direction of principal stresses with local minimum compliance. The algorithm is implemented as a built-in material model within Abaqus finite element suite. Both the optimal material density distribution and the distribution of fiber orientation vector are determined for three structural elements used as benchmarks: the bending of simply supported 2D beam under central point load, the loading of 3D cube by vertical load, and the bending of 3D cantilever beam.