import pyvista
from dolfinx import mesh, fem, plot, io, default_scalar_type
from dolfinx.fem.petsc import LinearProblem
from mpi4py import MPI
from ufl import FiniteElement, MixedElement # Ensure this import is correct
import numpy as np
L = 72
W = 4.0
mu = 1
rho = 1
delta = W / L
gamma = 0.4 * delta**2
beta = 1.25
lambda_ = beta
g = gamma
domain = mesh.create_box(MPI.COMM_WORLD, [np.array([0, 0, 0]), np.array([L, W, W])],
[20, 6, 6], cell_type=mesh.CellType.hexahedron)
# Create two elements: one for displacements (3 DOFs) and one for rotations (3 DOFs)
displacement_element = FiniteElement("Lagrange", domain.ufl_cell(), 1, dim=3)
rotation_element = FiniteElement("Lagrange", domain.ufl_cell(), 1, dim=3)
# Combine the elements into a mixed element
element = MixedElement([displacement_element, rotation_element])
# Create a function space with the specified element
V = fem.FunctionSpace(domain, element)
Traceback (most recent call last):
File "/home/prusso/dolfinx-beam_w/dolfinx_beam_w.py", line 2, in <module>
from dolfinx import mesh, fem, plot, io, default_scalar_type
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/dolfinx/__init__.py", line 25, in <module>
from dolfinx import fem, geometry, graph, io, jit, la, log, mesh, nls, plot, utils
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/dolfinx/fem/__init__.py", line 23, in <module>
from dolfinx.fem.assemble import (
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/dolfinx/fem/assemble.py", line 22, in <module>
from dolfinx.fem.forms import Form
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/dolfinx/fem/forms.py", line 15, in <module>
import ufl
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/__init__.py", line 349, in <module>
from ufl.formoperators import replace, derivative, action, energy_norm, rhs, lhs,\
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/formoperators.py", line 46, in <module>
from ufl.algorithms import compute_form_adjoint, compute_form_action
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/algorithms/__init__.py", line 95, in <module>
from ufl.algorithms.compute_form_data import compute_form_data
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/algorithms/compute_form_data.py", line 31, in <module>
from ufl.algorithms.formtransformations import compute_form_arities
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/algorithms/formtransformations.py", line 39, in <module>
from ufl.algorithms.transformer import Transformer
File "/home/prusso/spack/var/spack/environments/fenicsx-env/.spack-env/view/lib/python3.11/site-packages/ufl/algorithms/transformer.py", line 26, in <module>
from inspect import getargspec
ImportError: cannot import name 'getargspec' from 'inspect' (/home/prusso/spack/opt/spack/linux-ubuntu22.04-skylake/gcc-11.4.0/python-3.11.9-3k6gyosv23rjrhjsrblavq5la3tdu2z5/lib/python3.11/inspect.py)
What I would like to do is set up the function space for a stuctural element of 6 d.o.f.'s… A tetrahedron, or hexahedron would be great… Any way to do this with dolfinx?