Skip to main content

Finite Element Analysis

Project description

Finite element analysis for continuum mechanics of solid bodies.

FElupe PyPI version shields.io PyPI - Python Version Documentation Status License: GPL v3 Project Status: Active – The project has reached a stable, usable state and is being actively developed. codecov DOI Codestyle black PyPI - Downloads lite-badge Open In Colab

FElupe is a Python 3.9+ 🐍 finite element analysis package 📦 focusing on the formulation and numerical solution of nonlinear problems in continuum mechanics of solid bodies 🔧. This package is intended for scientific research 💻, but is also suitable for running nonlinear simulations 🚂 in general 🏎️. In addition to the transformation of general weak forms into sparse vectors and matrices, FElupe provides an efficient high-level abstraction layer for the simulation of the deformation of solid bodies.

Highlights

Efficient NumPy-based math is realized by element-wise operating trailing axes [1]. The finite element method, as used in FElupe, is based on [2], [3] and [4]. Related scientific articles are listed in the sections of the API reference.

[!NOTE] FElupe is a combination of FE (finite element) and the german word Lupe 🔍 (magnifying glass) as a synonym for getting an insight 📖 how a finite element analysis code 🧮 looks like under the hood 🕳️.

Installation

Install Python, fire up 🔥 a terminal and run 🏃

pip install felupe[all]

The documentation covers more details, like required and optional dependencies and how to install the latest development version.

Getting Started

Expand the description...

This tutorial covers the essential high-level parts of creating and solving problems with FElupe. As an introductory example 👨‍🏫, a quarter model of a solid cube with hyperelastic material behaviour is subjected to a uniaxial elongation applied at a clamped end-face.

First, let’s import FElupe and create a meshed cube out of hexahedron cells with a given number of points per axis. A numeric region, pre-defined for hexahedrons, is created on the mesh. A vector-valued displacement field is initiated on the region. Next, a field container is created on top of this field.

A uniaxial load case is applied on the displacement field stored inside the field container. This involves setting up symmetry planes as well as the absolute value of the prescribed displacement at the mesh-points on the right-end face of the cube. The right-end face is clamped 🛠️: only displacements in direction x are allowed. The dict of boundary conditions for this pre-defined load case are returned as boundaries and the partitioned degrees of freedom as well as the external displacements are stored within the returned dict loadcase.

An isotropic pseudo-elastic Ogden-Roxburgh Mullins-softening model formulation in combination with an isotropic hyperelastic Neo-Hookean material formulation is applied on the displacement field of a nearly-incompressible solid body.

A step generates the consecutive substep-movements of a given boundary condition. The step is further added to a list of steps of a job 👩‍💻 (here, a characteristic-curve 📈 job is used). During evaluation ⏳, each substep of each step is solved by an iterative Newton-Rhapson procedure ⚖️. The solution is exported after each completed substep as a time-series ⌚ XDMF file. Finally, the result of the last completed substep is plotted.

For more details beside this high-level code snippet, please have a look at the 📝 documentation.

import felupe as fem

mesh = fem.Cube(n=6)
region = fem.RegionHexahedron(mesh)
field = fem.FieldContainer([fem.Field(region, dim=3)])

boundaries, loadcase = fem.dof.uniaxial(field, clamped=True)

umat = fem.OgdenRoxburgh(material=fem.NeoHooke(mu=1), r=3, m=1, beta=0.1)
solid = fem.SolidBodyNearlyIncompressible(umat, field, bulk=5000)

move = fem.math.linsteps([0, 1, 0, 1, 2, 1], num=5)
step = fem.Step(items=[solid], ramp={boundaries["move"]: move}, boundaries=boundaries)

job = fem.CharacteristicCurve(steps=[step], boundary=boundaries["move"])
job.evaluate(filename="result.xdmf")
fig, ax = job.plot(
    xlabel="Displacement $u$ in mm $\longrightarrow$",
    ylabel="Normal Force $F$ in N $\longrightarrow$",
)

solid.plot("Principal Values of Cauchy Stress").show()

Curve Solid Body

Documentation

The documentation is located here.

Extension Packages

The capabilities of FElupe may be enhanced with extension packages created by the community.

Package Description
hyperelastic Constitutive hyperelastic material formulations
matadi Material Definition with Automatic Differentiation (AD)
tensortrax Differentiable Tensors based on NumPy Arrays
feplot A visualization tool for FElupe

Testing

To run the FElupe unit tests, check out this repository and type

tox

Scientific Publications

A list of articles in which FElupe is involved. If you use FElupe in your scientific work, please star this repository, cite it DOI and add your publication to this list.

Expand the list...

Changelog

All notable changes to this project will be documented in this file. The format is based on Keep a Changelog, and this project adheres to Semantic Versioning.

References

[1] T. Gustafsson and G. McBain, "scikit-fem: A Python package for finite element assembly", Journal of Open Source Software, vol. 5, no. 52. The Open Journal, p. 2369, Aug. 21, 2020. DOI:10.21105/joss.02369.

[2] J. Bonet and R. D. Wood, "Nonlinear Continuum Mechanics for Finite Element Analysis". Cambridge University Press, Mar. 13, 2008. DOI:10.1017/cbo9780511755446.

[3] K. J. Bathe, "Finite Element Procedures". 2006, isbn: 978-0-9790049-0-2.

[4] O. C. Zienkiewicz, R. L. Taylor and J. Z. Zhu, "The Finite Element Method: its Basis and Fundamentals". Elsevier, 2013. DOI:10.1016/c2009-0-24909-9.

License

FElupe - finite element analysis (C) 2021-2024 Andreas Dutzler, Graz (Austria).

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see https://www.gnu.org/licenses/.

Project details


Release history Release notifications | RSS feed

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

felupe-9.2.0.tar.gz (270.2 kB view details)

Uploaded Source

Built Distribution

felupe-9.2.0-py3-none-any.whl (374.6 kB view details)

Uploaded Python 3

File details

Details for the file felupe-9.2.0.tar.gz.

File metadata

  • Download URL: felupe-9.2.0.tar.gz
  • Upload date:
  • Size: 270.2 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.12.9

File hashes

Hashes for felupe-9.2.0.tar.gz
Algorithm Hash digest
SHA256 581fb4ed76a00b3c129ecc5a7c69104a5497d7858a0c00117c418f8dea73c93c
MD5 b8653a90c2e6ebd10e39b2b528e11860
BLAKE2b-256 47728c280b0d753c8f7d78301a0db5fe3ed389f220bf34d0abde635c68abb96d

See more details on using hashes here.

File details

Details for the file felupe-9.2.0-py3-none-any.whl.

File metadata

  • Download URL: felupe-9.2.0-py3-none-any.whl
  • Upload date:
  • Size: 374.6 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.12.9

File hashes

Hashes for felupe-9.2.0-py3-none-any.whl
Algorithm Hash digest
SHA256 6e2c06bec8199202d65b3678b63ffa0b6c217e05fff12667c12bdf1fe9bef445
MD5 bd08020821e7b1cd19fcf14fdc70982d
BLAKE2b-256 b21f82b8acf344bbde23f7343a80fa30710c330eea9eb3cd3a790f983a8fb9f3

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page