Surrogate optimization of a lattice foot orthotic
Article [Accepted Manuscript]
Is part of
Computers in biology and medicine ; vol. 155.Publisher(s)
ElsevierAffiliation
- Université de Montréal. Faculté de médecine. École de kinésiologie et des sciences de l'activité physique
- École polytechnique (Montréal, Québec). Laboratoire de mécanique multi-échelles (LM2)
- Université de Montréal. Laboratoire de simulation et modélisation du mouvement
- Hôpital Sainte-Justine. Centre de recherche
Abstract(s)
Background:
Additive manufacturing enables to print patient-specific Foot Orthotics (FOs). In FOs featuring lattice structures, the variation of the cell’s dimensions provides a locally variable stiffness to meet the therapeutic needs of each patient. In an optimization problem, however, using explicit Finite Element (FE) simulation of lattice FOs with converged 3D elements is computationally prohibitive. This paper presents a framework to efficiently optimize the cell’s dimensions of a honeycomb lattice FO for flat foot condition.
Methods:
We built a surrogate based on shell elements whose mechanical properties were computed by the numerical homogenization technique. The model was submitted to a static pressure distribution of a flat foot and it predicted the displacement field for a given set of geometrical parameters of the honeycomb FO. This FE simulation was considered as a black-box and a derivative-free optimization solver was employed. The cost function was defined based on the difference between the predicted displacement by the model against a therapeutic target displacement.
Results:
Using the homogenized model as a surrogate significantly accelerated the stiffness optimization of the lattice FO. The homogenized model could predict the displacement field 78 times faster than the explicit model. When 2000 evaluations were required in an optimization problem, the computational time was reduced from 34 days to 10 hours using the homogenized model rather than explicit model. Moreover, in the homogenized model, there was no need to re-create and re-mesh the insole’s geometry in each iteration of the optimization. It was only required to update the effective properties.
Conclusion:
The presented homogenized model can be used as a surrogate within an optimization framework to customize cell’s dimensions of honeycomb lattice FO in a computationally efficient manner.