parafem

The goal of this project is to have a working fem code supporting membrane and truss elements. The theory is taken from the paper EXPLICIT AlGORITHMS FOR THE NONLINEAR DYNAMICS OF SHELLS

Tasks/Todos

[ ] A working c++ source
[ ] Comparison of elements
[ ] Test framework
[ ] JSON inputs
[ ] Python binding with pybind11
[ ] Make faster (fixed size arrays, ...)
[ ] line-forces can be computed with a simplified geometry.
- Using a implicit integration will make the computation faster.

Tests and use-cases

Pillow-test

Compare klatsch-stange

Deformation

Deformation of glider due to aerodynamic forces acting on membrane-elements

Computing line-forces

Computing the top line forces which influence the line-geometry and line-lengths:

Python bindings

import parafem
import numpy as np

mat = parafem.TrussMaterial(1000)
mat.rho = 1
mat.d_structural = 0.0
mat.d_velocity = 1

node1 = parafem.Node(-1, 0, 0)
node2 = parafem.Node(0, 0, 0)
node3 = parafem.Node(1, 0, 0)

node1.fixed = np.array([0, 0, 0])
node3.fixed = np.array([0, 0, 0])
node2.add_external_force(np.array([0, 1, 0]))

truss1 = parafem.Truss([node1, node2], mat)
truss2 = parafem.Truss([node2, node3], mat)

case = parafem.Case([truss1, truss2])

writer = parafem.vtkWriter("/tmp/parafem/truss1_py/output")

for i in range(10000):
    case.explicit_step(0.01)
    if (i % 10) == 0:
        writer.writeCase(case, 0.3)

License

GPLv3 (see LICENSE]

booya-at / parafem

readme