[9.3.0] - 2025-06-29

Added

Add Mesh.add_points(points) to update the mesh with additional points.
Add Mesh.clear_points_without_cells() to clear the list of points without cells (useful for center-points of multi-point constraints).
Release FElupe on conda-forge, starting with v9.2.0.
Add ConstitutiveMaterial.is_stable() which returns a boolean mask of stability for isotropic material model formulations. Note that this will require an additional volumetric part of the strain energy density function for hyperelastic material model formulations without a volumetric part.
Add the linear-elastic material formulation constitution.LinearElastic1D() and a truss-body mechanics.TrussBody().
Add RegionTruss with a truss element. This is a line element with a GaussLobatto(order=0) quadrature, i.e. with two quadrature points located at the boundaries.

Changed

Change the required setuptools-version in the build-system table of pyproject.toml to match PEP639 (setuptools>=77.0.3).
Change the labels to well-known labels for the URLs in pyproject.toml.
Change the first return value of ViewMaterial.uniaxial(), ViewMaterial.planar(), ViewMaterial.biaxial(), ViewMaterialIncompressible.uniaxial(), ViewMaterialIncompressible.planar(), ViewMaterialIncompressible.biaxial() from the stretch to a list of all three principal stretches.

Fixed

Fix the declaration of the (spdx identifier) license and license-file in pyproject.toml.
Fix missing import of TriQuadraticHexahedron in the top-level namespace.
Fix the path to docs/_static/logo_without_text.svg in docs/conf.py.
Fix a typo in the docstring of MeshContainer.from_unstructured_grid().
Fix CompositeMaterial for input lists of length 1.

What's Changed

Add Mesh.add_points(points) and Mesh.clear_points_without_cells() by @adtzlr in https://github.com/adtzlr/felupe/pull/947
Update license in pyproject.toml by @adtzlr in https://github.com/adtzlr/felupe/pull/954
Make license-files more explicit in pyproject.toml by @adtzlr in https://github.com/adtzlr/felupe/pull/955
pyproject.toml ensure PEP639 support by @adtzlr in https://github.com/adtzlr/felupe/pull/956
Fix missing import of TriQuadraticHexahedron by @adtzlr in https://github.com/adtzlr/felupe/pull/958
Fix the path to the logo in docs/conf.py by @adtzlr in https://github.com/adtzlr/felupe/pull/959
Fix a typo in the docstring of MeshContainer.from_unstructured_grid() by @adtzlr in https://github.com/adtzlr/felupe/pull/960
Enhance some docstrings (only minor changes) by @adtzlr in https://github.com/adtzlr/felupe/pull/961
Docs: Use matplotlib-directive for plots by @adtzlr in https://github.com/adtzlr/felupe/pull/962
Fix typo from seperated to separated by @tkoyama010 in https://github.com/adtzlr/felupe/pull/963
Add ConstitutiveMaterial.is_stable() and fix CompositeMaterial by @adtzlr in https://github.com/adtzlr/felupe/pull/965
Add mechanics.Truss and constitution.LinearElastic1D by @adtzlr in https://github.com/adtzlr/felupe/pull/966
Add RegionTruss by @adtzlr in https://github.com/adtzlr/felupe/pull/967

Full Changelog: https://github.com/adtzlr/felupe/compare/v9.2.0...v9.3.0

- Python
Published by adtzlr 11 months ago

Added

Add SolidBody.assemble.mass(density=1.0) and SolidBodyNearlyIncompressible.assemble.mass(density=1.0) to assemble the mass matrix.
Add SolidBody.evaluate.stress(field) to evaluate the (first Piola-Kirchhoff) stress tensor (engineering stress in linear elasticity).
Add a free-vibration modal analysis Step/Job FreeVibration(items, boundaries) with methods to evaluate FreeVibration.evaluate() and to extract field, frequency = FreeVibration.extract(n) its n-th result.
Add Boundary.plot() to plot the points and prescribed directions of a boundary.
Add BoundaryDict as a subclassed dict with methods to plot(), screenshot() and imshow().
Add a new argument to apply a callable on the assembled vector/matrix of a solid body, SolidBody(..., apply=None). This may be used to sum the list of sub-blocks instead of stacking them together, SolidBody(..., block=False, apply=None). This is useful for mixed formulations where both the deformation gradient and the displacement values are required.
Add support for non-symmetric bilinear mixed forms in IntegralForm.
Add element.Element to the top-level package namespace.
Raise a TypeError if a mesh container is used as the mesh-argument in a region. The error message explains that the container is not supported. A mesh must be used instead.
Add the Vertex element formulation.
Add a vertex-region RegionVertex.
Add MeshContainer.as_vertex_mesh() to create a merged vertex mesh from the meshes of the mesh container.
Add Field.from_mesh_container(container) to create a top-level field on a vertex mesh.
Add GaussLobatto quadrature. The order-argument is equal to the 1d-sample points minus two. This ensures identical minimum order-arguments for Gauss-Legendre and Gauss-Lobatto schemes.

Changed

The first Piola-Kirchhoff stress tensor is evaluated if ViewSolid(stress_type=None).
Autodetect the stress-type in SolidBody.plot(name) from name.
Enhance the hello_world(axisymmetric=False, planestrain=False, curve=False, xdmf=False, container=False) function with new arguments to customize the generated template script.
Enhance Boundary with added support for multiaxial prescribed values.
Enhance math.linsteps(..., values=0) with default values except for the column axis if axis is not None.
Link all field-values to the values of the first field if no other field is given in FieldContainer.link().
Change the default arguments from block=True to block=None in SolidBody.assemble.vector(block=None) and SolidBody.assemble.matrix(block=None) and define block on creation, SolidBody(..., block=True) instead.
Integrate and assemble None-items in the fun-argument of IntegralForm. On integration, None is returned. None will be assembled to an emtpy sparse matrix.
Enhance mesh.interpolate_line(mesh, xi, axis=None, ...) which by default uses a curve-progress variable when axis=None. The curve progress evaluation points xi must be within 0 <= xi <= 1.

Fixed

Fix Boundary(..., mode="and") by ignoring any undefined axis.
Fix tools.hello_world(planestrain=True) with the correct region RegionQuad for the plane-strain template.
Fix the cells-array in RegionQuadraticQuadBoundary for midpoints of edges.
Fix FieldContainer.plot() for regions without the shape-function gradients available.

Removed

Remove the unused shape-argument in element.Element(shape). Adopt the arbitrary-lagrange element to use its own dim-argument. This simplifies the creation of custom finite element formulations.

What's Changed

Add SolidBody.evaluate.mass(density=1.0) by @adtzlr in https://github.com/adtzlr/felupe/pull/906
Add SolidBody.evaluate.stress() and autodetect the stress-type in SolidBody.plot(name) from name by @adtzlr in https://github.com/adtzlr/felupe/pull/908
Add FreeVibration(items, boundaries) by @adtzlr in https://github.com/adtzlr/felupe/pull/909
Enhance Boundary: Add support for multiaxial prescribed value-argument by @adtzlr in https://github.com/adtzlr/felupe/pull/912
Add Boundary.plot() and fix Boundary(..., mode="and") by @adtzlr in https://github.com/adtzlr/felupe/pull/913
Add BoundaryDict by @adtzlr in https://github.com/adtzlr/felupe/pull/915
Add PointLoad.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/916
Enhance math.linsteps(..., values=0.0) by @adtzlr in https://github.com/adtzlr/felupe/pull/917
Enhance FieldContainer.link() by @adtzlr in https://github.com/adtzlr/felupe/pull/918
Add SolidBody(..., block=True, apply=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/920
Enhance IntegralForm: Non-symmetric mixed forms and None by @adtzlr in https://github.com/adtzlr/felupe/pull/922
Simplify init of base-class element.Element by @adtzlr in https://github.com/adtzlr/felupe/pull/924
Fix hello_word(planestrain=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/928
Fix RegionQuadraticQuadBoundary by @adtzlr in https://github.com/adtzlr/felupe/pull/933
Don't create a figure if error in umat.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/934
Raise a TypeError if a container is used in a region by @adtzlr in https://github.com/adtzlr/felupe/pull/939
Add MeshContainer.as_vertex_mesh() and Field.from_mesh_container() by @adtzlr in https://github.com/adtzlr/felupe/pull/940
Fix FieldContainer.plot() for RegionVertex (no gradient available) by @adtzlr in https://github.com/adtzlr/felupe/pull/942
Add GaussLobatto quadrature scheme by @adtzlr in https://github.com/adtzlr/felupe/pull/943
Enhance mesh.interpolate_line(axis=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/945
Set version tag to 9.2.0 by @adtzlr in https://github.com/adtzlr/felupe/pull/946

Full Changelog: https://github.com/adtzlr/felupe/compare/v9.1.0...v9.2.0

- Python
Published by adtzlr about 1 year ago

[9.1.0] - 2024-11-23

[!TIP] This release adds JAX as alternative backend to the default tensortrax-backend for automatic differentiation. Hyperelastic and Material classes are available at felupe.constutition.jax, but only some models are included for JAX. JAX is great for computationally expensive models. However, it runs on single-precision by default. Hence, either the the tolerance of the Newton-Rhapson method must be relaxed or JAX must be enforced to use 64bit floats.

Added

Add the hessian of the element shape functions for a quadratic quad element QuadraticQuad.hessian().
Add the order-argument to FieldContainer.extract(order="C") as well as for Field, FieldAxisymmetric, FieldPlaneStrain to return C-contiguous arrays by default.
Add an optional multiplier to Laplace(multiplier=1.0).
Add optional keyword-arguments to math.transpose(**kwargs) to support optional out and order-keywords.
Add the attribute RegionBoundary.tangents, which contains a list of tangent unit vectors. For quad cell-types the length of this list is one and for hexahedron cell-types it is of length two.
Add math.inplane(A, vectors) to return the in-plane components of a symmetric tensor A, where the plane is defined by its standard unit vectors.
Add constitution.jax.Hyperelastic as a feature-equivalent alternative to Hyperelastic with jax as backend.
Add constitution.jax.Material(..., jacobian=None) with JAX as backend. A custom jacobian-callable may be passed to switch between forward- and backward-mode automatic differentiation.
Add material models for JAX-based materials. Hyperelastic models available at constitution.jax.models.hyperelastic: extended_tube(), miehe_goektepe_lulei(), mooney_rivlin(), neo_hooke(), yeoh(), third_order_deformation(), van_der_waals(). Lagrange (stress-based) models available at constitution.jax.models.lagrange: morph(), morph_representative_directions().
Add constitution.jax.total_lagrange(), constitution.jax.updated_lagrange() and constitution.jax.isochoric_volumetric_split() function decorators for the JAX hyperelastic material class.
Add an optional keyword-argument Region.astype(copy=True) to modify the data types of the arrays of the region in-place if copy=False.
Add FieldContainer.evaluate.right_cauchy_green_deformation() and math.right_cauchy_green_deformation(field) to evaluate the right Cauchy-Green deformation tensor.
Add math.strain(field, C=None, ..., **kwargs) to use a given right Cauchy-Green deformation tensor for the evaluation of the strain tensor.
Add the hyperelastic foam models storakers() and blatz_ko() for both AD-backends.
Add the saint_venant_kirchhoff_orthotropic() hyperelastic model (tensortrax).

Changed

Change default np.einsum(..., order="K") to np.einsum(..., order="C") in the methods of Field, FieldAxisymmetric, FieldPlaneStrain and FieldContainer.
Change supported Python versions to 3.9 - 3.12.
Change the dtype-argument in Region.astype(dtype) from an optional to a required argument.
Make tensortrax an optional dependency (again). Now FElupe does only depend on NumPy and SciPy, all other extras are optional.

Fixed

Fix the number of points for non-disconnected dual meshes. This reduces the assembled (sparse) vector- and matrix-shapes, which are defined on mixed-fields.
Fix wrong results of ConstitutiveMaterial.plot() if any stretch is non-physical, i.e. lower or equal zero. This raises an error now.

Removed

Remove the unused, deprecated aliases UserMaterial, UserMaterialStrain and UserMaterialHyperelastic.

What's Changed

Add QuadraticQuad.hessian() by @adtzlr in https://github.com/adtzlr/felupe/pull/851
Update issue templates by @adtzlr in https://github.com/adtzlr/felupe/pull/855
Fix pyOpenSci basic checks by @adtzlr in https://github.com/adtzlr/felupe/pull/854
Extract C-contiguous arrays by default from fields by @adtzlr in https://github.com/adtzlr/felupe/pull/856
Add an optional multiplier in Laplace(multiplier=1.0) by @adtzlr in https://github.com/adtzlr/felupe/pull/864
Add optional kwargs in math.transpose(**kwargs) by @adtzlr in https://github.com/adtzlr/felupe/pull/865
Fix the number of points for non-disconnected dual-meshes by @adtzlr in https://github.com/adtzlr/felupe/pull/867
Fix wrong result of ConstitutiveMaterial.plot() if any stretch <= 0 by @adtzlr in https://github.com/adtzlr/felupe/pull/872
Add RegionBoundary.tangents and math.inplane() by @adtzlr in https://github.com/adtzlr/felupe/pull/873
Add constitution.jax.Hyperelastic by @adtzlr in https://github.com/adtzlr/felupe/pull/876
Group the constitution module by autodiff-backends by @adtzlr in https://github.com/adtzlr/felupe/pull/878
Add felupe.constitution.autodiff.jax.Material by @adtzlr in https://github.com/adtzlr/felupe/pull/879
Add morph() for JAX-based Material by @adtzlr in https://github.com/adtzlr/felupe/pull/880
Simplify the constitution module by @adtzlr in https://github.com/adtzlr/felupe/pull/881
Add Total-Lagrange and Updated-Lagrange decorators for JAX by @adtzlr in https://github.com/adtzlr/felupe/pull/883
Add more JAX-based models by @adtzlr in https://github.com/adtzlr/felupe/pull/884
Add isochoric_volumetric_split() for JAX-hyperelastic material by @adtzlr in https://github.com/adtzlr/felupe/pull/885
Move constitution.jax.isochoric_volumetric_split() by @adtzlr in https://github.com/adtzlr/felupe/pull/887
Add optional in-place cast Region.astype(dtype, copy=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/888
Make the examples run faster on CI/CD by @adtzlr in https://github.com/adtzlr/felupe/pull/890
Add miehe_goektepe_lulei model for JAX by @adtzlr in https://github.com/adtzlr/felupe/pull/892
Make tensortrax optional by @adtzlr in https://github.com/adtzlr/felupe/pull/894
Test the import in an environment with no autodiff-backend present by @adtzlr in https://github.com/adtzlr/felupe/pull/895
Enable jax.Hyperelastic(..., parallel=True) with jax.pmap on quadrature-points axis by @adtzlr in https://github.com/adtzlr/felupe/pull/896
Update GitHub-Actions by @adtzlr in https://github.com/adtzlr/felupe/pull/897
Add math-function for the right Cauchy-Green deformation tensor by @adtzlr in https://github.com/adtzlr/felupe/pull/898
Add the Storåkers hyperelastic foam model by @adtzlr in https://github.com/adtzlr/felupe/pull/900
Add the Blatz-Ko hyperelastic foam model by @adtzlr in https://github.com/adtzlr/felupe/pull/901
Add the Van-der-Waals model for JAX by @adtzlr in https://github.com/adtzlr/felupe/pull/902
Add more hyperelastic models by @adtzlr in https://github.com/adtzlr/felupe/pull/903
Add constitution.lame_converter_orthotropic() by @adtzlr in https://github.com/adtzlr/felupe/pull/904
Add an optional strain exponent for the SVK-material by @adtzlr in https://github.com/adtzlr/felupe/pull/905

Full Changelog: https://github.com/adtzlr/felupe/compare/v9.0.0...v9.1.0

- Python
Published by adtzlr over 1 year ago

[9.0.0] - 2024-09-06

Added

Add Region.astype(dtype=None) to copy and cast the region arrays to a specified type.
Add Field(..., dtype=None) to cast the array with field values to a specified type.
Add Field.extract(dtype=None) to cast the extracted field gradient/interpolated values to a specified type.
Add hello_world() to print the lines of a minimal-working-example to the console.
Add mesh.interpolate_line(mesh, xi, axis) to interpolate a line mesh. The attribute points_derivative holds the derivatives of the independent variable w.r.t. the dependent variable(s). The column of the independent variable is filled with zeros.
Add optional keyword-argument SolidBody.assemble.matrix(block=True), also for SolidBody.assemble.vector(block=True) and the assemble-methods of SolidBodyNearlyIncompressible. If block=False, these methods will assemble a list of the upper-triangle sub block-vectors/-matrices instead of the combined block-vector/-matrix.
Add SolidBodyForce as a replacement for SolidBodyGravity with more general keyword arguments ("values" instead of "gravity", "scale" instead of "density").
Add Laplace to be used as user-material in a solid body. Works with scalar- and vector-valued fields.
Add an optional mechanics.Assemble(..., multiplier=None) argument which is used in external items like SolidBodyForce, SolidBodyGravity and PointLoad and is applied in newtonrhapson(items, ...).
Add a new submodule view which contains the View... classes like ViewSolid or ViewField, previously located at tools._plot.
Add the grad- and hess-arguments to the reload()- and copy()-methods of a Region, i.e. Region.reload(grad=None, hess=None).
Add LinearElasticOrthotropic.
Add SolidBodyCauchyStress in addition to SolidBodyPressure.
Add mesh.cell_types() which returns an object-array with cell-type mappings for FElupe and PyVista.
Add MeshContainer.from_unstructured_grid(grid, dim=None, **kwargs) to create a mesh-container from an unstructured grid (PyVista).

Changed

Change the internal initialization of field = Field(region, values=1, dtype=None) values from field.values = np.ones(shape) * values to field = np.full(shape, fill_value=values, dtype=dtype). This enforces field = Field(region, values=1) to return the gradient array with data-type int which was of type float before.
Initialize empty matrices of SolidBodyForce, SolidBodyGravity and PointLoad with dtype=float.
Don't multiply the assembled vectors of SolidBodyForce, SolidBodyGravity and PointLoad by -1.0. Instead, mechanics.Assemble(multiplier=-1.0) is used in the solid bodies.
Change the visibility of the internal helpers mechanics.Assemble, mechanics.Evaluate and mechanics.Results from private to public.
Import the assembly module to the global namespace.
Isolate the submodules, i.e. a submodule only uses the public API of another submodule. If necessary, this will help to change one or more modules to a future extension package.
Enforce contiguous arrays for the region shape-function and -gradient arrays h and dhdX. This recovers the integral-form assembly performance from v8.6.0.
Make the private basis classes public (assembly.expression.Basis, assembly.expression.BasisField and assembly.expression.BasisArray) as especially their docstrings are useful to understand how a Basis is created on a field.
Remove material parameter keyword-arguments in the function()-, gradient()- and hessian()-methods of the core constitutive material formulations. This removes the ability for temporary material parameters in LinearElastic(E=None, nu=0.3).gradient(x, E=1.0), use LinearElastic(E=1.0, nu=0.3).gradient(x) instead. This affects the material formulations LinearElastic, LinearElasticPlaneStrain, LinearElasticPlaneStress, LinearElasticLargeStrain, NeoHooke and NeoHookeCompressible. None is still supported for the material parameters of NeoHooke and NeoHookeCompressible.
Remove LinearElasticPlaneStrain from the top-level package namespace because this should only be used with Field(region, dim=2). The preferred method is to use FieldPlaneStrain(region, dim=2) and the default LinearElastic. LinearElasticPlaneStrain remains available in constitution.LinearElasticPlaneStrain.
Rename Mesh.as_pyvista() to Mesh.as_unstructured_grid() and add Mesh.as_pyvista() as alias.

Deprecated

Deprecate SolidBodyGravity, SolidBodyForce should be used instead.

What's Changed

Add Region.astype(dtype=None) and Field(..., dtype=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/797
Add hello_world() by @adtzlr in https://github.com/adtzlr/felupe/pull/799
Add mesh.interpolate_line(mesh, xi, axis) by @adtzlr in https://github.com/adtzlr/felupe/pull/811
Change the derivative storage in mesh.interpolate_line() by @adtzlr in https://github.com/adtzlr/felupe/pull/812
Add optional SolidBody.assemble.matrix(block=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/813
Add SolidBodyForce and Laplace by @adtzlr in https://github.com/adtzlr/felupe/pull/814
Add optional Assemble(multiplier=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/817
Change the visibility of the helpers in the mechanics-module from private to public by @adtzlr in https://github.com/adtzlr/felupe/pull/818
Fix Ex.09: Take care of changed sign of external form items by @adtzlr in https://github.com/adtzlr/felupe/pull/820
Add submodule view by @adtzlr in https://github.com/adtzlr/felupe/pull/823
Enhance the imports of the quadrature submodule by @adtzlr in https://github.com/adtzlr/felupe/pull/824
Enhance the imports of the mechanics submodule by @adtzlr in https://github.com/adtzlr/felupe/pull/825
Enhance (isolate) the cross-imports between submodules by @adtzlr in https://github.com/adtzlr/felupe/pull/826
Fix the assembly performance by @adtzlr in https://github.com/adtzlr/felupe/pull/828
Test if region attribute arrays are C-contiguous by @adtzlr in https://github.com/adtzlr/felupe/pull/829
Add an example for the third medium contact method by @adtzlr in https://github.com/adtzlr/felupe/pull/834
Make the classes Basis, BasisField and BasisArray public by @adtzlr in https://github.com/adtzlr/felupe/pull/835
Add hess argument to Region(hess=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/836
Add LinearElasticOrthotropic by @adtzlr in https://github.com/adtzlr/felupe/pull/840
Add SolidBodyCauchyStress by @adtzlr in https://github.com/adtzlr/felupe/pull/842
Remove temporary kwargs in core constitutive materials by @adtzlr in https://github.com/adtzlr/felupe/pull/843
Don't import LinearElasticPlaneStrain in top-level package namespace by @adtzlr in https://github.com/adtzlr/felupe/pull/844
Add felupe.mesh.cell_types() by @adtzlr in https://github.com/adtzlr/felupe/pull/846
Add MeshContainer.from_unstructured_grid() by @adtzlr in https://github.com/adtzlr/felupe/pull/847
Update README.md by @adtzlr in https://github.com/adtzlr/felupe/pull/850

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.8.0...v9.0.0

- Python
Published by adtzlr over 1 year ago

[8.8.0] - 2024-06-16

Added

Add Region(uniform=False), a flag to invoke a reduced evaluation of the element shape functions and their gradients. If True, only the first cell is evaluated. This drastically speeds up linear-elasticity on uniform grid meshes (e.g. voxel-like quads and hexahedrons).
Add an update-method in FormItem(..., ramp_item=0).update(value). This enables a FormItem to be used as a ramped-item in a Step.

What's Changed

Add Region(uniform=False): A performance feature for uniform grid meshes by @adtzlr in https://github.com/adtzlr/felupe/pull/792
Add FormItem().update(values) by @adtzlr in https://github.com/adtzlr/felupe/pull/796

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.7.0...v8.8.0

- Python
Published by adtzlr almost 2 years ago

[8.7.0] - 2024-06-07

[!NOTE]
This is the first release with support for NumPy 2.0.

Added

Add math.solve_nd(A, b, n=1) as a generalized function of math.solve_2d(A, b) with supported broadcasting on the elementwise-operating trailing axes.

Changed

Rebase math.solve_2d(A, b) on math.solve_nd(A, b, n=2) with a batched-rhs supported in NumPy 2.0.
Change ogden_roxburgh() and OgdenRoxburgh to use the Gauss error function erf instead of tanh as internal sigmoid function.
Flatten the returned inverse indices-array of np.unique(..., return_inverse=True) in mesh.merge_duplicate_points() to ensure compatibility with NumPy 2.0.

Fixes

Reset state variables in PlotMaterial.evaluate() after each completed load case.

What's Changed

Reset state variables in PlotMaterial.evaluate() after each completed load case by @adtzlr in https://github.com/adtzlr/felupe/pull/784
Add math.solve_nd(A, b, n=1) by @adtzlr in https://github.com/adtzlr/felupe/pull/785
Support broadcasting in matrix-axes of math.solve_nd() by @adtzlr in https://github.com/adtzlr/felupe/pull/787
Fix small typo focusses -> focuses by @tkoyama010 in https://github.com/adtzlr/felupe/pull/786
Change ogden_roxburgh() and OgdenRoxburgh to use erf instead of tanh by @adtzlr in https://github.com/adtzlr/felupe/pull/788
Add support for NumPy 2.0 by @adtzlr in https://github.com/adtzlr/felupe/pull/789

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.6.0...v8.7.0

- Python
Published by adtzlr almost 2 years ago

[8.6.0] - 2024-05-29

Added

Add the isotropic-hyperelastic alexander(C1, C2, C2, gamma, k) material model formulation to be used in Hyperelastic().
Add the isotropic-hyperelastic Micro-Sphere miehe_goektepe_lulei(mu, N, U, p, q) material model formulation to be used in Hyperelastic().
Add the pseudo-elastic ogden_roxburgh(r, m, beta, material, **kwargs) material model formulation to be used in Hyperelastic().
Add an optional relative-residuals argument to ConstitutiveMaterial.optimize(relative=False).
Add a class-decorator constitutive_material(Msterial, name=None).
Add the Total-Lagrange MORPH material formulation implemented by the concept of representative directions morph_representative_directions(p) to be used in MaterialAD().
Add the Total-Lagrange (original) MORPH material formulation morph(p) to be used in Hyperelastic().
Add decorators @total_lagrange and @updated_lagrange for Total / Updated Lagrange material formulations to be used in MaterialAD.
Add the isotropic-hyperelastic anssari_benam_bucchi(mu, N) material model formulation to be used in Hyperelastic().
Add the isotropic-hyperelastic lopez_pamies(mu, alpha) material model formulation to be used in Hyperelastic().

Changed

Recfactor the constitution module.

Fixed

Fix plotting the keyword-arguments of a constitutive material ConstitutiveMaterial.plot(show_kwargs=True). For list-based material parameters of length 1, the brackets aren't shown now. E.g., this affects optimized material parameters.
Don't update the material parameters in-place in ConstitutiveMaterial.optimize().
Don't convert material parameter scalars to arrays in ConstitutiveMaterial.optimize().

What's Changed

Add the isotropic-hyperelastic Alexander model by @adtzlr in https://github.com/adtzlr/felupe/pull/756
Re-organize the constitution module by @adtzlr in https://github.com/adtzlr/felupe/pull/758
Delete modelshyperelasticity_ad.py by @adtzlr in https://github.com/adtzlr/felupe/pull/759
Fix the Alexander model by @adtzlr in https://github.com/adtzlr/felupe/pull/761
Add ogden_roxburgh using AD for Hyperelastic by @adtzlr in https://github.com/adtzlr/felupe/pull/762
Add ConstitutiveMaterial.optimize(relative=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/764
Add the Micro-Sphere hyperelastic material model framework by @adtzlr in https://github.com/adtzlr/felupe/pull/766
Add constitutive_material(Material) decorator by @adtzlr in https://github.com/adtzlr/felupe/pull/768
Add the MORPH material model (by the method of representative directions) by @adtzlr in https://github.com/adtzlr/felupe/pull/770
Add the MORPH material model formulation by @adtzlr in https://github.com/adtzlr/felupe/pull/771
Add @total_lagrange and @updated_lagrange decorators for MaterialAD by @adtzlr in https://github.com/adtzlr/felupe/pull/773
Rebase morph_representative_directions() as Total-Lagrange material for MaterialAD by @adtzlr in https://github.com/adtzlr/felupe/pull/774
Add the Anssari-Benam Bucchi isotropic-hyperelastic material model formulation by @adtzlr in https://github.com/adtzlr/felupe/pull/775
Don't update the material parameters in-place and don't convert them to arrays in ConstitutiveMaterial.optimize() by @adtzlr in https://github.com/adtzlr/felupe/pull/778
Add the Lopez-Pamies hyperelastic model lopez_pamies(mu, alpha) by @adtzlr in https://github.com/adtzlr/felupe/pull/779
Fix small typo cartesion -> cartesian by @tkoyama010 in https://github.com/adtzlr/felupe/pull/781
add example for nonlinear viscoelasticity by @bhaveshshrimali in https://github.com/adtzlr/felupe/pull/780

New Contributors

@tkoyama010 made their first contribution in https://github.com/adtzlr/felupe/pull/781
@bhaveshshrimali made their first contribution in https://github.com/adtzlr/felupe/pull/780

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.5.1...v8.6.0

- Python
Published by adtzlr about 2 years ago

[8.5.1] - 2024-05-08

Fixed

Fix the gradient evaluation in NeoHookeCompressible(mu=1, lmbda=None).
Fix MaterialStrain.plot().

What's Changed

Fix gradient-evaluation in NeoHookeCompressible(lmbda=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/754
Fix MaterialStrain.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/755

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.5.0...v8.5.1

- Python
Published by adtzlr about 2 years ago

[8.5.0] - 2024-04-27

Added

Add umat_new, res = ConstitutiveMaterial.optimize(ux=[stretches, stresses], ps=None, bx=None) to optimize the material parameters on given experimental data for incompressible uniaxial, biaxial and / or planar tension / compression data by scipy.optimize.least_squares().
Add initial default material parameters for the hyperelastic material model formulations.

What's Changed

Add ConstitutiveMaterial.optimize() by @adtzlr in https://github.com/adtzlr/felupe/pull/750

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.4.0...v8.5.0

- Python
Published by adtzlr about 2 years ago

[8.4.0] - 2024-04-12

Added

Add math.solve_2d(A, b, solver=np.linalg.solve, **kwargs) to be used in newtonrhapson(solve=solve_2d, ...) for two-dimensional unknowns. This is useful for local Newton-iterations related to viscoelastic evolution equations inside constitutive material formulations.
Add x- and y-offsets in Job.plot(xoffset=0.0, yoffset=0.0).

Changed

Wrap the ax-title with the parameters of the material model in ConstitutiveMaterial.plot().

Fixed

Sort array of principal values in descending order before plotting in Scene.plot("Principal Values of ..."). This ensures that the labels are matching user-defined arrays of principal values.

What's Changed

Sort array of principal values in Scene.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/739
Add math.solve_2d(A, b) by @adtzlr in https://github.com/adtzlr/felupe/pull/742
Wrap the ax-title ConstitutiveMaterial.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/743

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.3.1...v8.4.0

- Python
Published by adtzlr about 2 years ago

[8.3.1] - 2024-04-06

Fixed

Set the default verbosity level to None in newtonrhapson(verbose=None) and Job.evaluate(verbose=None). If None, this defaults to True (as before) but evaluates the environmental variable FELUPE_VERBOSE if present with FELUPE_VERBOSE == "true". This does not ignore custom verbosity levels.

What's Changed

Change default verbose=None in Job.evaluate() and newtonrhapson() by @adtzlr in https://github.com/adtzlr/felupe/pull/736

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.3.0...v8.3.1

- Python
Published by adtzlr about 2 years ago

[8.3.0] - 2024-04-02

Added

Add a method to convert a mesh to a PyVista unstructured grid Mesh.as_pyvista(cell_type=None).

Changed

Change default line-width from 1.0 to 2.0 in ViewMesh.plot(line_width=2.0).
Enforce a 3d points-array in Mesh.as_meshio().

Fixed

Fix previously ignored line-width in the plot of a mesh ViewMesh.plot(line_width=1.0).
Fix math.tovoigt() for one-dimensional tensors.

What's Changed

Add Mesh.as_pyvista(), fix and change default ViewMesh.plot(line_width=1.0) by @adtzlr in https://github.com/adtzlr/felupe/pull/733

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.2.1...v8.3.0

- Python
Published by adtzlr about 2 years ago

[8.2.1] - 2024-03-30

Fixed

Fix FieldContainer.evaluate.strain(fun=lambda stretch: stretch) for custom strain-stretch callables. The fun-argument was previously ignored.

What's Changed

Add field.EvaluateFieldContainer(field) by @adtzlr in https://github.com/adtzlr/felupe/pull/728
Outsource EvaluateFieldContainer to field/_evaluate.py by @adtzlr in https://github.com/adtzlr/felupe/pull/729
Fix FieldContainer.evaluate.strain(fun=lambda stretch: stretch) for custom fun-callables by @adtzlr in https://github.com/adtzlr/felupe/pull/731

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.2.0...v8.2.1

- Python
Published by adtzlr about 2 years ago

[8.2.0] - 2024-03-25

Added

Add methods to evaluate different strain measures of a field container, i.e. FieldContainer.evaluate.strain(tensor=True, asvoigt=False, k=0), FieldContainer.evaluate.log_strain(tensor=True, asvoigt=False) or FieldContainer.evaluate.green_lagrange_strain(tensor=True, asvoigt=False). These methods refer to math.strain(k=0) which uses math.strain_stretch_1d(k=0) by default with the strain exponent k of the Seth-Hill strain formulation.

Changed

Change the return type of math.eig() which returns a namedtuple with attributes eigenvalues and eigenvectors. Now consistent with NumPy. This also affects math.eigh().
Change the shape of the array of eigenvectors returned by math.eig() from (a, i, ...) to (i, a, ...). Now consistent with NumPy. This also affects math.eigh().
Switch from a hard-coded logarithmic strain evaluation to a more generalized Seth-Hill strain-stretch relation in math.strain(fun=math.strain_stretch_1d, k=0) with a default strain-exponent of k=0 which returns the logarithmic strain.

What's Changed

Enhance plots of quadrature schemes by @adtzlr in https://github.com/adtzlr/felupe/pull/723
Enhance the plots of GaussLegendreBoundary by @adtzlr in https://github.com/adtzlr/felupe/pull/724
Add FieldContainer.evaluate.strain() by @adtzlr in https://github.com/adtzlr/felupe/pull/726

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.1.0...v8.2.0

- Python
Published by adtzlr about 2 years ago

[8.1.0] - 2024-03-23

Added

Add an argument to disable the (default) expansion of the points-array of a mesh in mesh.expand(expand_dim=True) and mesh.revolve(expand_dim=True). E.g., this allows the expansion and / or revolution of a quad-mesh with points in 3d-space.
Add MultiPointContact.plot(offset=0, **kwargs) to plot the rigid contact plane(s) or line(s) at a given offset.
Add MultiPointConstraint.plot(**kwargs) to plot the lines of a multi-point constraint.

Changed

Don't raise an error if the total angle of revolution is greater than 360 degree in mesh.revolve(phi=361).

What's Changed

Don't raise an error for mesh-revolutions with >360 degree, e.g. mesh.revolve(phi=361) by @adtzlr in https://github.com/adtzlr/felupe/pull/715
Enhance mesh.expand(expand_dim=True) and mesh.revolve(expand_dim=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/716
Add MultiPointContact.plot(**kwargs) method by @adtzlr in https://github.com/adtzlr/felupe/pull/718
Add MultiPointConstraint.plot(**kwargs) by @adtzlr in https://github.com/adtzlr/felupe/pull/719
Enhance the plot()-methods of region, element and quadrature by @adtzlr in https://github.com/adtzlr/felupe/pull/721
Enhance element plots by @adtzlr in https://github.com/adtzlr/felupe/pull/722

Full Changelog: https://github.com/adtzlr/felupe/compare/v8.0.0...v8.1.0

- Python
Published by adtzlr about 2 years ago

[8.0.0] - 2024-03-18

Added

Add axis of expansion in mesh.expand(axis=-1) (ignored for n=1).
Add an optional mask-argument to select points for rotation in mesh.rotate(mask=None).
Add Lagrange quad/hex cell-types in ViewMesh.
Add optional projection of stresses from quadrature-points to mesh.points in SolidBody.plot(project=None), where project has to be a callable like project(values, region).
Add optional projection of internel cell-data (Displacement, Logarithmic Strain and Deformation Gradient) from quadrature-points to mesh-points in FieldContainer.plot(project=None), where project has to be a callable like project(values, region).

Changed

The internal BasisField.basis is now a subclassed array BasisArray with a grad-attribute.
math.grad(x, **kwargs) is enhanced to return gradients of fields (like before) and the gradient-attribute of basis-arrays (added).
The grad_v and grad_u arguments are removed from the form-expression decorator Form. This changes the required function signature of the weakform-callable to weakform(v, u, **kwargs). The tuple of optional arguments is also removed. Gradients of v and u are now obtained by math.grad(v) or v.grad.
Enforce quadrature schemes with minimal order for projections in project() for Triangle, Tetra as well as their MINI- and Quadratic-variants.
Fall-back to extrapolate(mean=True) in project(mean=True).
Don't ravel the results of res = extrapolate(values, region), i.e. values.shape = (3, 3, 4, 100) will be returned as res.shape = (121, 3, 3) instead of res.shape = (121, 9).
Stack only a selection of meshes in MeshContainer.stack([idx]).
Enable list-based indexing in MeshContainer[idx].
Add the opacity=0.99 argument to MeshContainer.plot() and MeshContainer.screenshot().
Pass the dpi-argument to the matplotlib figure in imshow(dpi=None) for solids, field- and mesh-containers.
Permute GaussLegendre(order=2, dim=2) according to the points of the BiQuadraticQuad element by default.
Permute the 2- and 3-dimensional GaussLegendre quadrature schemes for order > 2 according to the VTK-Lagrange element formulations. That means for linear and quadratic quads and hexahedrons, the points of GaussLegendre are sorted according to the default VTK elements and for all higher-order elements according to the Lagrange-elements.
Enable default point-permutations in RegionLagrange(permute=True) by default.
Hide internal edges of higher-order cell-types in ViewScene.plot() by default.
Simplify tools.topoints(values, region, average=True, mean=False). Remove all other arguments. If values of single quadrature-point per cells is given, then the values are broadcasted to the number of points-per-cell. If values are provided on more quadrature points than the number of points-per-cell, then the values are trimmed. E.g., this is required for QuadraticHexahedron with 20 points and 27 quadrature-points.

Fixed

Fix mesh-expansion with one layer mesh.expand(n=1). This expands the dimension of the points-array.
Fix VTK-compatible cells in CubeArbitraryOrderHexahedron.
Fix Cauchy-stress evaluation of SolidBody and SolidBodyNearlyIncompressible on a 2d-Field (plane stress): Automatic fall-back to Kirchhoff-stress and print a warning.

Removed

Remove the deprecated old-style argument move in dof.biaxial().
Remove the deprecated old-style arguments move, axis_compression, axis_shear and compression in dof.shear().

What's Changed

Simplify Form: Remove required grad_v and grad_u args and obtain the gradient by v.grad or grad(v) by @adtzlr in https://github.com/adtzlr/felupe/pull/688
Remove deprecated arguments in dof.biaxial() and dof.shear() by @adtzlr in https://github.com/adtzlr/felupe/pull/689
Fall-back to tools.extrapolate(mean=True) for project(mean=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/691
Fix single-layer mesh.expand(n=1) and add axis of expansion mesh.expand(axis=-1) by @adtzlr in https://github.com/adtzlr/felupe/pull/697
Enhance imshow()-methods of solid, field and mesh by @adtzlr in https://github.com/adtzlr/felupe/pull/699
Add optional mask in mesh.rotate(mask=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/700
Make ArbitraryOrderLagrangeElement VTK-compatible by @adtzlr in https://github.com/adtzlr/felupe/pull/702
Enhance Scene.plot(): Hide internal edges of quadratic celltypes by @adtzlr in https://github.com/adtzlr/felupe/pull/705
Add optional projection of stresses in SolidBody.plot(project=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/707
Add FieldContainer.plot(project=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/710
Simplify tools.topoints(values, region) by @adtzlr in https://github.com/adtzlr/felupe/pull/712
Enhance topoints( values, region, average=True, mean=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/713
Fix Cauchy-stress for 2d-Fields (plane stress) by @adtzlr in https://github.com/adtzlr/felupe/pull/714

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.19.1...v8.0.0

- Python
Published by adtzlr about 2 years ago

[7.19.1] - 2024-03-08

Fixed

Fix tools.project() for meshes where some points are not connected to cells.

What's Changed

Fix project() for meshes with points without cells by @adtzlr in https://github.com/adtzlr/felupe/pull/687

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.19.0...v7.19.1

- Python
Published by adtzlr about 2 years ago

[7.19.0] - 2024-03-08

Added

Add FieldDual(disconnect=True) for a dual (secondary) field with an optionally disconnected mesh. This also enables FieldsMixed(disconnect=True) in mixed fields.
Add a quadrature scheme for integrating the surface of a unit hemisphere BazantOh(n=21).
Add NearlyIncompressible as a simplified version of ThreeFieldVariation. A constitutive material formulation on the distortional part of a strain energy function in terms of the deformation gradient has to be provided, e.g. by umat = NearlyIncompressible(NeoHooke(mu=1), bulk=5000).
Add optional kwargs to a job-callback Job(callback=lambda stepnumber, substepnumber, substep, **kwargs: None, **kwargs) and CharacteristicCurve(callback=lambda stepnumber, substepnumber, substep, **kwargs: None, **kwargs).
Add DiscreteGeometry properties x, y and z to access the columns of the points-array.
Add a new math-function math.equivalent_von_mises(A) for three-dimensional second-order tensors.
Add the evaluation of the equivalent von Mises Cauchy stress as cell-data in ViewSolid, available as Solid.plot("Equivalent of Cauchy Stress").
Add mesh.stack(meshes) as method to MeshContainer.stack(). Note that this only supports mesh containers with meshes of same cell-types.
Add NeoHooke.gradient(out=None) and NeoHooke.hessian(out=None) for a location to store the results. Also for NeoHookeCompressible.
Add out-keyword to gradient() and hessian of NearlyIncompressible and ThreeFieldVariation.
Add optional initial state variables in ViewMaterial(statevars=None) and ViewMaterialIncompressible(statevars=None).
Add the L2-projection as tools.project(values, region, average=True, mean=False, dV=None, solver=scipy.sparse.linalg.spsolve) to project given values at quadrature points to mesh-points. This replaces the old tools.project(values, region, average=True, mean=False) in a backward-compatible way. The new method is computationally more expensive but is also much more flexible.
Add fifth-order quadrature schemes quadrature.Triangle(order=5) and quadrature.Tetrahedron(order=5).
Add Region.copy(mesh=None, element=None, quadrature=None) to copy a region and re-evaluate this copy if necessary.

Changed

Rename Mesh.save() to Mesh.write() and add Mesh.save() as an alias to Mesh.write().
Enhance the performance of NeoHooke, NeoHookeCompressible, SolidBody and SolidBodyNearlyIncompressible.
Enhance the performance of math.inv(out=None) and math.det(out=None).
Use only the offical API of tensortrax. A workaround is used to ensure compatibility with tensortrax <= v0.17.1.
Pass optional keyword-arguments in the plot-methods ViewMaterial.plot(**kwargs) and ViewMaterialIncompressible.plot(**kwargs) to the matplotlib axes object ax.plot(**kwargs).
Only add off_screen and notebook keyword-arguments to pyvista.Plotter(**kwargs) if they are True. This is needed for not ignoring a global variable like pyvista.OFF_SCREEN = True.
Enforce verbose=0 if the environmental variable "FELUPE_VERBOSE" is "false". This is useful for running the examples when building the documentation.
Don't require a bilinearform in FormItem(bilinearform=None). An empty FormItem is now a valid item in a Step. For empty vectors/matrices, the shape is inferred from sum(FieldContainer.fieldsizes) instead of FieldContainer.fields[0].values.size.
Rename the old-project method to tools.extrapolate(values, region, average=True, mean=False) which extrapolates values at quadrature points to mesh-points.
Change the sorting of quadrature points for triangles and tetrahedrons (due to internal code simplifications).
The reload-method of a region does only re-evaluate it if at least one of the arguments are not None Region.reload(mesh, element, quadrature).

Fixed

Fix missing support for third-order- and second-order tensor combinations to math.dot(A, B, mode=(2,3)) and math.ddot(A, B, mode=(2,3)).
Fix error if FieldDual is in the fields of a FieldContainer for IntegralForm.
Fix math.inv(A) for arrays with shape A.shape = (1, 1, ...). Also raise an error if shape[:2] not in [(3, 3), (2, 2), (1, 1)].
Raise an error in math.det(A) if A.shape[:2] not in [(3, 3), (2, 2), (1, 1)].
Fix mutable keyword-arguments in SolidBody._vector(kwargs={}) by SolidBody._vector(kwargs=None). Also for ._matrix() and for SolidBodyNearlyIncompressible.
Fix wrong shape and the resulting error during assembly in fem.assembly.expression.Form for the integration of a linear form with different mesh- and field-dimensions.

What's Changed

Rename Mesh.save() to Mesh.write() and add Mesh.save() as an alias of it by @adtzlr in https://github.com/adtzlr/felupe/pull/626
Refactor constitution by @adtzlr in https://github.com/adtzlr/felupe/pull/628
Add FieldDual for the dual (secondary) fields in mixed fields by @adtzlr in https://github.com/adtzlr/felupe/pull/629
Add quadrature scheme BazantOh() by @adtzlr in https://github.com/adtzlr/felupe/pull/632
Add NearlyIncompressible(material, bulk) by @adtzlr in https://github.com/adtzlr/felupe/pull/633
Add math.dot(A, B, mode=(2,3)) and math.ddot(A, B, mode=(2,3)) by @adtzlr in https://github.com/adtzlr/felupe/pull/641
Add optional kwargs to Job(callback=lambda i, j, res, **kwargs: None, **kwargs) by @adtzlr in https://github.com/adtzlr/felupe/pull/642
Add support for FieldDual in IntegralForm by @adtzlr in https://github.com/adtzlr/felupe/pull/643
Fix math.inv(A) for A.shape[:2] == (1, 1) by @adtzlr in https://github.com/adtzlr/felupe/pull/644
Add FieldDual(disconnect=None). by @adtzlr in https://github.com/adtzlr/felupe/pull/646
Add DiscreteGeometry properties x, y and z by @adtzlr in https://github.com/adtzlr/felupe/pull/648
Add math.equivalent_von_mises(A) by @adtzlr in https://github.com/adtzlr/felupe/pull/650
Add MeshContainer.stack() by @adtzlr in https://github.com/adtzlr/felupe/pull/652
Enhance the performance by @adtzlr in https://github.com/adtzlr/felupe/pull/654
Enhance the performance of `math.inv(out=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/655
Improve the performance of math.det(A, out=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/659
Use only the official API of tensortrax.math by @adtzlr in https://github.com/adtzlr/felupe/pull/667
Only add off_screen and notebook keyword-arguments to pyvista.Plotter(**kwargs) if they are True by @adtzlr in https://github.com/adtzlr/felupe/pull/670
Docs: Examples are generated by Sphinx-Gallery by @adtzlr in https://github.com/adtzlr/felupe/pull/671
Fix error during assembly in Form for different mesh- and field-dimensions by @adtzlr in https://github.com/adtzlr/felupe/pull/675
Add the L2-projection and replace it with the old tools.project() by @adtzlr in https://github.com/adtzlr/felupe/pull/683

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.18.0...v7.19.0

- Python
Published by adtzlr about 2 years ago

[7.18.0] - 2024-02-16

Added

Create a FieldContainer by the &-operator between fields and field containers, i.e. field = displacement & pressure, where displacement = Field(region, dim=2) and pressure = Field(region). This also works for field & pressure as well as pressure & field.
Add a method to create a field container from a field, i.e. Field(region, dim=3).as_container() is equal to FieldContainer([Field(region, dim=3)]).
Add ViewMaterial(umat) to view force-stretch curves for uniaxial tension/compression, planar shear and equi-biaxial tension.
Add ViewMaterialIncompressible(umat) to view force-stretch curves for incompressible uniaxial tension/compression, planar shear and equi-biaxial tension.
Add a base class for constitutive materials with methods ConstitutiveMaterial.view(incompressible=False), ConstitutiveMaterial.plot(incompressible=False) and ConstitutiveMaterial.screenshot(incompressible=False).
Add a dict-attribute with material parameters to all built-in materials, e.g. NeoHooke.kwargs = {"mu": self.mu, "bulk": self.bulk}.
Add umat = CompositeMaterial(material, other_material).
Add &-operator to combine constitutive materials umat = material & other_material. Note that only the first material must contain state variables.

Changed

Don't disconnect the dual mesh by default for regions RegionQuadraticTriangle and RegionQuadraticTetra in FieldsMixed.

Fixed

Fix Mesh.flip(mask=None): Take care of the mask (it wasn't applied to the cells-array of the mesh).

What's Changed

Create FieldContainer by all logical-and combinations of Field and FieldContainer by @adtzlr in https://github.com/adtzlr/felupe/pull/618
Add Field.as_container() by @adtzlr in https://github.com/adtzlr/felupe/pull/619
Add ViewMaterialIncompressible by @adtzlr in https://github.com/adtzlr/felupe/pull/620
Add base class ConstitutiveMaterial by @adtzlr in https://github.com/adtzlr/felupe/pull/621
Add Volumetric, CompositeMaterial and &-operator for ConstitutiveMaterial by @adtzlr in https://github.com/adtzlr/felupe/pull/623
Fix Mesh.flip(mask=None): Take care of the mask by @adtzlr in https://github.com/adtzlr/felupe/pull/625

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.17...v7.18.0

- Python
Published by adtzlr over 2 years ago

[7.17.0] - 2024-02-15

Added

Add a mesh for a single vertex point vertex = Point(a=0).
Add expansion of a vertex point to a line mesh vertex.expand(n=11, z=1).
Add revolution of a vertex point to a line mesh vertex.revolve(n=11, phi=180).

Changed

Assume that no state variables are used in an umat if it has no attribute umat.x. Set the shape of the state variables by default to (0, q, c) in SolidBody and SolidBodyNearlyIncompressible.

What's Changed

Add mesh.Point(a=0) by @adtzlr in https://github.com/adtzlr/felupe/pull/616
Assume no state variables in umat if it has no attribute x by @adtzlr in https://github.com/adtzlr/felupe/pull/617

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.16...v7.17

- Python
Published by adtzlr over 2 years ago

[7.16.0] - 2024-02-13

Added

Add MeshContainer.plot(), img = MeshContainer.screenshot() and ax = MeshContainer.imshow(). The default list of colors contains PyVista's default color as first item and then the list of matplotlib's named colors C1, C2, etc (excluding C0).
Add Mesh.merge_duplicate_points(decimals=None) and make Mesh.sweep(decimals=None) an alias of it.
Add Mesh.merge_duplicate_cells().

Fixed

Fix Mesh.imshow(ax=None), FieldContainer.imshow(ax=None) and SolidBody.imshow(ax=None) from v7.15.

What's Changed

Add MeshContainer.plot(), img = MeshContainer.screenshot() and ax = MeshContainer.imshow() by @adtzlr in https://github.com/adtzlr/felupe/pull/611
Add mesh.merge_duplicate_points() and mesh.merge_duplicate_cells() by @adtzlr in https://github.com/adtzlr/felupe/pull/612
Fix Mesh.imshow(ax=None), FieldContainer.imshow(ax=None) and SolidBody.imshow(ax=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/613

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.15.0...v7.16

- Python
Published by adtzlr over 2 years ago

[7.15.0] - 2024-02-11

Added

Add optional ax-arguments to Mesh.imshow(ax=None), FieldContainer.imshow(ax=None) and SolidBody.imshow(ax=None).

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.14.0...v7.15.0

- Python
Published by adtzlr over 2 years ago

[7.14.0] - 2024-02-11

Added

Add optional output location for FieldContainer.extract(out=None).
Add Mesh.update(callback=None). This is especially useful if the points array of a mesh is changed and an already existing instance of a region has to be reloaded: Mesh.update(points=new_points, callback=region.reload).
Add ax = FieldContainer.imshow() which acts as a wrapper on top of the img = FieldContainer.screenshot(filename=None) method. The image data is passed to a matplotlib figure and the ax object is returned.
Add ax = Mesh.imshow() which acts as a wrapper on top of the img = Mesh.screenshot(filename=None) method. The image data is passed to a matplotlib figure and the ax object is returned.
Add view-methods for SolidBody and SolidBodyNearlyIncompressible (same as already implemented for mesh and fields).
Add lists with norms of values and norms of the objective function in tools.NewtonResult(xnorms=None, fnorms=None).
Add lists of norms of the objective function as attribute to Job.fnorms.
Add new method to Mesh for getting point ids next to a given point coordinate Mesh.get_point_ids(value).
Add new method to Mesh for getting cells attached to a given point coordinate Mesh.get_cell_ids(point_ids).
Add new method to Mesh for getting neighbour cells Mesh.get_cell_ids_neighbours(cell_ids).
Add new method to Mesh for getting shared points between neighbour cells Mesh.get_point_ids_shared(cell_ids_neighbours).
Add new method to Mesh for getting points on regular grids which are located at corners Mesh.get_point_ids_corners().
Add new method to Mesh for modifying the cell connectivity at corners Mesh.modify_corners(point_ids=None), supported for regular quad and hexahedron meshes.

Changed

Pass optional keyword-arguments in math.dot(**kwargs) to the underlying einsum-calls.
Enhance the performance of math.cdya() by reducing the number of (intermediate) arrays to be created from 4 to 2.
Use fixed output locations for the extracted field-gradients and the integrated stiffness matrices in SolidBody and SolidBodyNearlyIncompressible. This enhances the performance.
Change default filename in Mesh.screenshot() from filename="field.png" to filename="mesh.png".
Change the return value on job-evaluation from None = Job.evaluate() to job = Job.evaluate().
Change implementation of LinearElasticLargeStrain from NeoHooke to NeoHookeCompressible.
Do not invoke pyvista.start_xvfb() on a posix-os. If required, run it manually.
Rename tools._newton.Result to tools._newton.NewtonResult and add it to the public API as tools.NewtonResult because this class is returned as a result of Newton's method.
Rename the r-arguments of tools.force() and tools.moment() to forces.
Rename the point-argument of tools.moment() to centerpoint.
Rename the r-argument in tools.save() to forces. Remove the unused argument converged.
Change the default file-extension from .vtk to .vtu in tools.save(filename="result.vtu").
Change the default values of the gravity vector to zeros if gravity=None in SolidBodyGravity(field, gravity=None, density=1.0).

Fixed

Fix tools.moment(). Use math.cross(). The old implementation was completely wrong!

What's Changed

Enhance performance of cdya() and pass kwargs in dot() by @adtzlr in https://github.com/adtzlr/felupe/pull/593
Add output location for field.extract(out=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/594
Add out-locations to the integral forms and use them in the solid bodies by @adtzlr in https://github.com/adtzlr/felupe/pull/595
Add Mesh.update(callback=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/596
Do not invoke pyvista.start_xvfb() on a posix-os by @adtzlr in https://github.com/adtzlr/felupe/pull/597
Add Mesh.imshow() and FieldContainer.imshow() by @adtzlr in https://github.com/adtzlr/felupe/pull/598
Change return value in job = Job.evaluate() by @adtzlr in https://github.com/adtzlr/felupe/pull/599
Add view-methods for solid bodies by @adtzlr in https://github.com/adtzlr/felupe/pull/600
Change implementation of LinearElasticLargeStrain from NeoHooke to NeoHookeCompressible by @adtzlr in https://github.com/adtzlr/felupe/pull/601
Rename tools._newton.Result to tools.NewtonResult by @adtzlr in https://github.com/adtzlr/felupe/pull/603
Fix tools.moment(field, forces, boundary, centerpoint=np.zeros(3)) by @adtzlr in https://github.com/adtzlr/felupe/pull/605
Add tools.NewtonResult(xnorms=None, fnorms=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/606
Rename tools.save()-argument r to forces by @adtzlr in https://github.com/adtzlr/felupe/pull/607
Add new methods to Mesh by @adtzlr in https://github.com/adtzlr/felupe/pull/608
Rename arguments of new Mesh-methods by @adtzlr in https://github.com/adtzlr/felupe/pull/609
Remove fig-argument in CharacteristicCurve.plot(ax=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/610

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.13.0...v7.14.0

- Python
Published by adtzlr over 2 years ago

[7.13.0] - 2023-12-22

Added

Add NeoHookeCompressible for compressible hyperelastic materials or even as a fast alternative for NeoHooke when used in SolidBodyIncompressible.

Changed

Vectorize mesh.expand() which enhances the performance of quad/hex mesh-generators like Rectangle() and Cube().

Fixed

Fix logarithmic strain tensor evaluation in Job.evaluate(filename="result.xdmf") and in field.plot("Logarithmic Strain", component=0).

What's Changed

Vectorize mesh.expand() by @adtzlr in https://github.com/adtzlr/felupe/pull/584
Add compressible Neo-Hookean material formulation by @adtzlr in https://github.com/adtzlr/felupe/pull/586
Fix logarithmic strain tensor evaluation by @adtzlr in https://github.com/adtzlr/felupe/pull/589

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.12.0...v7.13.0

- Python
Published by adtzlr over 2 years ago

[7.12.0] - 2023-12-05

Added

Add plot- and screenshot-methods to Region and Scheme (base class for quadratures).
Add item = FormItem(bilinearform, linearform=None) to be used as an item in a Step(items=[item]).
Add a new method Boundary.apply_mask(mask). This simplifies re-definitions of boundary conditions with a custom mask.
Add support for two-dimensional dof-based masks in Boundary(mask) with mask.shape of (mesh.npoints, field.dim) in addition to point-based masks with mask.size of mesh.npoints.
Add a bubble-multiplier argument for RegionTriangleMINI(mesh, bubble_multiplier=0.1) and RegionTetraMINI(mesh, bubble_multiplier=0.1).
Add region.reload(mesh, element, quadrature) to re-evaluate a region, already linked to a field, with a modified mesh or element class or quadrature.

Changed

Refactor the assembly-submodule. Move the weak-form expression-related classes to the assembly.expression submodule.
Move Basis to the new assembly.expression submodule.
Make the field-submodule public.
Always import felupe as fem in docs and tests.
Change default optional (keyword) arguments of a weak-form expression decorator from Form(args=(), kwargs={}) to Form(args=None, kwargs=None).
Change default value of the skip-argument Boundary(skip=None). This will be set to (False, False, False) during initialization if mask=None.
Change the default bubble-multiplier in RegionTriangleMINI and RegionTetraMINI from 1.0 to 0.1. This affects only the template regions and not the element formulations TriangleMINI and TetraMINI, which still default to a bubble-multiplier of 1.0.
Pass optional keyword-arguments to math.einsum(**kwargs). This enables support for the out-argument.
Don't broadcast math.identity().
Rename quadrature/_base.py to quadrature/_scheme.py.

Fixed

Fix Boundary and subsequently also dof.symmetry() for different dimensions of the mesh and the field.
Fix negative cell-volumes error in RegionTriangleMINI for meshes like Rectangle(n=11).triangulate().add_midpoints_faces() by scaling down the (arbitrary) bubble-multiplier from 1.0 to 0.1.

Removed

Don't import Basis to the global namespace (not necessary as it is used only internally by the weak-Form expression decorator).
Remove unused internal assemble-methods from assembly.expression._linear.LinearForm and assembly.expression._bilinear.BilinearForm.
Remove extra-index-url https://wheels.vtk.org as they are now available on PyPI for Python 3.12.

What's Changed

Add Region.plot_scheme() by @adtzlr in https://github.com/adtzlr/felupe/pull/554
Add Region.plot() and quadrature.Scheme.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/556
Refactor the assembly submodule by @adtzlr in https://github.com/adtzlr/felupe/pull/558
Change private _field to public field submodule by @adtzlr in https://github.com/adtzlr/felupe/pull/560
Docs and Tests: Always import felupe as fem by @adtzlr in https://github.com/adtzlr/felupe/pull/565
Add FormItem(bilinearform, linearform=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/566
Add Boundary.apply_mask() by @adtzlr in https://github.com/adtzlr/felupe/pull/571
Boundary: Add support for a dof-based mask by @adtzlr in https://github.com/adtzlr/felupe/pull/573
VTK-wheels: Remove extra-index-url https://wheels.vtk.org by @adtzlr in https://github.com/adtzlr/felupe/pull/574
Add argument and change default RegionTriangleMINI(bubble_multiplier=0.1), also for RegionTetraMINI by @adtzlr in https://github.com/adtzlr/felupe/pull/576
Don't broadcast math.identity() and add optional kwargs inmath.einsum(**kwargs)` by @adtzlr in https://github.com/adtzlr/felupe/pull/579
Rename quadrature/_base.py to quadrature/_scheme.py. by @adtzlr in https://github.com/adtzlr/felupe/pull/580
Add Region.reload(mesh=None, element=None, quadrature=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/582

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.11.0...v7.12.0

- Python
Published by adtzlr over 2 years ago

[7.11.0] - 2023-10-22

Added

Add cell-type argument to Mesh.add_midpoints_volumes(cell_type=None) and its variants for edges and faces.
Add element.Element.view(), element.Element.plot() and element.Element.screenshot(). This enables an interactive plot of the element in the reference configuration with its point ids, e.g. Hexahedron().plot().show().

Changed

Change function signature and enhance dof.biaxial(field, lefts=(None, None), rights=(None, None), moves=(0.2, 0.2), axes=(0, 1), clampes=(False, False), sym=True). Now with a full-featured docstring including an example. The old-style input move defaults to None, is deprecated and will be removed in the next major release.
Change function signature and enhance dof.shear(field, bottom=None, top=None, moves=(0.2, 0.0, 0.0), axes=(0, 1), sym=True). Now with a full-featured docstring including an example. The old-style move, axis_compression, axis_shear and compression all default to None, are deprecated and will be removed in the next major relase.
Merge keyword-arguments for the dual-regions with hard-coded arguments in FieldsMixed(region, **kwargs).
Replace np.product() (will be removed in NumPy 2.0) with the equivalent np.prod().

Fixed

Fix FieldsMixed() for regions with MINI-element formulations: Disable the disconnection of the dual mesh.
Fix dof.shear(sym=True) which was previously ignored due to a wrong setup of the symmetry boundaries.
Fix the install command on Python 3.12 by adding an extra-index-url for VTK wheels if they are not yet available on PyPI (the extra index is provided by Kitware).
Fix a warning because the timings of the Newton-Rhapson solver are printed from a one-dimensional array. Take the first item of the runtime-array to resolve this warning.

Removed

Remove the undocumented dof.planar() because this is a special case of the biaxial load case dof.biaxial(field, clampes=(True, False), moves=(0.2, 0), sym=False, axes=(0, 1)).

What's Changed

Enhance loadcase biaxial by @adtzlr in https://github.com/adtzlr/felupe/pull/534
Remove planar load case by @adtzlr in https://github.com/adtzlr/felupe/pull/535
Fix FieldsMixed() for MINI-elements by @adtzlr in https://github.com/adtzlr/felupe/pull/538
Fix dof.shear(sym=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/540
FieldsMixed(region, **kwargs): Merge kwargs for the dual region by @adtzlr in https://github.com/adtzlr/felupe/pull/541
Add support for Python 3.12 by @adtzlr in https://github.com/adtzlr/felupe/pull/543
Add element methods view, plot and screenshot by @adtzlr in https://github.com/adtzlr/felupe/pull/551
Deprecate old-style args of load cases by @adtzlr in https://github.com/adtzlr/felupe/pull/553

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.10.0...v7.11.0

- Python
Published by adtzlr over 2 years ago

[7.10.0] - 2023-09-14

Added

Add ViewField as method to a field container FieldContainer.view(). Now a field container provides the view sub-methods directly, i.e. add FieldContainer.plot() and FieldContainer.screenshot(filename="mesh.png").

Changed

Hide the undeformed mesh in Mesh.plot() by default (also affects Mesh.screenshot()).

Fixed

Fix taking screenshots of a mesh.

What's Changed

Add ViewField(field) as method to FieldContainer.view() by @adtzlr in https://github.com/adtzlr/felupe/pull/533

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.9.0...v7.10.0

- Python
Published by adtzlr over 2 years ago

[7.9.0] - 2023-09-13

Added

Add ViewMesh as method to a mesh Mesh.view(). Now a mesh provides the view sub-methods directly, i.e. add Mesh.plot() and Mesh.screenshot(filename="mesh.png").

What's Changed

Add plot-related methods from ViewMesh to Mesh: view, plot, screenshot by @adtzlr in https://github.com/adtzlr/felupe/pull/532

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.8.0...v7.9.0

- Python
Published by adtzlr over 2 years ago

[7.8.0] - 2023-09-05

Added

Add mesh.Triangle(a, b, c, n=2) for the creation of a quad-meshed triangle.

Changed

Change mesh.Circle(n=6) to the minimum mesh.Circle(n=2).
Enhance dof.uniaxial(axis=0, sym=(False, True, False)) by a user-defined axis and selective symmetries.
Enhance dof.shear(axis_compression=1, axis_shear=0, sym=True) by user-defined axes of compression and shear.

What's Changed

Add mesh.Triangle(a, b, c, n) by @adtzlr in https://github.com/adtzlr/felupe/pull/526
Fix mesh.concatenate() for a list of length one by @adtzlr in https://github.com/adtzlr/felupe/pull/528
Enhance uniaxial- and shear-loadcases by @adtzlr in https://github.com/adtzlr/felupe/pull/529

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.7.0...v7.8.0

- Python
Published by adtzlr over 2 years ago

[7.7.0] - 2023-08-31

Added

Add optional normalization of mesh runouts (which are then indents) by mesh.runouts(normalize=False).
Add LinearElasticLargeStrain(E=None, nu=None, parallel=False), suitable for large-rotation analyses. This is based on NeoHooke() with converted Lamé-constants.
Add a simple boundary-based quad- or hex mesher: A mesh tool for filling the face or volume between two line- or quad-meshes mesh.fill_between(mesh, other_mesh, n=11).
Add Circle(radius, centerpoint, n) for the creation of a quad-meshed circle.

Changed

Update the mesh also with a new points array: this changes the Mesh-update function mesh.update(cells, cell_type=None) to mesh.update(points=None, cells=None, cell_type=None). Note that this could break old scripts which use mesh.update(new_cells) instead of mesh.update(cells=new_cells).
Move the copy-method of a Mesh Mesh.copy() to its base-class and extend it to optionally update the attributes DiscreteGeometry.copy(points=None, cells=None, cell_type=None).

Removed

Remove tests on Python 3.7 (end of life).

What's Changed

Normalize runouts of meshes mesh.runouts(normalize=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/510
Add LinearElasticLargeStrain for linear-elastic materials by @adtzlr in https://github.com/adtzlr/felupe/pull/514
Add a minimalistic boundary-based quad- or hex-mesh tool by @adtzlr in https://github.com/adtzlr/felupe/pull/519
Update-method of Mesh: Add optional points-argument by @adtzlr in https://github.com/adtzlr/felupe/pull/521
Move Mesh.copy() to DiscreteGeometry.copy() by @adtzlr in https://github.com/adtzlr/felupe/pull/523
Add Circle(radius, centerpoint, n) for creating quad-meshed circles by @adtzlr in https://github.com/adtzlr/felupe/pull/524

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.6.1...v7.7.0

- Python
Published by adtzlr almost 3 years ago

What's Changed

Start virtual framebuffer for notebook plotting of views by @adtzlr in https://github.com/adtzlr/felupe/pull/508

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.6.0...v7.6.1

- Python
Published by adtzlr almost 3 years ago

What's Changed

Add View.plot(notebook=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/507

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.5.1...v7.6.0

- Python
Published by adtzlr almost 3 years ago

What's Changed

Don't set a pyvista-theme by default in View by @adtzlr in https://github.com/adtzlr/felupe/pull/505

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.5.0...v7.5.1

- Python
Published by adtzlr almost 3 years ago

[7.5.0] - 2023-07-20

Added

Add ViewSolid, which enables view = ViewSolid(field, solid=None) the view of cauchy stresses, e.g. view.plot("Principal Values of Cauchy Stress").show().
Add constitutive models to top-level namespace, e.g. yeoh() from constitution.yeoh(). This makes typing hyperelastic material formulations shorter: Hyperelastic(yeoh, C10=0.5, C20=-0.1, C30=0.02).
Add CharacteristicCurve.plot(swapaxes=False).
Add MaterialAD: A user-defined material definition with Automatic Differentiation. Only the first Piola-Kirchhoff stress tensor must be provided.

Changed

Add optional point- and cell-data args for ViewMesh(mesh, point_data=None, cell_data=None) like already implemented in ViewField.
Enforce contiguous arrays in UserMaterialHyperelastic (enhance performance).
View: Switch from ViewField to ViewSolid.
View: Always plot the undeformed mesh with opacity=0.2 and show_edges=False.
Rename UserMaterial to Material, UserMaterialStrain to MaterialStrain, UserMaterialHyperelastic to Hyperelastic (keep old alias names until next major release).
Use consistent indices in einsum() for (elementwise operating) trailing axes: q for quadrature point and c for cell.
Rename internal IntegralFormMixed to IntegralForm, which is now consistent internally and in the top-level namespace. The previous internal base-class for a single-field IntegralForm is renamed to WeakForm.
Don't plot x- and y-labels in CharacteristicCurve.plot(xlabel=None, ylabel=None) if they are not specified.

Fixed

Don't warp the mesh in ViewMesh.plot().
Warp the mesh in case no name is passed in View.plot(name=None).
Don't modify a given label in Scene.plot(label=None).
Fix the second invariant of the distortional part of the right Cauchy-Green deformation tensor in hyperelastic material formulations using tensortrax, i.e. fix the implementations of mooney_rivlin(), third_order_deformation() and van_der_waals().

Removed

Remove internal (unused) imports of the assembly submodule. ## What's Changed
Enhance View: Add optional point- and cell-data to ViewMesh by @adtzlr in https://github.com/adtzlr/felupe/pull/475
Add ViewSolid(field, solid=None) and set View = ViewSolid by @adtzlr in https://github.com/adtzlr/felupe/pull/479
Warp the mesh for View.plot(name=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/481
Rename UserMaterial... to shorter names by @adtzlr in https://github.com/adtzlr/felupe/pull/483
Add constitutive material formulations to top-level namespace by @adtzlr in https://github.com/adtzlr/felupe/pull/485
Docs: Add Example Inflation by @adtzlr in https://github.com/adtzlr/felupe/pull/484
Consistent idx in einsum(), Enhance FieldsMixed(n=1) by @adtzlr in https://github.com/adtzlr/felupe/pull/489
Basis: Rename idx in einsum() by @adtzlr in https://github.com/adtzlr/felupe/pull/490
Rename internal IntegralFormMixed to IntegralForm by @adtzlr in https://github.com/adtzlr/felupe/pull/491
Remove unused imports of the assembly submodule by @adtzlr in https://github.com/adtzlr/felupe/pull/492
Don't modify a given label in Scene.plot(label=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/495
Enhance plotting in a characteristic-curve job by @adtzlr in https://github.com/adtzlr/felupe/pull/498
Fix mooney_rivlin(), third_order_deformation() and van_der_waals() by @adtzlr in https://github.com/adtzlr/felupe/pull/502
Add MaterialAD: PK1-based Material w/ Automatic Differentiation by @adtzlr in https://github.com/adtzlr/felupe/pull/493

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.4.1...v7.5.0

- Python
Published by adtzlr almost 3 years ago

[7.4.1] - 2023-05-02

Changed

Change the logo.

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.4.0...v7.4.1

- Python
Published by adtzlr about 3 years ago

[7.4.0] - 2023-04-29

Added

Add ViewMesh(mesh) and rebase View (internally renamed to ViewField) on ViewMesh with additional point- and cell-data.

Changed

Change math.linsteps(axis=None, axes=None) to create optional multi-column arrays, where the steps are inserted at the given axis.
Make last gravity and density arguments of SolidBodyGravity(field, gravity=None, density=1.0) optional.

What's Changed

Add ViewMesh(mesh, cell_type=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/470
Multi-column math.linsteps(axis=None, axes=None) and optional args in SolidBodyGravity(field, gravity=None, density=1.0) by @adtzlr in https://github.com/adtzlr/felupe/pull/472

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.3.0...v7.4.0

- Python
Published by adtzlr about 3 years ago

[7.3.0] - 2023-04-28

Changed

Show a progress bar during Job.evaluate(verbose=True) (new optional dependency tqdm). The legacy detailed output is available with Job.evaluate(verbose=2).

Removed

Remove config files for MyBinder. They are now located in a different repository adtzlr/felupe-web.

What's Changed

Job: Add progress bar by @adtzlr in https://github.com/adtzlr/felupe/pull/467
Job.evaluate(verbose=True): Further reduce the output by @adtzlr in https://github.com/adtzlr/felupe/pull/469

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.2.0...v7.3.0

- Python
Published by adtzlr about 3 years ago

[7.2.0] - 2023-04-26

Added

Add environment.yml config file for MyBinder.
Add a timetrack-list as Job.timetrack which is updated incrementally on Job.evaluate().
Add View(field, point_data=None, cell_data=None), a result plotter powered by pyvista.
Add ViewXdmf(filename, time=0), a result plotter powered by pyvista.

Changed

Make everything in /src compliant with flake8.
Generalize the math-module so that all functions handle an arbitrary number of elementwise-operating trailing axes.
The special contraction modes of math.dot(mode=(2,2)) and math.ddot(mode=(2,2)) have to be specified by the mode-argument and are not detected by the shapes of the operands.
Enhance the overall performance by enforcing the identity matrix to a C-contiguous array.
Change point- and cell-data functions used in Job.evaluate(point_data=None, cell_data=None) from fun(substep) to fun(field, substep).

Fixed

Fix timings shown in newtonrhapson(verbose=True): The solve time was only related to one call of the solver while the assembly time referred to the whole runtime subtracted by the single-call solve time.

What's Changed

Flake8 Compliance (without tests) by @adtzlr in https://github.com/adtzlr/felupe/pull/448
Flake8: Remove E741 from ignore list by @adtzlr in https://github.com/adtzlr/felupe/pull/449
math: Handle arbitrary number of trailing axes by @adtzlr in https://github.com/adtzlr/felupe/pull/451
Fix timings in newtonrhapson() by @adtzlr in https://github.com/adtzlr/felupe/pull/453
Job: Add timetrack attribute by @adtzlr in https://github.com/adtzlr/felupe/pull/455
Add Result() for plotting XDMF-files by @adtzlr in https://github.com/adtzlr/felupe/pull/456
Enhance and simplify Result by @adtzlr in https://github.com/adtzlr/felupe/pull/458
Result: Switch to more default args of pyvista by @adtzlr in https://github.com/adtzlr/felupe/pull/460
Job.evaluate(): Change point- and cell-data function arguments by @adtzlr in https://github.com/adtzlr/felupe/pull/461
Result.plot(name=None): Plot only the undeformed mesh by @adtzlr in https://github.com/adtzlr/felupe/pull/462
Rename Result to View by @adtzlr in https://github.com/adtzlr/felupe/pull/463
View: Add cell_type-argument by @adtzlr in https://github.com/adtzlr/felupe/pull/464
Docs: Add Notebook support by @adtzlr in https://github.com/adtzlr/felupe/pull/465

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.1.0...v7.2.0

- Python
Published by adtzlr about 3 years ago

[7.1.0] - 2023-04-15

Added

Add string representations for Region and FieldContainer.
Add Job.evaluate(parallel=True) in addition to Job.evaluate(kwargs={"parallel": True}). If both are given, the key in the dict is overwritten by the user-defined value.
Add mesh.stack(meshes) for joining meshes with identical points-arrays and cell-types. Contrary to mesh.concatenate(meshes), the points are not stacked and no offsets are inserted into the cells-arrays.
Add mesh.translate(move, axis) for the translation of mesh-points.

Changed

Pass optional keyword-arguments in CharacteristicCurve.plot(**kwargs) to the figure.
Don't invoke CharacteristicCurve.evaluate() from CharacteristicCurve.plot(), raise an error if the current job is not evaluated instead.
Make the endpoint of math.linsteps(endpoint=True) optional.
Don't modify the mesh for the dual regions RegionConstantQuad() and RegionConstantHexahedron(). Instead, it is required to pass a dual (disconnected) mesh with one point per cell RegionConstantQuad(mesh.dual(points_per_cell=1)).
Make requirement einsumt optional again due to issues with JupyterLite.
Add matplotlib to optional requirements.

Fixed

Catch ModuleNotFoundError if from einsumt import einsumt fails (in JupyterLite) and fall back to from numpy import einsum as einsumt.

What's Changed

Add string representations for Region and FieldContainer by @adtzlr in https://github.com/adtzlr/felupe/pull/432
CharacteristicCurve: Pass kwargs to the plot by @adtzlr in https://github.com/adtzlr/felupe/pull/435
Add Job.evaluate(parallel=False) by @adtzlr in https://github.com/adtzlr/felupe/pull/437
Add Mesh.dual() and rebase Mesh.disconnect(), create dual-meshes in FieldsMixed() by @adtzlr in https://github.com/adtzlr/felupe/pull/441
Add mesh.stack(meshes) and mesh.translate(move, axis) by @adtzlr in https://github.com/adtzlr/felupe/pull/443
Threaded einsum: Fall-back to numpy if einsumt is not available by @adtzlr in https://github.com/adtzlr/felupe/pull/446

Full Changelog: https://github.com/adtzlr/felupe/compare/v7.0.0...v7.1.0

- Python
Published by adtzlr about 3 years ago

[7.0.0] - 2023-04-07

Added

Add boundary regions RegionQuadraticQuadBoundary and RegionBiQuadraticQuadBoundary for quadratic quads.
Add boundary regions RegionQuadraticHexahedronBoundary and RegionTriQuadraticHexahedronBoundary for quadratic hexahedrons.
Add mesh.flip(mask=None) to flip a mirrored or wrong indexed cells array, applied on a given boolean mask of cells.

Changed

Change einsumt from an optional to a required dependency.
Vectorize implementations of MultiPointConstraint and MultiPointContact and re-implement both as scipy.sparse.lil_matrix().
Rename old Mesh to DiscreteGeometry and rebase new Mesh on DiscreteGeometry.
Simplify the usage of explicit mesh-related tools by adding them as methods to Mesh, i.e. mesh.tools.expand(Rectangle()) is equivalent to Rectangle().expand().
Print runtimes for time spent on Assembly and Solve in newtonrhapson(verbose=True).
Check for negative cell-volumes at quadrature points and print a warning along with a possible fix.

Fixed

Fix tools.project() for higher-order quad- and hexahedron elements.
Fix transposed output of tools.project().
Fix failed scalar-value checks by using np.isscalar() in mesh.expand(z=1) and mesh.revolve(phi=180) where z or phi are of type np.int32.
Fix read a Mesh with no cells in mesh.read().

Removed

Remove jit-compilation of forms (parallel is the preferred method).
Remove unused tools.check().
Remove optional dependency sparse.

What's Changed

Remove jit and reduce complexity of newtonrhapson() by @adtzlr in https://github.com/adtzlr/felupe/pull/408
Enforce einsumt by @adtzlr in https://github.com/adtzlr/felupe/pull/409
Vectorize implementation of MultiPointConstraint by @adtzlr in https://github.com/adtzlr/felupe/pull/412
Add boundary regions for quadratic quads by @adtzlr in https://github.com/adtzlr/felupe/pull/414
Add boundary regions for quadratic hexahedrons by @adtzlr in https://github.com/adtzlr/felupe/pull/415
Fix TriQuadraticHexahedron().points and tools.project by @adtzlr in https://github.com/adtzlr/felupe/pull/419
Fix scalar-checks in mesh.expand() and mesh.revolve() by @adtzlr in https://github.com/adtzlr/felupe/pull/421
Fix read a mesh with no cells in mesh.read() by @adtzlr in https://github.com/adtzlr/felupe/pull/423
Enhance Mesh: Mesh tools as methods by @adtzlr in https://github.com/adtzlr/felupe/pull/425
Print newton-timings for time spent on Assembly and Solve by @adtzlr in https://github.com/adtzlr/felupe/pull/427
Fix MultiPointConstraint and MultiPointContact with active skips by @adtzlr in https://github.com/adtzlr/felupe/pull/429

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.4.0...v7.0.0

- Python
Published by adtzlr about 3 years ago

[6.4.0] - 2023-04-01

Added

Add a new argument to pass a mesh for the dual regions in FieldsMixed(mesh=None).
Add quadrature and grad arguments to RegionLagrange.
Add order attribute to RegionLagrange.
Add items-argument for custom slicing of characteristic curve plots in CharacteristicCurve.plot(items=None).

Changed

Enhance Domain integration in IntegralForm: Ensure C-contiguous arrays as fun-argument.
Enhance performance of hyperelastic material model formulations using automatic differentiation.
Perform reshape and broadcasting if user-defined Field-values with correct size are given.
Make symmetry-related boundary conditions in uniaxial loadcase optional dof.uniaxial(sym=True).
Merge custom point- and cell-data dicts with default dicts in Job instead of replacing them, also add Job(point_data_default=True, cell_data_default=True).
Allow to change cell-type in Mesh.update(cells, cell_type=None).
Enhance the creation of a disconnected mesh for mixed-field formulations by Mesh.disconnect(points_per_cell=None, calc_points=True).
Change required to optional arguments in Step(items, ramp=None, boundaries=None).

Fixed

Fix broadcast arrays for the geometric stiffness contribution of hyperelastic material model formulations using automatic differentiation.

What's Changed

IntegralForm: Ensure C-contiguous arrays as fun-argument by @adtzlr in https://github.com/adtzlr/felupe/pull/386
Perform reshape and broadcasting if user-defined Field-values with correct size are given. by @adtzlr in https://github.com/adtzlr/felupe/pull/388
Make symmetry boundary conditions optional in dof.uniaxial(sym=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/392
Job: Append point- and cell-data to default instead of replacing them by @adtzlr in https://github.com/adtzlr/felupe/pull/394
Add a custom mesh for the dual regions in FieldsMixed(mesh=None) by @adtzlr in https://github.com/adtzlr/felupe/pull/396
Enhance FieldsMixed(): Add support for RegionLagrange (family of Taylor-Hood elements) by @adtzlr in https://github.com/adtzlr/felupe/pull/400
Characteristic Curve - Job plot: Add items-argument for custom slicing by @adtzlr in https://github.com/adtzlr/felupe/pull/390
Enhance Mesh.disconnect(points_per_cell=None, calc_points=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/403
Step: Make arguments ramp and boundaries optional by @adtzlr in https://github.com/adtzlr/felupe/pull/405

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.3.0...v6.4.0

- Python
Published by adtzlr about 3 years ago

[6.3.0] - 2023-02-06

Added

Add more hyperelastic material formulations: Saint Venant-Kirchhoff, Mooney-Rivlin, Arruda-Boyce, Extended Tube and Van der Waals.
Add BiQuadraticQuad element (as a permuted version of ArbitraryOrderLagrangeElement).
Add Quadratic element.
Add new arguments for more flexible region templates, e.g. RegionQuad(mesh, quadrature=GaussLegendre(order=1, dim=2), grad=True).
Add support for triangle/tetra in FieldsMixed().
Add optional state variables to UserMaterialHyperelastic(model, nstatevars=0).
Add finite-strain viscoelastic material formulation.

Changed

Switch to a src/-layout.
Import base-element class as element.Element().
Auto-detect min/max-coordinates of mesh-points for loadcases if dof.uniaxial(right=None).
Ensure compatibility with tensortrax>0.6.0 (Van-der-Waals and Viscoelastic models).

Fixed

Fix rotation matrix for a rotation around the y-axis.

What's Changed

Add more hyperelastic material formulations by @adtzlr in https://github.com/adtzlr/felupe/pull/369
Fix rotation matrix around y-axis by @adtzlr in https://github.com/adtzlr/felupe/pull/370
Add QuadraticQuad() and BiQuadraticQuad() elements by @adtzlr in https://github.com/adtzlr/felupe/pull/371
Add support for dual triangle/tetra-regions by @adtzlr in https://github.com/adtzlr/felupe/pull/373
Switch to src/-layout by @adtzlr in https://github.com/adtzlr/felupe/pull/374
Add optional state variables to UserMaterialHyperelastic(nstatevars=0) by @adtzlr in https://github.com/adtzlr/felupe/pull/377
Import base-element class element.Element() by @adtzlr in https://github.com/adtzlr/felupe/pull/378
Change loadcases: Detect Min/Max Coordinates of mesh.points by @adtzlr in https://github.com/adtzlr/felupe/pull/381
Fix tensortrax>0.6.0 compatibility by @adtzlr in https://github.com/adtzlr/felupe/pull/382

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.2.5...v6.3.0

- Python
Published by adtzlr over 3 years ago

[6.2.5] - 2023-01-02

Fixed

Once again fix init if tensortrax is not installed.

What's Changed

Fix import of UserMaterialHyperelastic by @adtzlr in https://github.com/adtzlr/felupe/pull/367

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.2.4...v6.2.5

- Python
Published by adtzlr over 3 years ago

[6.2.4] - 2023-01-01

Changed

Update actions used in CI/CD workflows.

- Python
Published by adtzlr over 3 years ago

[6.2.3] - 2023-01-01

Changed

Remove setup.cfg, change pyproject.toml and store the version tag only once within the source code (__about__.py).

What's Changed

Change pyproject.toml by @adtzlr in https://github.com/adtzlr/felupe/pull/362

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.2.2...v6.2.3

- Python
Published by adtzlr over 3 years ago

[6.2.2] - 2022-12-20

Fixed

Fix init if tensortrax is not installed.

What's Changed

Fix init if tensortrax is not installed. by @adtzlr in https://github.com/adtzlr/felupe/pull/360

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.2.1...v6.2.2

- Python
Published by adtzlr over 3 years ago

[6.2.1] - 2022-12-19

Fixed

Fix version string.

What's Changed

fix version string by @adtzlr in https://github.com/adtzlr/felupe/pull/359

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.2.0...v6.2.1

- Python
Published by adtzlr over 3 years ago

[6.2.0] - 2022-12-16

Added

Add Total-Lagrange UserMaterialHyperelastic(fun, parallel=False, **kwargs) based on optional tensortrax. Only available if tensortrax is installed.
Add constitutive isotropic hyperelastic model formulations to be used with UserMaterialHyperelastic() (constitution.ogden(C, mu, alpha), etc.).

Auto-Generated Release-Notes

What's Changed

Add UserMaterialHyperelastic() by @adtzlr in https://github.com/adtzlr/felupe/pull/355
Make tensortrax optional by @adtzlr in https://github.com/adtzlr/felupe/pull/357
Add models for isotropic hyperelasticity UserMaterialHyperelastic() by @adtzlr in https://github.com/adtzlr/felupe/pull/358

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.1.0...v6.2.0

- Python
Published by adtzlr over 3 years ago

[6.1.0] - 2022-12-10

Changed

Enhance plotting with custom x- and y-data in CharacteristicCurve.plot(x, y) and allow a list of items for force evaluation in CharacteristicCurve(items=[...]) to be passed.
Enhance math.linsteps(points=[0, 5, 0], num=[5, 10]) by supporting a list of substeps.
Enhance compression of shear loadcase: Apply the compression on the bottom and the top dof.shear(compression=(0, 0)).

Auto-Generated Release-Notes

What's Changed

Enhance CharacteristicCurve(items=[...]).plot(x, y) by @adtzlr in https://github.com/adtzlr/felupe/pull/350
Enhance linsteps by supporting a list of substeps by @adtzlr in https://github.com/adtzlr/felupe/pull/352
Enhance dof.shear(compression=(0, 0) by @adtzlr in https://github.com/adtzlr/felupe/pull/354

Full Changelog: https://github.com/adtzlr/felupe/compare/v6.0.0...v6.1.0

- Python
Published by adtzlr over 3 years ago

[6.0.0] - 2022-11-20

Added

Add project(mean=True) to project cell mean-values to mesh-points. Now project() supports Triangles and Tetrahedrons.
Add RegionBoundary.mesh_faces() for a mesh with face-cells on the selected boundary of a region.
Add pseudo-elastic material OgdenRoxburgh() which may be used with SolidBodyNearlyIncompressible().
Add umat = UserMaterial(stress, elasticity, nstatevars=0, **kwargs) with user-defined functions for the (first Piola-Kirchhoff) stress tensor P, statevars_new = umat.gradient([F, statevars], **kwargs) and the according fourth-order elasticity tensor A = umat.hessian([F, statevars], **kwargs) based on the deformation gradient.
Add UserMaterialStrain() for small-strain based user-defined material formulations with an umat-interface suitable for elastic-plastic frameworks.
Add LinearElasticPlasticIsotropicHardening() which is based on UserMaterialStrain() and constitution.linear_elastic_plastic_isotropic_hardening().
Add new math helpers math.ravel() and math.reshape().
Add optional axis argument on which the norm is evaluated math.norm(axis=None).

Changed

Unify material definition with methods for the stress P, statevars_new = umat.gradient([F, statevars]) and the elasticity tensor A = umat.hessian([F, statevars]). This breaks support for materials defined by matadi<=0.1.10.
Do not broadcast the (constant) elasticity tensor for linear-elastic materials as einsumt>=0.9.3 supports broadcasting along the parallel-executed dimension.
Change not-updating attribute of FieldContainer(fields).values to a method FieldContainer(fields).values() which returns the current field values.

Removed

Remove unused SolidBodyTensor() and SolidBodyTensorNearlyIncompressible().
Remove unused region argument of LinearElastic().hessian().

Auto-Generated Release-Notes

What's Changed

Add project(mean=True) by @adtzlr in https://github.com/adtzlr/felupe/pull/329
Add RegionBoundary.mesh_faces() by @adtzlr in https://github.com/adtzlr/felupe/pull/334
Add pseudo-elastic material OgdenRoxburgh() by @adtzlr in https://github.com/adtzlr/felupe/pull/335
Unify umat for hyperelastic and tensor-based material formulations by @adtzlr in https://github.com/adtzlr/felupe/pull/339
Do not broadcast linear-elastic elasticity tensor by @adtzlr in https://github.com/adtzlr/felupe/pull/342
FieldContainer: Change attribute values to method values() by @adtzlr in https://github.com/adtzlr/felupe/pull/344
Add UserMaterialStrain() by @adtzlr in https://github.com/adtzlr/felupe/pull/346
Enhance and Fix the Documentation by @adtzlr in https://github.com/adtzlr/felupe/pull/347

Full Changelog: https://github.com/adtzlr/felupe/compare/v5.3.1...v6.0.0

- Python
Published by adtzlr over 3 years ago

Changelog since [5.3.0] - 2022-11-03

Fixed

Fix volume evaluation of (nearly) incompressible solids for axisymmetric fields.

Auto-Generated Release-Notes

What's Changed

Fix incompressible solid: Wrong volume evaluation for axisymmetric fields by @adtzlr in https://github.com/adtzlr/felupe/pull/327

Full Changelog: https://github.com/adtzlr/felupe/compare/v5.3.0...v5.3.1

- Python
Published by adtzlr over 3 years ago

Changelog since [5.2.0] - 2022-10-08

Added

Add optional pre-compression to shear-loadcase dof.shear(compression=0.0).
Add MeshContainer and string-representation for Mesh objects.
Add a mesh-reader using meshio mesh.read(filename, ...).
Add SolidBodyNearlyIncompressible(umat, field, bulk) for (nearly) incompressible solids and a given (distortional-part only) constitutive material formulation. This is a pure displacement-based alternative to the three-field-formulation technique.

Changed

Support an optional user-defined meshio-object in Job().evaluate(mesh=None, filename="result.xdmf").
Support a distortional-part only Neo-Hookean material formulation with no bulk modulus defined NeoHooke(mu=1.0).

Fixed

Fix missing ArbitraryOrderLagrangeElement.points attribute.
Fix ignored mask only_surface=True for RegionBoundary().mesh.cells_faces.
Set default pressure to zero in SolidBodyPressure().
Take the mesh from the global x0-field if x0 is passed to job.evaluate(x0=...).
Fix missing update of global field x0 in job.evaluate(x0) after each completed substep.

Auto-Generated Release-Notes

What's Changed

Add points attribute to lagrange element by @adtzlr in https://github.com/adtzlr/felupe/pull/301
Add pre-compression to shear loadcase by @adtzlr in https://github.com/adtzlr/felupe/pull/305
Fix ignored surface mask for cells_faces of RegionBoundary by @adtzlr in https://github.com/adtzlr/felupe/pull/307
Add MeshContainer and string-representation for Mesh objects by @adtzlr in https://github.com/adtzlr/felupe/pull/312
Enhance Mesh-Container by @adtzlr in https://github.com/adtzlr/felupe/pull/313
Add mesh.read(filename, ...) by @adtzlr in https://github.com/adtzlr/felupe/pull/316
job.evaluate(x0=...): Take mesh from global field x0 if present by @adtzlr in https://github.com/adtzlr/felupe/pull/320
Add SolidBodyNearlyIncompressible(umat, field, bulk) by @adtzlr in https://github.com/adtzlr/felupe/pull/323
Fix missing update of x0 in job.evaluate(x0) after each completed substep by @adtzlr in https://github.com/adtzlr/felupe/pull/324

Full Changelog: https://github.com/adtzlr/felupe/compare/v5.2.0...v5.3.0

- Python
Published by adtzlr over 3 years ago

Changelog since [5.1.0] - 2022-09-09

Added

Add xscale and yscale arguments to CharacteristicCurve.plot().
Add mesh.Grid(*xi) as generalized line, rectangle or cube with custom linspaces.
Add mesh.concatenate(meshes) to join a sequence of meshes with identical cell types.
Add x0 argument to Job.evaluate(x0=field).
Add mask argument to mesh.runouts(mask=slice(None)).
Add callback(stepnumber, substepnumber, substep) argument to CharacteristicCurve() (like in Job()).
Add an on-the-fly XDMF writer for a job (via meshio) Job.evaluate(filename="result.xdmf") with the possibility to add optional point_data and cell_data dicts.

Changed

Remove Warning if einsumt requirement is not found (switch to numpy without any warnings).
Requires Python 3.7+.

Fixed

Fix ignored axis argument of mesh.revolve(axis=1).

Auto-Generated Release-Notes

What's Changed

Add xscale and yscale arguments to CharacteristicCurve.plot() by @adtzlr in https://github.com/adtzlr/felupe/pull/281
Add mesh.Grid(*xi) and mesh.concatenate([mesh1, mesh2, ...]) by @adtzlr in https://github.com/adtzlr/felupe/pull/284
Add job.evaluate(x0=field) by @adtzlr in https://github.com/adtzlr/felupe/pull/286
Fix mesh.revolve(axis=1) by @adtzlr in https://github.com/adtzlr/felupe/pull/288
Add callback argument to CharacteristicCurve() by @adtzlr in https://github.com/adtzlr/felupe/pull/290
Remove Warning if einsumt is not installed by @adtzlr in https://github.com/adtzlr/felupe/pull/292
Add XDMF writer by @adtzlr in https://github.com/adtzlr/felupe/pull/294

Full Changelog: https://github.com/adtzlr/felupe/compare/v5.1.0...v5.2.0

- Python
Published by adtzlr over 3 years ago

Changelog since [5.0.0] - 2022-08-21

Changed

Enhance Boundary: Select Points by value in addition to a callable (fx=lambda x: x == 0 is equivalent to fx=0), also add mode="and" and mode="or" argument.
Support line elements within the revolution function mesh.revolve().
Import previously hidden functions fun_items() and jac_items() as tools.fun() and tools.jac(), respectively (useful for numeric continuation).
Add step- and substep-numbers as arguments to the callback(stepnumber, substepnumber, substep)-function of a Job.

Auto-Generated Release-Notes

What's Changed

Enhance Boundary class by @adtzlr in https://github.com/adtzlr/felupe/pull/274
Support Line-elements within mesh.revolve() by @adtzlr in https://github.com/adtzlr/felupe/pull/276
Import fun_items() and jac_items() from felupe.tools._newton by @adtzlr in https://github.com/adtzlr/felupe/pull/277
Docs: Add numeric continuation tutorial with contique by @adtzlr in https://github.com/adtzlr/felupe/pull/278
Callback: Add args (step- and substep-numbers) by @adtzlr in https://github.com/adtzlr/felupe/pull/279

Full Changelog: https://github.com/adtzlr/felupe/compare/v5.0.0...v5.1.0

- Python
Published by adtzlr over 3 years ago

Changelog since [4.0.0] - 2022-08-07

Added

Add SolidBodyGravity for body forces acting on a solid body.
Support list of linked fields in Newton-Rhapson solver newtonrhapson(fields=[field_1, field_2]).
Automatic init of state variables in SolidBodyTensor.
Add mesh.runouts() for the creation of runouts of rubber-blocks of rubber-metal structures.
Add FieldPlaneStrain which is a 2d-field and returns gradients of shape (3, 3) (for plane strain problems with 3d user materials).
Add PointLoad for the creation of external force vectors.
Add Step with a generator for substeps, Job and CharacteristicCurve.

Changed

Move MultiPointConstraint to mechanics module and unify handling with SolidBody.
Rename bodies argument of Newton-Rhapson solver to items (now supports MPC).
Return partitioned system as dict from loadcases loadcase=dict(dof0=dof0, dof1=dof1, ext0=ext0).
Check function residuals norm in newtonrhapson() instead of incremental field-values norm.

Fixed

Fix assembled vectors and results of SolidBodyPressure for initially defined pressure values.
Fix verbose=0 option of newtonrhapson().
Fix wrong assembly of axisymmetric mixed-fields due to introduced plane strain field-trimming.

Auto-Generated Release-Notes

What's Changed

Add SolidBodyGravity an Fix SolidBodyPressure by @adtzlr in https://github.com/adtzlr/felupe/pull/237
Modernize MultiPointConstraint by @adtzlr in https://github.com/adtzlr/felupe/pull/238
Improve loadcase by @adtzlr in https://github.com/adtzlr/felupe/pull/240
Newton-Rhapson: Support multiple fields by @adtzlr in https://github.com/adtzlr/felupe/pull/242
Newton remove fields arg by @adtzlr in https://github.com/adtzlr/felupe/pull/243
Fix multiple mixed subfields by @adtzlr in https://github.com/adtzlr/felupe/pull/245
Enhance SolidBodyTensor: Automatic init of statevars by @adtzlr in https://github.com/adtzlr/felupe/pull/247
Fix verbose of newton by @adtzlr in https://github.com/adtzlr/felupe/pull/249
Check function residuals norm in newton by @adtzlr in https://github.com/adtzlr/felupe/pull/251
Fix update of state variables in Newton-Rhapson by @adtzlr in https://github.com/adtzlr/felupe/pull/253
Add Mesh-tool mesh.runouts() by @adtzlr in https://github.com/adtzlr/felupe/pull/255
Add FieldPlaneStrain by @adtzlr in https://github.com/adtzlr/felupe/pull/257
Field Integration: Fix field-trimming for plane strain by @adtzlr in https://github.com/adtzlr/felupe/pull/260
Multiply force vector of gravity by -1 by @adtzlr in https://github.com/adtzlr/felupe/pull/262
Add PointLoad by @adtzlr in https://github.com/adtzlr/felupe/pull/263
Rename mask argument to points of point-load by @adtzlr in https://github.com/adtzlr/felupe/pull/266
Add Step, Job and CharacteristicCurve by @adtzlr in https://github.com/adtzlr/felupe/pull/270
Add missing import FieldPlaneStrain by @adtzlr in https://github.com/adtzlr/felupe/pull/272

Full Changelog: https://github.com/adtzlr/felupe/compare/v4.0.0...v5.0.0###

- Python
Published by adtzlr almost 4 years ago

Changelog since [3.1.0] - 2022-05-02

Added

Add SolidBody.evaluate.kirchhoff_stress() method. Contrary to the Cauchy stress method, this gives correct results in incompressible plane stress.
Add SolidBodyTensor for tensor-based material definitions with state variables.
Add bodies argument to newtonrhapson().
Add a container class for fields, FieldContainer (renamed from FieldMixed).
Add len(field) method for FieldContainer (length = number of fields).

Changed

Unify handling of Field and FieldMixed.
Constitutive models use lists as in- and output (consistency between single- and mixed-formulations).
Allow field updates directly from 1d sparse-solved vector without splitted by field-offsets.

Fixed

Fix tovoigt() helper for data with more or less than two trailing axes and 2D tensors.
Fix errors for force() and moment() helpers if the residuals are sparse.

Removed

Remove wrapper for matADi-materials (not necessary with field containers).
Remove IntegralFormMixed and IntegralFormAxisymmetric from global namespace.

Auto-Generated Release-Notes

What's Changed

Fix tovoigt() helper by @adtzlr in https://github.com/adtzlr/felupe/pull/219
Add SolidBody kirchhoff stress by @adtzlr in https://github.com/adtzlr/felupe/pull/221
Add SolidBodyTensor by @adtzlr in https://github.com/adtzlr/felupe/pull/224
Enhance newtonrhapson(): Add body argument by @adtzlr in https://github.com/adtzlr/felupe/pull/226
Newton-Rhapson: Remove body and add bodies argument by @adtzlr in https://github.com/adtzlr/felupe/pull/227
Unify handling of Field and FieldMixed by @adtzlr in https://github.com/adtzlr/felupe/pull/229
Add field-container length method by @adtzlr in https://github.com/adtzlr/felupe/pull/233
Docs: Fix tutorial "shear" by @adtzlr in https://github.com/adtzlr/felupe/pull/235

Full Changelog: https://github.com/adtzlr/felupe/compare/v3.1.0...v4.0.0

- Python
Published by adtzlr almost 4 years ago

Changelog since [3.0.0] - 2022-04-28

Added

Add optional parallel (threaded) basis evaluation and add Form(v, u, parallel=True).
Add mechanics submodule with SolidBody and SolidBodyPressure.

Fixed

Fix matADi materials for (mixed) axisymmetric analyses.
Fix missing radius in axisymmetric integral forms.

Auto-Generated Release-Notes

What's Changed

Enhance Basis: Parallel evaluation by @adtzlr in https://github.com/adtzlr/felupe/pull/212
Add mechanics submodule: SolidBody and SolidBodyPressure by @adtzlr in https://github.com/adtzlr/felupe/pull/216

Full Changelog: https://github.com/adtzlr/felupe/compare/v3.0.0...v3.1.0

- Python
Published by adtzlr about 4 years ago

Changelog since [2.0.1] - 2022-01-11

Added

Add sym argument to Bilinearform.integrate() and Bilinearform.assemble().
Add FieldsMixed which creates a FieldMixed of length n based on a template region.
Add function to mirror a Mesh mesh.mirror().
Add a new parallel assembly that uses a threaded version of np.einsum instead (einsumt).
Add parallel versions of math helpers (dya, cdya, dot, ddot) using einsumt.
Add parallel keyword to constitutive models (NeoHooke, LinearElasticTensorNotation and ThreeFieldVariation).
Add RegionBoundary along with template regions for Quad and Hexahedron and GaussLegendreBoundary.
Add optional normal vector argument for function and gradient methods of AreaChange.
Add a new Mesh-tool triangulate(), applicable on Quad and Hexahedron meshes.
Add a new Mesh-method Mesh.as_meshio().
Add a function decorator @Form(...) for linear and bilinear form objects.

Changed

Enforce consistent arguments for functions inside mesh (points, cells, cell_data or Mesh).
Rename Numba-parallel assembly to jit.
Move single element shape functions and their derivatives from region.h to region.element.h and region.dhdr to region.element.dhdr.
Repeat element shape functions and their derivatives for each cell (as preparation for an upcoming RegionBoundary).
Improve mesh.convert() by using the function decorator @mesh_or_data.
Allow an array to be passed as the expansion arguments of mesh.expand() and mesh.revolve().
Allow optional keyword args to be passed to Mesh.save(**kwargs), acts as a wrapper for Mesh.as_meshio(**kwargs).write().

Fixed

Fix area normal vectors of RegionBoundary.
Fix integration and subsequent assembly of BilinearForm if field and mesh dimensions are not equal.

Auto-Generated Release-Notes

What's Changed

Add sym argument to Bilinearform.integrate() and Bilinearform.assemble() by @adtzlr in https://github.com/adtzlr/felupe/pull/186
Add FieldsMixed by @adtzlr in https://github.com/adtzlr/felupe/pull/190
Update coverage.yml: Add Python 3.10 by @adtzlr in https://github.com/adtzlr/felupe/pull/182
Rename keyword parallel to jit, introduce new parallel, improved Mesh-tools by @adtzlr in https://github.com/adtzlr/felupe/pull/195
Preparations for RegionBoundary by @adtzlr in https://github.com/adtzlr/felupe/pull/196
Add a region on the boundary faces/edges on a mesh by @adtzlr in https://github.com/adtzlr/felupe/pull/197
Improve AreaChange by @adtzlr in https://github.com/adtzlr/felupe/pull/198
Add axisymmetric argument to FieldsMixed by @adtzlr in https://github.com/adtzlr/felupe/pull/199
Add Mesh-Tool: triangulate() by @adtzlr in https://github.com/adtzlr/felupe/pull/202
Improve Mesh-tool convert() by @adtzlr in https://github.com/adtzlr/felupe/pull/203
Improve Mesh by @adtzlr in https://github.com/adtzlr/felupe/pull/204
Fix dimensions of value-array in assembly of BilinearForm by @adtzlr in https://github.com/adtzlr/felupe/pull/206
Add function decorator for weakforms by @adtzlr in https://github.com/adtzlr/felupe/pull/208
Change Forms How-To by @adtzlr in https://github.com/adtzlr/felupe/pull/209
Simplify How-To on Forms by @adtzlr in https://github.com/adtzlr/felupe/pull/210

Full Changelog: https://github.com/adtzlr/felupe/compare/v2.0.1...v3.0.0

- Python
Published by adtzlr about 4 years ago

Changelog since [2.0.0] - 2022-01-10

Fixed

Fixed wrong result of assembly generated by a parallel loop with prange.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.6.0] - 2021-12-02

Added

Add a new method to deepcopy a Mesh with Mesh.copy()
Add broadcasting capability for trailing axes inside the parallel form integrators.
Add Basis on top of a field for virtual fields used in linear and bilinear forms.
Add LinearForm and BilinearForm (including mixed variants) for vector/matrix assembly out of weak form expressions.
Add parallel keyword for threaded integration/assembly of LinearForm and BilinearForm.

Changed

Enhance Boundary for the application of prescribed values of any user-defined Field which is part of FieldMixed.
The whole mixed-field has to be passed to dof.apply() along with the offsets returned from dof.partition for mixed-field formulations.
Set default value shape=(1, 1) for hessian() methods of linear elastic materials.

Fixed

Fixed einstein summation of math.dot() for two vectors with trailing axes.

Removed

Remove dof.extend because dof.partition does not need it anymore.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.5.0] - 2021-11-29

Added

Add LinearElasticPlaneStress and LinearElasticPlaneStrain material formulations.
Add region argument for LinearElastic.hessian().

Changed

Re-formulate LinearElastic materials in terms of the deformation gradient.
Re-formulate LinearElastic material in matrix notation (Speed-up of ~10 for elasticity matrix compared to previous implementation.)
Move previous LinearElastic to constitution.LinearElasticTensorNotation.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.4.0] - 2021-11-15

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.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.3.0]

[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.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.2.0]

[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.

- Python
Published by adtzlr over 4 years ago

Changelog since [1.1.0]

[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

- Python
Published by adtzlr over 4 years ago

Changelog since [1.0.1]

[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.

- Python
Published by adtzlr over 4 years ago

felupe - v1.0.1 Fix install with pip if `sparse` is not installed

Changelog

[1.0.1] - 2021-10-19

Fixed

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

- Python
Published by adtzlr over 4 years ago

felupe - v1.0.0 Code simplifications and faster indices evaluation of Fields

Changelog of [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.

- Python
Published by adtzlr over 4 years ago

felupe - v0.1.1 First BETA release of FElupe

First BETA release of FElupe 🥇

A lot of unused, undocumented or deprecated features have been removed. Test coverage is now > 97%. Docs have been simplified, re-tested and enhanced. The code structure is now easier to understand as all sub-modules are located in subfolders. Constitution has been simplified, Automatic Differentation capability is now outsourced to a separate package matadi.

Have fun 🏖️

- Python
Published by adtzlr over 4 years ago

felupe - v0.0.11 Bugfixes and Improvements

Changelog

Bugfixes: - fix MINI element by using a non-normalized bubble function r*s*t*(1-r-s-t) [1]. This prevents convergence problems.

Enhancements: - allow user-defined bubble-multiplier for MINI element - add point_data and cell_data to save option - add MultiPointContact (no friction), similar to MultiPointConstraint - add new three-field (Grad(u), P, F) template for Materials with automatic differentiation using CasADi - IntegralForm (lineaform): allow arbitrary dimension of fun if grad_v=None (non-parallel only) - add grad argument to Region (default is True) - add FieldMixed class for mixed-field problems - split huge source files into several smaller files

- Python
Published by adtzlr over 4 years ago

felupe - v0.0.10 Enhance IntegralForm classes

Now all IntegralForms have the same arguments.

- Python
Published by adtzlr almost 5 years ago

felupe - v0.0.9 Introduce Automatic Differentation Capability (via casADi)

Introduce Automatic Differentation Capability (via casADi)

felupe.constitution now contains ad.Material, ad.Materialup and ad.MaterialupJ for single-, two- and three-field formulations. Only the strain energy density function has to be defined - stress and elasticity tensors are created by automatic differentiation.

Example usage: Neo-Hookean material model

```python import felupe as fe import casadi as ca

def W(F, mu, bulk): "Strain energy density function for Neo-Hookean material formulation." J = ca.det(F) C = F.T @ F IC = ca.trace(C) return mu / 2 * (IC * J ** (-2 / 3) - 3) + bulk / 2 * (J - 1) ** 2

umat = fe.constitution.ad.Material(W, mu=1.0, bulk=2.0)

P = umat.f(F) A = umat.A(F) ```

- Python
Published by adtzlr almost 5 years ago

Recent Releases of felupe

felupe - v9.3.0

[9.3.0] - 2025-06-29

Added

Changed

Fixed

What's Changed

felupe - v9.2.0

Added

Changed

Fixed

Removed

What's Changed

felupe - v9.1.0

[9.1.0] - 2024-11-23

Added

Changed

Fixed

Removed

What's Changed

felupe - v9.0.0

[9.0.0] - 2024-09-06

Added

Changed

Deprecated

What's Changed

felupe - v8.8.0

[8.8.0] - 2024-06-16

Added

What's Changed

felupe - v8.7.0

[8.7.0] - 2024-06-07

Added

Changed

Fixes

What's Changed

felupe - v8.6.0

[8.6.0] - 2024-05-29

Added

Changed

Fixed

What's Changed

New Contributors

felupe - v8.5.1

[8.5.1] - 2024-05-08

Fixed

What's Changed

felupe - v8.5.0

[8.5.0] - 2024-04-27

Added

What's Changed

felupe - v8.4.0

[8.4.0] - 2024-04-12

Added

Changed

Fixed

What's Changed

felupe - v8.3.1

[8.3.1] - 2024-04-06

Fixed

What's Changed

felupe - v8.3.0

[8.3.0] - 2024-04-02

Added

Changed

Fixed

What's Changed

felupe - v8.2.1

[8.2.1] - 2024-03-30

Fixed

What's Changed

felupe - v8.2.0

[8.2.0] - 2024-03-25

Added

Changed

What's Changed

felupe - v8.1.0

[8.1.0] - 2024-03-23

Added

Changed