Topology optimization based on the Moving Node Algorithm using EFG, FEM or IIEFG methods.
This set of Matlab files are used for topology optimization using the Moving Node Approach (MNA) for 2D plane strain problems. In this approach, the material distribution is decoupled from the discretization.
The material distribution is used to specify where the material is. It is based on mass nodes. The density at one point can be computed thanks to a kernel approximation using cubic spline shape functions.
Some corrections are added in order to avoid numerical problems.
The governing linear elasticity equations have to be discretized to solve numerically the problem. The discretization methods can be :
- A meshless method called Element-Free Galerkin (EFG)
- A Finite Element Method (FEM)
The main file topologyOptimization.m is a Matlab script that launches the optimizers. The user can change the script to set
- The discretization method (EFG or FEM)
- The optimization variables (mass nodes, undeformable structural members or deformable structural members)
- The optimization algorithm
The problem constants are defined in the Constants\ directory. This includes
- Problem geometry
- Boundary conditions
- Material distribution constants
The discretization methods use functions from the Discretization\, EFG\ and FEM\ directories.
The files InitEFGMesh.m and InitFEMMesh.m contain the discretization constants. This includes
- Number of nodes or element per dimension
- Shape functions degree
- Number of Gauss points per cell/element
- Relative smoothing length (EFG)
The other files of these directories should not be modified directly by the user.
The optimizers as well as some useful functions are placed in the Optimization\ directory. Depending on the algorithm, some specific constants can be found in this directory files.