BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Jan Valdman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We formulate a large-strain model of single-slip crystal elastoplasticity in the framework of energetic solutions. Numerical performance of the model is compared with lab experiments on the compression of a stack of note papers.
Nonlinear energy functionals appearing in the calculus of variations can be discretized by the finite element (FE) method and formulated as a sum of energy contributions from local elements. A fast evaluation of energy functionals containing the first order gradient terms is a central part of this contribution. We describe a vectorized implementation using the simplest linear nodal (P1) elements in which all energy contributions are evaluated all at once without the loop over triangular or tetrahedral elements. Furthermore, in connection to the first-order optimization methods, the discrete gradient of energy functional is assembled in a way that the gradient components are evaluated over all degrees of freedom all at once. The key ingredient is the vectorization of exact or approximate energy gradients over nodal patches. It leads to a time-efficient implementation at higher memory-cost. Provided codes in MATLAB related to 2D/3D hyperelasticity and 2D p-Laplacian problem are available for download and structured in a way it can be easily extended to other types of vector or scalar forms of energies.
A simple MATLAB implementation of hierarchical shape functions on 2D rectangles is explained and available for download. Global shape functions are ordered for a given polynomial degree according to the indices of the nodes, edges, or elements to which they belong. For a uniform p-refinement, the hierarchical structure enables an effective assembly of mass and stiffness matrices. A solution of a boundary value problem is approximated for various levels of uniform h and p refinements.
Minimization of energy functionals is based on a discretization by the finite element method and optimization by the trust-region method. A key tool is a local evaluation of the approximated gradients together with sparsity of the resulting Hessian matrix. We describe a vectorized MATLAB implementation of the p-Laplace problem in one and two space-dimensions, however it is easily applicable to other energy formulations.
An equilibrium of a linear elastic body subject to loading and satisfying the friction and contact conditions can be described by a variational inequality of the second kind and the respective discrete model attains the form of a generalized equation. To its numerical solution we apply the semismooth* Newton method by Gfrerer and Outrata (2019) in which, in contrast to most available Newton-type methods for inclusions, one approximates not only the single-valued but also the multi-valued part. This is performed on the basis of limiting (Morduchovich) coderivative. In our case of the Tresca friction, the multi-valued part amounts to the subdifferential of a convex function generated by the friction and contact conditions. The full 3D discrete problem is then reduced to the contact boundary. Implementation details of the semismooth* Newton method are provided and numerical tests demonstrate its superlinear convergence and mesh independence.
Rahman and Valdman (2013) introduced a new vectorized way to assemble finite element matrices. We utilize underlying vectorization concepts and extend MATLAB codes to implementation of Bogner-Fox-Schmit C1 rectangular elements in 2D. Our focus is on the detailed construction of elements and simple computer demonstrations including energies evaluations and their visualizations.
Rahman and Valdman (2013) introduced a vectorized way to assemble finite element stiffness and mass matrices in MATLAB. Local element matrices are computed all at once by array operations and stored in multi-dimentional arrays (matrices). We build some iterative solvers on available multi-dimentional structures completely avoiding the use of a sparse matrix.
The paper starts with a concise description of the recently developed semismooth* Newton method for the solution of general inclusions. This method is then applied to a class of variational inequalities of the second kind. As a result, one obtains an implementable algorithm exhibiting a local superlinear convergence. Thereafter we suggest several globally convergent hybrid algorithms in which one combines the semismooth* Newton method with selected splitting algorithms for the solution of monotone variational inequalities. Their efficiency is documented by extensive numerical experiments.
We give some theoretical as well as computational results on Laplace and Maxwell constants. Besides the classical de Rham complex we investigate the complex of elasticity and the complex related to the biharmonic equation and general relativity as well using the general functional analytical concept of Hilbert complexes. We consider mixed boundary conditions and bounded Lipschitz domains of arbitrary topology. Our numerical aspects are presented by examples for the de Rham complex in 2D and 3D which not only confirm our theoretical findings but also indicate some interesting conjectures.
This paper is concerned with the two--phase obstacle problem, a type of a variational free boundary problem. We recall the basic estimates of Repin and Valdman (2015) and verify them numerically on two examples in two space dimensions. A solution algorithm is proposed for the construction of the finite element approximation to the two--phase obstacle problem. The algorithm is not based on the primal (convex and nondifferentiable) energy minimization problem but on a dual maximization problem formulated for Lagrange multipliers. The dual problem is equivalent to a quadratic programming problem with box constraints. The quality of approximations is measured by a functional a posteriori error estimate which provides a guaranteed upper bound of the difference of approximated and exact energies of the primal minimization problem. The majorant functional in the upper bound contains auxiliary variables and it is optimized with respect to them to provide a sharp upper bound. A space density of the nonlinear related part of the majorant functional serves as an indicator of the free boundary.
We study computational behavior of a mesoscopic model describing temperature/external magnetic field-driven evolution of magnetization. Due to nonconvex anisotropy energy describing magnetic properties of a body, magnetization can develop fast spatial oscillations creating complicated microstructures. These microstructures are encoded in Young measures, their first moments then identify macroscopic magnetization. Our model assumes that changes of magnetization can contribute to dissipation and, consequently, to variations of the body temperature affecting the length of magnetization vectors. In the ferromagnetic state, minima of the anisotropic energy density depend on temperature and they tend to zero as we approach the so-called Curie temperature. This brings the specimen to a paramagnetic state. Such a thermo-magnetic model is fully discretized and tested on two-dimensional examples. Computational results qualitatively agree with experimental observations. The own MATLAB code used in our simulations is available for download.
We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time.
The quasistatic, Prandtl-Reuss perfect plasticity at small strains is combined with a gradient, reversible (i.e. admitting healing) damage which influences both the elastic moduli and the yield stress. Existence of weak solutions of the resulted system of variational inequalities is proved by a suitable fractional-step discretisation in time with guaranteed numericalstability and convergence. After finite-element approximation, this scheme is computationally implemented and illustrative 2-dimensional simulations are performed. The model allows e.g. for application in geophysical modelling of re-occurring rupture of lithospheric faults. Resulted incremental problems are solved in MATLAB by quasi-Newton method to resolve elastoplasticity component of the solution while damage component is obtained by solution of a quadratic programming problem.
The quasistatic rate-independent damage combined with linearized plasticity with hardening at small strains is investigated. The fractional-step time discretisation is devised with the purpose to obtain a numerically efficient scheme converging possibly to a physically relevant stress-driven solutions, which however is to be verified a-posteriori by using a suitable integrated variant of the maximum-dissipation principle. Gradient theories both for damage and for plasticity are considered to make the scheme numerically stable with guaranteed convergence within the class of weak solutions. After finite-element approximation, this scheme is computationally implemented and illustrative 2-dimensional simulations are performed.
We verify functional a posteriori error estimate proposed by S. Repin for a class of obstacle problems. The obstacle problem is formulated as a quadratic minimization problem with constrains equivalently formulated as a variational inequality. New benchmarks with known analytical solutions in 2D are constructed based on 1D benchmark introduced by P. Harasim and J. Valdman. Numerical approximation of the obstacle problem is obtained by the finite element method using bilinear elements on a rectangular mesh. Error of the approximation is meassured in the energy norm and bounded from above by a functional majorant, whose value is minimized with respect to unknown gradient field discretized by Raviart-Thomas elements and Lagrange multipliers field discretized by piecewise constant functions.
We verify functional a posteriori error estimate for obstacle problem proposed by Repin. Simplification into 1D allows for the construction of a nonlinear benchmark for which an exact solution of the obstacle problem can be derived. Quality of a numerical approximation obtained by the finite element method is compared with the exact solution and the error of approximation is bounded from above by a majorant error estimate. The sharpness of the majorant error estimate is discussed.