finite element software developers

Abstract: Nearly incompressible materials arise in soft tissues, elastomers, gels, and other soft solids. Yet, numerical modeling remains challenging because displacement-based formulations can exhibit volumetric locking near incompressibility. Mixed displacement-pressure formulations address this difficulty but require compatible approximation spaces satisfying the inf-sup condition. In mixed formulations, discontinuous pressure spaces are attractive since pressure degrees of freedom can be condensed locally. In high-order isogeometric analysis, however, their Galerkin stability requires pressure-grid coarsening that increases with polynomial degree. Here, we develop a variational multiscale (VMS) stabilized mixed isogeometric formulation that retains this locality while substantially relaxing the restriction. The key observation is that VMS stabilization controls nonconstant elementwise pressure modes and is unconditionally stable with respect to its parameter, provided the elementwise constant pressure component satisfies the inf-sup condition. Consequently, for maximally continuous displacement splines, a pressure grid coarsened by a factor of two satisfies the constant-mode inf-sup criterion independently of polynomial degree. Local condensation then yields a displacement-only system, facilitating use in existing displacement solvers. The formulation is further extended to immersed isogeometric analysis with extended B-splines, pressure-cell aggregation, subdivision-based cut-cell quadrature, and weak enforcement of essential boundary conditions. Discrete inf-sup studies and convergence analyses in the Kirsch and Cook's membrane benchmarks demonstrate stability and show that the immersed treatment preserves mixed stability properties. The proposed framework therefore provides a local, high-order, and robust approach for nearly incompressible problems in mesh-conforming and immersed settings.

Abstract: Cumulative energy fields in irreversible evolution depend on the loading path, even when terminal states and global work agree. We derive a posteriori bounds for these fields in scalar gradient damage with a prescribed driving density. A perturbation estimate for the incremental obstacle map yields a continuous loading history, which is enclosed by nondecreasing residual supersolutions. A primitive identity transfers this enclosure to the cumulative field. An exact endpoint-rule defect and outward arithmetic certify the output delivered by the numerical method. The analysis separates spatial, loading, and algebraic errors and gives sufficient conditions for certificate decay. Computable realizations use splines and verified supersolutions in one dimension and on two-dimensional rectangles. The production bounds cover the independently enclosed errors of all 417 saved histories from 28 adaptive calculations and all 25 fixed-parameter controls. In a narrow-localization case, the verified relative error decreases from approximately 158.78\% to below 2.125\% at a 5\% certified target. Source-oriented marking reduces accumulated matrix--vector work by 35.9\% and 51.9\% in two comparisons. A mechanically generated two-dimensional benchmark shows how auxiliary supersolutions sharpen constant barriers and enable acceptance of both loading orders.