felupe 1.5.0

Finite Element Analysis

FElupe - Finite Element Analysis

FElupe is a Python 3.6+ finite element analysis package focussing on the formulation and numerical solution of nonlinear problems in continuum mechanics of solid bodies. Its name is a combination of FE (finite element) and the german word Lupe (magnifying glass) as a synonym for getting a little 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]

where [all] installs all optional dependencies. By default, FElupe does not require numba and sparse. In order to make use of all features of FElupe, it is suggested to install all optional dependencies.

Hello, FElupe!

A quarter model of a solid cube with hyperelastic material behavior is subjected to a uniaxial elongation applied at a clamped end-face. This involves the creation of a mesh, a region and a displacement field. Furthermore, the boundary conditions are created by a template for a uniaxial loadcase. The material behavior is defined through a FElupe-built-in Neo-Hookean material formulation. Inside a Newton-Rhapson procedure, the internal force vector and the tangent stiffness matrix are generated by assembling both linear and bilinear forms of static equilibrium. Finally, the solution of the incremental displacements is calculated und updated until convergence is reached. For more details beside this high-level code snippet, please have a look at the documentation.

import felupe as fe

# create a hexahedron-region on a cube
region = fe.RegionHexahedron(fe.Cube(n=11))

# add a displacement field and apply a uniaxial elongation on the cube
displacement = fe.Field(region, dim=3)
boundaries, dof0, dof1, ext0 = fe.dof.uniaxial(displacement, move=0.2, clamped=True)

# define the constitutive material behavior
umat = fe.NeoHooke(mu=1.0, bulk=2.0)

# newton-rhapson procedure
res = fe.newtonrhapson(displacement, umat=umat, dof1=dof1, dof0=dof0, ext0=ext0)

# save result
fe.save(region, res.x, filename="result.vtk")

Documentation

The documentation is located here.

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.

[Unreleased]

[1.5.0] - 2021-11-29

Added

• Add kwargs of field.extract() to fun and jac of newtonrhapson.

Changed

• Set default number of threads in MatadiMaterial to multiprocessing.cpu_count().
• Moved documentation to Read the Docs (Sphinx).

Fixed

• Fix dim in calculation of reaction forces (tools.force) for FieldMixed.
• Fix calculation of reaction moments (tools.moment) for FieldMixed.

[1.4.0] - 2021-11-15

Added

• Add mask argument to Boundary for the selection of user-defined points.
• Add shear loadcase.
• Add a wrapper for matadi materials as MatadiMaterial.
• Add verbose and timing arguments to newtonrhapson.

Fixed

• Obtain internal dim from Field in calculation of reaction force tools.force.
• Fix math.dot for combinations of rank 1 (vectors), rank 2 (matrices) and rank 4 tensors.

[1.3.0] - 2021-11-02

Changed

• Rename mesh.as_discontinous() to mesh.disconnect().
• Rename constitution.Mixed to constitution.ThreeFieldVariation.
• Rename unstack to offsets as return of dof-partition and all subsequent references.
• Import tools (newtonrhapson, project, save) and constitution (NeoHooke, LinearElastic and ThreeFieldVariation) to FElupe's namespace.
• Change minimal README-example to a high-level code snippet and refer to docs for details.

[1.2.0] - 2021-10-31

Added

• Add template regions, i.e. a region with a Hexahedron() element and a quadrature scheme GaussLegendre(order=1, dim=3) as RegionHexahedron, etc.
• Add biaxial and planar loadcases (like uniaxial).
• Add a minimal README-example (Hello FElupe!).

Changed

• Deactivate clamped boundary (clamped=False) as default option for uniaxial loading dof.uniaxial.

[1.1.0] - 2021-10-30

Added

• Add inverse quadrature method quadrature.inv() for Gauss-Legendre schemes.
• Add discontinous representation of a mesh as mesh method mesh.as_discontinous().
• Add tools.project() to project (and average) values at quadrature points to mesh points.

Changed

• Removed quadpy dependency and use built-in polynomials of numpy for Gauss-Legendre calculation.

Fixed

• Fix typo in first shear component of math.tovoigt() function.
• Fix wrong stress projection in tools.topoints() due to different quadrature and cell ordering.

[1.0.1] - 2021-10-19

Fixed

• Fix import of dof-module if sparse is not installed.

[1.0.0] - 2021-10-19

Added

• Start using a Changelog.
• Added docstrings for essential classes, methods and functions.
• Add array with point locations for all elements.

Changed

• Rename element methods (from basis to function and from basisprime to gradient).
• Make constitutive materials more flexible (allow material parameters to be passed at stress and elasticity evaluation umat.gradient(F, mu=1.0)).
• Rename ndim to dim.
• Simplify element base classes.
• Speed-up calculation of indices (rows, cols) for Fields and Forms (about 10x faster now).
• Update test_element.py according to changes in element methods.

Removed

• Automatic check if the gradient of a region can be calculated based on the dimensions. The grad argument in region(grad=False) has to be enforced by the user.

License

FElupe - finite element analysis (C) 2021 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/.