A MATLAB Script for Solving 2D/3D Minimum Compliance Problems using Anisotropic Mesh Adaptation

Kristian Ejlebjerg Jensen*

*Corresponding author for this work

Research output: Contribution to journalConference articlepeer-review

320 Downloads (Pure)


We present a pure MATLAB implementation for solving 2D/3D compliance minimization problems using the density method. A filtered design variable with a minimum length is computed using a Helmholtz-type differential equation. The optimality criteria is used as optimizer and to avoid local minima we apply continuation of an exponent that controls the stiffness associated with intermediate design variables. We constrain the volume from above and use the implementation to show that optimizations with dynamic meshes can save significant amounts of computational time compared to fixed meshes without introducing mesh dependence for the mesh topology. This is despite the fact that the dynamic meshes cause oscillations of the objective function, particular for coarse meshes in 3D. The meshes are generated using anisotropic mesh adaptation based on local mesh modifications and we extent these modifications to preserve the information required for interpolating the design variables between meshes. We exploit symmetry boundaries in 3D, but not in 2D. Dirichlet boundary conditions are used to prevent non-zero filtered design variables on free boundaries. Mesh adaptation involves substantial book keeping, so the implementation totals some 5,000 lines of MAT-LAB code, but the functions associated with the forward analysis, geometry/mesh setup and optimization are concise and well documented, so the implementation can be used as a starting point for research on related topics.
Original languageEnglish
JournalProcedia Engineering
Pages (from-to)102-114
Publication statusPublished - 2017
Event26th International Meshing Roundtable - Crowne Plaza Fira Center, Barcelona, Spain
Duration: 18 Sep 201721 Sep 2017


Conference26th International Meshing Roundtable
LocationCrowne Plaza Fira Center


  • Adaption
  • Elasticity
  • Compliance
  • Optimization


Dive into the research topics of 'A MATLAB Script for Solving 2D/3D Minimum Compliance Problems using Anisotropic Mesh Adaptation'. Together they form a unique fingerprint.

Cite this