Skip to content

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

Notifications You must be signed in to change notification settings

sedaoturak/linear_elasticity_3D_fenics

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

12 Commits
 
 
 
 

Repository files navigation

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.

def

Visualizations

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