[Getfem-commits] r5109 - in /trunk/getfem: doc/sphinx/source/ doc/sphinx/source/userdoc/ interface/tests/matlab/

[Yves . Renard]

Author: renard
Date: Tue Oct 27 17:19:53 2015
New Revision: 5109

URL: http://svn.gna.org/viewcvs/getfem?rev=5109&view=rev
Log:
small improvements of Xfem doc

Modified:
    trunk/getfem/doc/sphinx/source/biblio.rst
    trunk/getfem/doc/sphinx/source/userdoc/xfem.rst
    trunk/getfem/interface/tests/matlab/demo_fictitious_domains.m

Modified: trunk/getfem/doc/sphinx/source/biblio.rst
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/doc/sphinx/source/biblio.rst?rev=5109&r1=5108&r2=5109&view=diff
==============================================================================
--- trunk/getfem/doc/sphinx/source/biblio.rst   (original)
+++ trunk/getfem/doc/sphinx/source/biblio.rst   Tue Oct 27 17:19:53 2015
@@ -54,6 +54,9 @@
 .. [Duan2014] H. Duan,
    *A finite element method for Reissner-Mindlin plates*,
    Math. Comp., 83:286, 701-733, 2014.
+
+.. [Fa-Po-Re2015] M. Fabre, J. Pousin, Y. Renard.
+   *A fictitious domain method for frictionless contact problems in elasticity 
using Nitsche's method*. preprint, 
https://hal.archives-ouvertes.fr/hal-00960996v1
 
 .. [Fa-Pa2003] F. Facchinei and J.-S. Pang,
    *Finite-Dimensional Variational Inequalities and Complementarity Problems, 
Vol. II*,

Modified: trunk/getfem/doc/sphinx/source/userdoc/xfem.rst
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/doc/sphinx/source/userdoc/xfem.rst?rev=5109&r1=5108&r2=5109&view=diff
==============================================================================
--- trunk/getfem/doc/sphinx/source/userdoc/xfem.rst     (original)
+++ trunk/getfem/doc/sphinx/source/userdoc/xfem.rst     Tue Oct 27 17:19:53 2015
@@ -36,7 +36,8 @@
 The tools for Xfem have been then enriched by the PhD works
 of J. Larsy (see for instance [LA-RE-SA2010]_) the one
 of E. Chahine (see for instance [CH-LA-RE2011]_, [NI-RE-CH2011]_),
-and of S. Amdouni  (see for instance [AM-MO-RE2014]_, [AM-MO-RE2014b]_).
+of S. Amdouni  (see for instance [AM-MO-RE2014]_, [AM-MO-RE2014b]_)
+and of M. Fabre (see for instance  [Fa-Po-Re2015]_).
 
 
 .. important::
@@ -171,12 +172,11 @@
 Cut-fem
 -------
 
-The previous tools (description of the level-sets, cut mesh and cut 
integration method) are already sufficient to build cut-fem applications.
-
-+ exemple (:file:`interface/tests/python/demo_fictitious_domain.py`)
-
-A few word on stabilizations methods ... 
-
+The implementation of a cut finite element method such as described in 
[bu-ha2010]_, i.e. a finite element on a fictitious domain restricted to a 
smaller real domain, is possible just using the previous tools and mainly the 
adapted integration method. Several examples are available on |gf| test 
programs. See for instance 
:file:`interface/tests/python/demo_fictitious_domain.py` or 
:file:`interface/tests/matlab/demo_fictitious_domain.m`.
+
+In this context, one often needs to restrict the unknown finite element field 
to the degrees of freedom whose corresponding shape function supports have an 
intersection with the real domain. This can be done using the 
``partial_mesh_fem`` object. See for instance 
:file:`interface/tests/matlab/demo_structural_optimization.m`.
+
+Note that often, a stabilization technique have to be considered in order to 
treat eventual locking phenomena due to element with very small intersection 
with the real domain for example when applying a Dirichlet condition. See for 
instance [bu-ha2010]_,  [HA-RE2009]_ and [Fa-Po-Re2015]_.
 
 
 Discontinuous field across some level-sets
@@ -197,16 +197,31 @@
 To initialize the object or to actualize it when the value of the level-set
 function is modified, one has to call the method ``mfls.adapt()``.
 
-To represent discontinuous fields, the finite element method is enriched with
-discontinuous functions which are the product of a Heaviside function by the 
base
-functions of the finite element method represented by ``mf`` (see [Xfem]_ for
-more details).
-
-
-Fictitious domain approach with Xfem
-------------------------------------
-
-An example of a Poisson problem with a Dirichlet condition posed on a boundary
-independant of the mesh is present on the ``tests`` directory of the 
distribution.
-
-See :file:`contrib/xfem_contact/xfem_dirichlet.cc` file.
+To represent discontinuous fields, the finite element method is enriched
+with discontinuous functions which are the product of a Heaviside function
+by the base functions of the finite element method represented by ``mf``
+(see [HA-HA2004]_ and [Xfem]_ for more details).
+
+
+Xfem
+----
+
+mesh_fem_product, mesh_fem_sum, Xfem_plus, Xfem_minus ...
+
+
+
+
+Note that some procedures are available in the file 
:file:`src/getfem/getfem_crack_sif.h` to compute the stress intensity factors 
in 2D (restricted to homogeneous isotropic linearized elasticity).
+
+
+
+Post treatment
+--------------
+
+Several tools are available to represent the solution only on a side of a 
levels-set or on both taking into account the discontinuity (for Xfem 
approximation).
+
+
+Slices ...
+
+
+cut-mesh, interpolation

Modified: trunk/getfem/interface/tests/matlab/demo_fictitious_domains.m
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/tests/matlab/demo_fictitious_domains.m?rev=5109&r1=5108&r2=5109&view=diff
==============================================================================
--- trunk/getfem/interface/tests/matlab/demo_fictitious_domains.m       
(original)
+++ trunk/getfem/interface/tests/matlab/demo_fictitious_domains.m       Tue Oct 
27 17:19:53 2015
@@ -105,7 +105,7 @@
 cv_out = get(mim, 'convex_index');
 cv_in = setdiff(gf_mesh_get(m, 'cvid'), cv_out);
 
-mfu = gfMeshFem('partial', mfu0, dof_out, cv_in);
+% mfu = gfMeshFem('partial', mfu0, dof_out, cv_in);
 
 
 md=gf_model('real');