GetFEM++  5.3
Go to the documentation of this file.
1 /* -*- c++ -*- (enables emacs c++ mode) */
2 /*===========================================================================
4  Copyright (C) 2004-2017 Yves Renard, Jeremie Lasry
6  This file is a part of GetFEM++
8  GetFEM++ is free software; you can redistribute it and/or modify it
9  under the terms of the GNU Lesser General Public License as published
10  by the Free Software Foundation; either version 3 of the License, or
11  (at your option) any later version along with the GCC Runtime Library
12  Exception either version 3.1 or (at your option) any later version.
13  This program is distributed in the hope that it will be useful, but
14  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15  or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
16  License and GCC Runtime Library Exception for more details.
17  You should have received a copy of the GNU Lesser General Public License
18  along with this program; if not, write to the Free Software Foundation,
19  Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
21  As a special exception, you may use this file as it is a part of a free
22  software library without restriction. Specifically, if other files
23  instantiate templates or use macros or inline functions from this file,
24  or you compile this file and link it with other files to produce an
25  executable, this file does not by itself cause the resulting executable
26  to be covered by the GNU Lesser General Public License. This exception
27  does not however invalidate any other reasons why the executable file
28  might be covered by the GNU Lesser General Public License.
30 ===========================================================================*/
31 /**@file getfem_linearized_plates.h
32  @author Yves Renard <>
33  @date November 1, 2004.
34  @brief Reissner-Mindlin plate model brick.
35 */
40 #include "getfem_models.h"
43 namespace getfem {
45  /** Add the elementary transformation corresponding to the projection
46  on rotated RT0 element for two-dimensional elements to the model.
47  The name is the name given to the elementary transformation.
48  */
49  void add_2D_rotated_RT0_projection(model &md, std::string name);
52  /** Add a term corresponding to the classical Reissner-Mindlin plate
53  model for which `u3` is the transverse displacement,
54  `Theta` the rotation of
55  fibers normal to the midplane, 'param_E' the Young Modulus,
56  `param_nu` the poisson ratio,
57  `param_epsilon` the plate thickness,
58  `param_kappa` the shear correction factor. Note that since this brick
59  uses the high level generic assembly language, the parameter can
60  be regular expression of this language.
61  There are three variants.
62  `variant = 0` corresponds to the an
63  unreduced formulation and in that case only the integration
64  method `mim` is used. Practically this variant is not usable since
65  it is subject to a strong locking phenomenon.
66  `variant = 1` corresponds to a reduced integration where `mim` is
67  used for the rotation term and `mim_reduced` for the transverse
68  shear term. `variant = 2` (default) corresponds to the projection onto
69  a rotated RT0 element of the transverse shear term. For the moment, this
70  is adapted to quadrilateral only (because it is not sufficient to
71  remove the locking phenomenon on triangle elements). Note also that if
72  you use high order elements, the projection on RT0 will reduce the order
73  of the approximation.
74  Returns the brick index in the model.
75  */
77  (model &md, const mesh_im &mim, const mesh_im &mim_reduced,
78  const std::string &u3,
79  const std::string &Theta, const std::string &param_E,
80  const std::string &param_nu, const std::string &param_epsilon,
81  const std::string &param_kappa, size_type variant = size_type(2),
82  size_type region = size_type(-1));
84 } /* end of namespace getfem. */
size_t size_type
used as the common size type in the library
Definition: bgeot_poly.h:49
Model representation in Getfem.
GEneric Tool for Finite Element Methods.
size_type add_Mindlin_Reissner_plate_brick(model &md, const mesh_im &mim, const mesh_im &mim_reduced, const std::string &u3, const std::string &Theta, const std::string &param_E, const std::string &param_nu, const std::string &param_epsilon, const std::string &param_kappa, size_type variant=size_type(2), size_type region=size_type(-1))
Add a term corresponding to the classical Reissner-Mindlin plate model for which u3 is the transverse...
void add_2D_rotated_RT0_projection(model &md, std::string name)
Add the elementary transformation corresponding to the projection on rotated RT0 element for two-dime...