Source author record

Ondřej Rokoš

Ondřej Rokoš appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci physics.comp-ph Computational Engineering, Finance, and Science cond-mat.soft Machine Learning math.AP

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Homogenization with Guaranteed Bounds via Primal-Dual Physically Informed Neural Networks

Physics-informed neural networks (PINNs) have shown promise in solving partial differential equations (PDEs) relevant to multiscale modeling, but they often fail when applied to materials with discontinuous coefficients, such as media with piecewise constant properties. This paper introduces a dual formulation for the PINN framework to improve the reliability of the homogenization of periodic thermo-conductive composites, for both strong and variational (weak) formulations. The dual approach facilitates the derivation of guaranteed upper and lower error bounds, enabling more robust detection of PINN failure. We compare standard PINNs applied to smoothed material approximations with variational PINNs (VPINNs) using both spectral and neural network-based test functions. Our results indicate that while strong-form PINNs may outperform VPINNs in controlled settings, they are sensitive to material discontinuities and may fail without clear diagnostics. In contrast, VPINNs accommodate piecewise constant material parameters directly but require careful selection of test functions to avoid instability. Dual formulation serves as a reliable indicator of convergence quality, and its integration into PINN frameworks enhances their applicability to homogenization problems in micromechanics.

preprint2022arXiv

Learning constitutive models from microstructural simulations via a non-intrusive reduced basis method

In order to optimally design materials, it is crucial to understand the structure-property relations in the material by analyzing the effect of microstructure parameters on the macroscopic properties. In computational homogenization, the microstructure is thus explicitly modeled inside the macrostructure, leading to a coupled two-scale formulation. Unfortunately, the high computational costs of such multiscale simulations often render the solution of design, optimization, or inverse problems infeasible. To address this issue, we propose in this work a non-intrusive reduced basis method to construct inexpensive surrogates for parametrized microscale problems; the method is specifically well-suited for multiscale simulations since the coupled simulation is decoupled into two independent problems: (1) solving the microscopic problem for different (loading or material) parameters and learning a surrogate model from the data; and (2) solving the macroscopic problem with the learned material model. The proposed method has three key features. First, the microscopic stress field can be fully recovered. Second, the method is able to accurately predict the stress field for a wide range of material parameters; furthermore, the derivatives of the effective stress with respect to the material parameters are available and can be readily utilized in solving optimization problems. Finally, it is more data efficient, i.e. requiring less training data, as compared to directly performing a regression on the effective stress. For the microstructures in the two test problems considered, the mean approximation error of the effective stress is as low as 0.1% despite using a relatively small training dataset. Embedded into the macroscopic problem, the reduced order model leads to an online speed up of approximately three orders of magnitude while maintaining a high accuracy as compared to the FE$^2$ solver.

preprint2020arXiv

Reduced integration schemes in micromorphic computational homogenization of elastomeric mechanical metamaterials

Exotic behaviour of mechanical metamaterials often relies on an internal transformation of the underlying microstructure triggered by its local instabilities, rearrangements, and rotations. Depending on the presence and magnitude of such a transformation, effective properties of a metamaterial may change significantly. To capture this phenomenon accurately and efficiently, homogenization schemes are required that reflect microstructural as well as macro-structural instabilities, large deformations, and non-local effects. To this end, a micromorphic computational homogenization scheme has recently been developed, which employs the particular microstructural transformation as a non-local mechanism, magnitude of which is governed by an additional coupled partial differential equation. Upon discretizing the resulting problem it turns out that the macroscopic stiffness matrix requires integration of macro-element basis functions as well as their derivatives, thus calling for a higher-order integration rules. Because evaluation of constitutive law in multiscale schemes involves an expensive solution of a non-linear boundary value problem, computational efficiency can be improved by reducing the number of integration points. Therefore, the goal of this paper is to investigate reduced-order schemes in computational homogenization, with emphasis on the stability of the resulting elements. In particular, arguments for lowering the order of integration from the expensive mass-matrix to a cheaper stiffness-matrix equivalent are first outlined. An efficient one-point integration quadrilateral element is then introduced and proper hourglass stabilization discussed. Performance of the resulting set of elements is finally tested on a benchmark bending example, showing that we achieve accuracy comparable to the full quadrature rules.

preprint2016arXiv

A Variational Formulation of Dissipative Quasicontinuum Methods

Lattice systems and discrete networks with dissipative interactions are successfully employed as meso-scale models of heterogeneous solids. As the application scale generally is much larger than that of the discrete links, physically relevant simulations are computationally expensive. The QuasiContinuum (QC) method is a multiscale approach that reduces the computational cost of direct numerical simulations by fully resolving complex phenomena only in regions of interest while coarsening elsewhere. In previous work (Beex et al., J. Mech. Phys. Solids 64, 154-169, 2014), the originally conservative QC methodology was generalized to a virtual-power-based QC approach that includes local dissipative mechanisms. In this contribution, the virtual-power-based QC method is reformulated from a variational point of view, by employing the energy-based variational framework for rate-independent processes (Mielke and Roubíček, Rate-Independent Systems: Theory and Application, Springer-Verlag, 2015). By construction it is shown that the QC method with dissipative interactions can be expressed as a minimization problem of a properly built energy potential, providing solutions equivalent to those of the virtual-power-based QC formulation. The theoretical considerations are demonstrated on three simple examples. For them we verify energy consistency, quantify relative errors in energies, and discuss errors in internal variables obtained for different meshes and two summation rules.

preprint2012arXiv

Localization analysis of variationally based gradient plasticity model

The paper presents analytical or semi-analytical solutions for the formation and evolution of localized plastic zone in a uniaxially loaded bar with variable cross-sectional area. A variationally based formulation of explicit gradient plasticity with linear softening is used, and the ensuing jump conditions and boundary conditions are discussed. Three cases with different regularity of the stress distribution are considered, and the problem is converted to a dimensionless form. Relations linking the load level, size of the plastic zone, distribution of plastic strain and plastic elongation of the bar are derived and compared to another, previously analyzed gradient formulation.