sedaoturak / linear_elasticity_3d_fenics

Finite element modeling for linear elasticity problem in 3D by using FEniCS software

Linear Elasticity Problem in 3D by FEniCS

This repository shows how to model a three dimensional (3D) part and simulate for a linear elasticity problem. FEniCS, Python based open-source partial differential equation software, is used to simulate the part, while 3D cylindrical part is modeled by pygmsh library -python package for gmsh software.

The code is written in Google Colab. The dependencies or installation details may change dor local machines or other platforms.

This example is expected to be helpful to students, researchers or anyone who starts to use FEniCS.

Any contribution is more than welcome!

Geometry

The part is a cylinder with radius of 20 mm and height of 20 mm. It was created by using pygmsh library which is the Python interface of Gmsh mesh generator software.

Material

Material is chosen as aluminum to show the linear elasticity behaviour.

Problem Definition

The part is simulated under applied displacement in compression mode. There is also a mode for applying force application in the problem definition.

As boundary conditions, the bottom surface of the part is fixed in x, y and z directions, while on top, the surface is free in z direction and fixed in x and y directions.

Visualizations

Plotly library is used to vusalize the part, mesh and the deformed body after the simulation.