Hyperelas2

A hyper elastic cube in compression.

Static Scene

Interactive Scene

Static Scene

Interactive Scene

===== hyperelastic problem at iteration 0 =====
At Newton iteration 1 norm is 3.873448821042e-01
At Newton iteration 2 norm is 9.538664605952e-03
At Newton iteration 3 norm is 8.884193859163e-04
At Newton iteration 4 norm is 5.678328312556e-05
At Newton iteration 5 norm is 4.665333617495e-07

import numpy as np

from EasyFEA import Display, ElemType, Models, Simulations, PyVista
from EasyFEA.Geoms import Domain

if __name__ == "__main__":
    Display.Clear()

    # ----------------------------------------------
    # Configuration
    # ----------------------------------------------
    dim = 3

    # geom
    L = 1
    h = 1

    # ----------------------------------------------
    # Mesh
    # ----------------------------------------------
    meshSize = h / 10

    contour = Domain((0, 0), (L, h), h / 10)

    if dim == 2:
        mesh = contour.Mesh_2D([], ElemType.QUAD8, isOrganised=True)
    else:
        mesh = contour.Mesh_Extrude(
            [], [0, 0, h], [h / meshSize], ElemType.HEXA8, isOrganised=True
        )

    nodesX0 = mesh.Nodes_Conditions(lambda x, y, z: x == 0)
    nodesXL = mesh.Nodes_Conditions(lambda x, y, z: x == L)

    # ----------------------------------------------
    # Simulation
    # ----------------------------------------------

    isot = Models.Elastic.Isotropic(dim, E=1, v=0.3)
    lmbda = isot.get_lambda()
    mu = isot.get_mu()
    mat = Models.HyperElastic.SaintVenantKirchhoff(dim, lmbda, mu, thickness=h)

    simu = Simulations.HyperElastic(mesh, mat)

    uc = -0.3
    simu.add_dirichlet(nodesX0, [0] * dim, simu.Get_unknowns())
    values = np.zeros_like(simu.Get_unknowns(), dtype=float)
    values[0] = uc
    simu.add_dirichlet(nodesXL, values, simu.Get_unknowns())

    simu.Solve()

    # ----------------------------------------------
    # Results
    # ----------------------------------------------

    PyVista.Plot_BoundaryConditions(simu).show()
    PyVista.Plot(simu, "ux", 1, plotMesh=True).show()