[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] r5109 - in /trunk/getfem: doc/sphinx/source/ doc/sphinx
From: |
Yves . Renard |
Subject: |
[Getfem-commits] r5109 - in /trunk/getfem: doc/sphinx/source/ doc/sphinx/source/userdoc/ interface/tests/matlab/ |
Date: |
Tue, 27 Oct 2015 16:19:54 -0000 |
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');
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-commits] r5109 - in /trunk/getfem: doc/sphinx/source/ doc/sphinx/source/userdoc/ interface/tests/matlab/,
Yves . Renard <=