[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

If you use this software, please cite it using these metadata.