From: frohne Date: Mon, 11 Feb 2013 18:18:45 +0000 (+0000) Subject: add a section adaptive mesh refinement and in implementation a passage for the Input... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6b8ba4140609513fa78852feb5c41cd7f4a6debe;p=dealii-svn.git add a section adaptive mesh refinement and in implementation a passage for the Input class git-svn-id: https://svn.dealii.org/trunk@28314 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-42/doc/intro.dox b/deal.II/examples/step-42/doc/intro.dox index beb68024ea..3b221f2348 100644 --- a/deal.II/examples/step-42/doc/intro.dox +++ b/deal.II/examples/step-42/doc/intro.dox @@ -315,6 +315,15 @@ using a primal-dual active set strategy and algebraic multigrid, Int. J. Numer. Meth. Engng, 2007, 69, pp. 524-543. +

Adaptive mesh refinement

+ +Since we run our program in 3d, there is a good reason to use adaptive +mesh refined. To make life a bit easier we are choosing the +KellyErrorEstimator that is already implemented in deal.II. We hand the +solution vector to it which contains the displacement $u$. As we will see in the +results it yields a quite reasonable adaptive mesh for the contact zone. + +

Implementation

This tutorial is essentailly a mixture of step-40 and step-41 but instead of @@ -327,5 +336,6 @@ Another difficulty is the handling of the different constraints from arises from the contact. For this purpose we create three objects of type ConstraintMatrix. - - +Beside the ConstitutiveLaw class there is another new class called Input. This +class allows us to read in an obstacle from a file. In our example the file +'obstacle\_file.dat' containts data which describe the Chinese symbol for force.