Linear elasticity with a non-penetration surface condition

バージョン 1.0.0 (4.41 MB) 作成者: Jan
A 3D FEM solver based on an energy minimization formulated in C^0 deformation including a new surface non-penetration term.
ダウンロード: 9
更新 2023/2/1


A complementary code of the article
Stefan Krömer, Jan Valdman - Surface penalization of self-interpenetration in linear and nonlinear elasticity
Description: a 3D linear elasticity is formulated as an energy minimization and combined with the (nonlocal) self-interpenetration term defined on the deforming domain boundary
The code is based on modification of benchmark 4 from the article
1) A. Moskovka1 and J. Valdman - Fast MATLAB evaluation of nonlinear energies using FEM in 2D and 3D: Nodal elements
Vectorization/Solver of linear elasticity is taken from the article
2) T. Rahman, J. Valdman - Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
Mesh visualizations are taken from the article
3) M. Čermák, S. Sysala, J. Valdman - Efficient and flexible MATLAB implementation of 2D and 3D elastoplastic problems
Evaluation of the self-interpenetration term is modified from the article
4) S. Krömer, J. Valdman, Global injectivity in second-gradient Nonlinear Elasticity and its approximation with penalty terms
3D pincers mesh was created by A. Moskovka.
Please cite the paper if you find the code useful. A link to the paper can be found at the author web page located at
Available scripts:


Jan (2024). Linear elasticity with a non-penetration surface condition (, MATLAB Central File Exchange. 取得済み .

MATLAB リリースの互換性
作成: R2022b
R2021a 以降 R2022b 以前と互換性あり
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!










バージョン 公開済み リリース ノート