Computer Methods in Applied Mechanics and Engineering - 2025

Self-support structure topology optimization for multi-axis additive manufacturing incorporated with curved layer slicing

A differentiable framework that concurrently optimizes structural topology and curved fabrication layers, enabling self-supporting designs for multi-axis additive manufacturing without sacrificing unnecessary structural performance.

Shuzhi Xu1, Jikai Liu2,*, Dong He3, Kai Tang4, Kentaro Yaji1,*
1 Department of Mechanical Engineering, Graduate School of Engineering, Osaka University, Japan
2 School of Mechanical Engineering, Shandong University, China
3 Department of Mechanical and Aerospace Engineering, Hong Kong University of Science and Technology, Hong Kong
4 Smart Manufacturing Thrust, Hong Kong University of Science and Technology (Guangzhou), China
* Corresponding authors
Workflow from concurrent optimization to curved layers and manufacturing paths

The optimized structural density and process fields are converted into curved layers, slice geometry, toolpaths, G-code, and a manufacturable part.

Abstract

Multi-axis additive manufacturing expands the range of printable geometries, but the slicing strategy directly affects support requirements and nozzle accessibility. This work embeds a numerically tractable process-planning model into density-based topology optimization. Heat diffusion first supplies a fabrication sequence, and a Poisson problem transforms it into a geodesic-like process field with nearly uniform gradient magnitude. Curved layer normals and structural boundary gradients then define a differentiable self-support measure. By optimizing the structural and pseudo-density fields together, the method balances physical performance, continuous curved slicing, and support-free manufacturability.

Curved layer slicing

A scalar process field provides continuous layer sequence and spatial distance information.

Self-support control

Layer normals and part-boundary gradients identify regions that violate the overhang requirement.

Concurrent design

Structural density and process pseudo-density evolve together in one optimization framework.

Differentiable pipeline

Linear PDEs and adjoint sensitivities allow the use of mature gradient-based optimizers.

Method

The method begins with a heat-diffusion problem defined between the printing start and end boundaries. Normalized temperature gradients guide a Poisson equation whose solution is the process field. Iso-values of this field define successive curved layers. In parallel, filtered and projected structural densities describe the part boundary. Comparing the curved-layer normal with the local boundary gradient yields a non-support field that is aggregated as an explicit optimization constraint.

Curved layer generation from a scalar process field
A process field with nearly uniform gradient magnitude is sampled at constant intervals to create curved fabrication layers.
Non-support area analysis workflow
Structural and process design fields are filtered, projected, and coupled to evaluate unsupported boundary regions.
Concurrent topology and process optimization framework
The optimization loop couples structural analysis, heat diffusion, the Poisson process field, self-support evaluation, sensitivities, and MMA updates.

Numerical Results

Two-dimensional cantilever and bridge examples, together with three-dimensional bracket and heat-sink cases, demonstrate that the process field can evolve with structural topology while retaining smooth, continuous layer patterns. The comparison uses conventional topology optimization, planar-layer self-support optimization, and the proposed multi-axis formulation.

6.4% 2D performance improvement over planar-layer optimization
4.7% 3D performance improvement over planar-layer optimization
84.566 J Compliance of the proposed multi-axis method in the 2D comparison
330.866 J Compliance of the proposed multi-axis method in the 3D comparison
Optimized short cantilever with curved slicing fields
Short cantilever example: boundary conditions, optimized topology and layers, and the coupled structural and process fields.
Optimized bracket, curved slices, and displacement simulation
The three-dimensional bracket combines a self-supporting topology with curved slicing and displacement verification.
Optimized heat sink, curved slices, and thermal simulation
The heat-sink case demonstrates the same concurrent strategy for a thermal topology optimization problem.
Performance comparison of three topology optimization strategies
The proposed Multi-Axis-OPT preserves more geometric freedom and improves performance relative to planar-layer self-support design in both 2D and 3D.

Manufacturability and Printing

Directly slicing an unconstrained topology can create unavoidable nozzle collisions or unstable deposition around horizontal connections. The proposed concurrent formulation modifies both topology and layer geometry to reduce these conflicts. For the demonstrated heat-sink workflow, normalized process iso-surfaces are converted into 175 STL layers with an average thickness of 0.50 mm plus or minus 0.02 mm. Boundary paths are extracted, converted to G-code, verified in CAM software, and finally used for printing.

Manufacturability comparison showing nozzle collision and powder leakage
Manufacturability analysis shows how incompatible layer geometry can cause nozzle collision or an unstable melt pool, while the optimized design avoids these conflicts.
Conversion of optimized paths to G-code and final printing validation
Additive-axis vector planning, toolpath generation, G-code simulation, and the final manufacturing validation.

Citation

@article{xu2025selfsupport,
  title   = {Self-support structure topology optimization for multi-axis
             additive manufacturing incorporated with curved layer slicing},
  author  = {Xu, Shuzhi and Liu, Jikai and He, Dong and Tang, Kai and Yaji, Kentaro},
  journal = {Computer Methods in Applied Mechanics and Engineering},
  volume  = {438},
  pages   = {117841},
  year    = {2025},
  doi     = {10.1016/j.cma.2025.117841}
}

Acknowledgements

This work acknowledges support from the National Natural Science Foundation of China under Grant 52105462, the Aeronautical Science Foundation of China under Grant 2022Z008001001, and JSPS KAKENHI under Grant 23H03799.