From 45f738ff5c4a2321f807c961f80d9c1c8741f216 Mon Sep 17 00:00:00 2001 From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de> Date: Wed, 2 Dec 2020 10:30:54 +0100 Subject: [PATCH] Correct parset file to calculate the deformation for film-on-substrate Corrected the formula in the code before; now this also matches the value we publish in the paper. --- problems/film-on-substrate-deformation.parset | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/problems/film-on-substrate-deformation.parset b/problems/film-on-substrate-deformation.parset index 0a1c371..f4f147e 100644 --- a/problems/film-on-substrate-deformation.parset +++ b/problems/film-on-substrate-deformation.parset @@ -97,18 +97,17 @@ mooneyrivlin_12 = 0 mooneyrivlin_03 = 0 # volume-preserving parameter -mooneyrivlin_k = 75e+6 +mooneyrivlin_k = 90e+6 # How to choose the volume-preserving parameter? # We need a stretch of 30% (45e4 10e4 2e4 in x-direction, so a stretch of 45e4*0.3 = 13.5e4) -# 184 2:1, mooneyrivlin_k = .. and mooneyrivlin_energy = square, neumannValues = 27e4 0 0 -# 184 2:1, mooneyrivlin_k = 57e+6 and mooneyrivlin_energy = log, neumannValues = 27e4 0 0 +# 184 2:1, mooneyrivlin_k = 90e+6 approximately (depends also on the number of grid levels!) and mooneyrivlin_energy = square or log, neumannValues = 27e4 0 0 mooneyrivlin_energy = square # log, square or ciarlet; different ways to compute the Mooney-Rivlin-Energy # ciarlet: Fomula from "Ciarlet: Three-Dimensional Elasticity", here no penalty term is # log: Generalized Rivlin model or polynomial hyperelastic model, using 0.5*mooneyrivlin_k*log(det(∇φ)) as the volume-preserving penalty term -# square: Generalized Rivlin model or polynomial hyperelastic model, using mooneyrivlin_k*(det(∇φ)-1)² as the volume-preserving penalty term +# square: Generalized Rivlin model or polynomial hyperelastic model, using 0.5*mooneyrivlin_k*(det(∇φ)-1)² as the volume-preserving penalty term [] -- GitLab