GetFEM  5.4.2
getfem_linearized_plates.h
Go to the documentation of this file.
1 // /* -*- c++ -*- (enables emacs c++ mode) */
2 /*===========================================================================
3 
4  Copyright (C) 2004-2020 Yves Renard, Jeremie Lasry
5 
6  This file is a part of GetFEM
7 
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.
20 
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.
29 
30 ===========================================================================*/
31 /**@file getfem_linearized_plates.h
32  @author Yves Renard <Yves.Renard@insa-lyon.fr>
33  @date November 1, 2004.
34  @brief Reissner-Mindlin plate model brick.
35 */
36 
37 #ifndef GETFEM_LINEARIZED_PLATES_H__
38 #define GETFEM_LINEARIZED_PLATES_H__
39 
40 #include "getfem_models.h"
41 
42 
43 namespace getfem {
44 
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);
50 
51  /** Add the elementary transformation corresponding to the projection
52  on P0 element.
53  The name is the name given to the elementary transformation.
54  */
55  void add_P0_projection(model &md, std::string name);
56 
57 
58  /** Add a term corresponding to the classical Reissner-Mindlin plate
59  model for which `u3` is the transverse displacement,
60  `Theta` the rotation of
61  fibers normal to the midplane, 'param_E' the Young Modulus,
62  `param_nu` the poisson ratio,
63  `param_epsilon` the plate thickness,
64  `param_kappa` the shear correction factor. Note that since this brick
65  uses the high level generic assembly language, the parameter can
66  be regular expression of this language.
67  There are three variants.
68  `variant = 0` corresponds to the an
69  unreduced formulation and in that case only the integration
70  method `mim` is used. Practically this variant is not usable since
71  it is subject to a strong locking phenomenon.
72  `variant = 1` corresponds to a reduced integration where `mim` is
73  used for the rotation term and `mim_reduced` for the transverse
74  shear term. `variant = 2` (default) corresponds to the projection onto
75  a rotated RT0 element of the transverse shear term. For the moment, this
76  is adapted to quadrilateral only (because it is not sufficient to
77  remove the locking phenomenon on triangle elements). Note also that if
78  you use high order elements, the projection on RT0 will reduce the order
79  of the approximation.
80  Returns the brick index in the model.
81  */
83  (model &md, const mesh_im &mim, const mesh_im &mim_reduced,
84  const std::string &u3,
85  const std::string &Theta, const std::string &param_E,
86  const std::string &param_nu, const std::string &param_epsilon,
87  const std::string &param_kappa, size_type variant = size_type(2),
88  size_type region = size_type(-1));
89 
90  /** Add a term corresponding to the enriched Reissner-Mindlin plate
91  model for which `varname_ua` is the membrane displacements,
92  `varname_u3` is the transverse displacement,
93  `varname_theta` the rotation of
94  fibers normal to the midplane,
95  `varname_theta3` the pinching,
96  'param_E' the Young Modulus,
97  `param_nu` the poisson ratio,
98  `param_epsilon` the plate thickness. Note that since this brick
99  uses the high level generic assembly language, the parameter can
100  be regular expression of this language.
101  There are four variants.
102  `variant = 0` corresponds to the an
103  unreduced formulation and in that case only the integration
104  method `mim` is used. Practically this variant is not usable since
105  it is subject to a strong locking phenomenon.
106  `variant = 1` corresponds to a reduced integration where `mim` is
107  used for the rotation term and `mim_reduced1` for the transverse
108  shear term and `mim_reduced2` for the pinching term.
109  `variant = 2` (default) corresponds to the projection onto
110  a rotated RT0 element of the transverse shear term and a reduced integration for the pinching term.
111  For the moment, this is adapted to quadrilateral only (because it is not sufficient to
112  remove the locking phenomenon on triangle elements). Note also that if
113  you use high order elements, the projection on RT0 will reduce the order
114  of the approximation.
115  `variant = 3` corresponds to the projection onto
116  a rotated RT0 element of the transverse shear term and the projection onto P0 element of the pinching term.
117  For the moment, this is adapted to quadrilateral only (because it is not sufficient to
118  remove the locking phenomenon on triangle elements). Note also that if
119  you use high order elements, the projection on RT0 will reduce the order
120  of the approximation.
121  Returns the brick index in the model.
122  */
124  (model &md, const mesh_im &mim, const mesh_im &mim_reduced1, const mesh_im &mim_reduced2,
125  const std::string &ua,const std::string &Theta,
126  const std::string &u3,const std::string &Theta3,
127  const std::string &param_E, const std::string &param_nu, const std::string &param_epsilon,
128  size_type variant = size_type(3), size_type region = size_type(-1));//size_type(2)
129 
130 
131 
132 
133 
134 
135 } /* end of namespace getfem. */
136 
137 
138 #endif /* GETFEM_LINEARIZED_PLATES_H__ */
getfem::add_Mindlin_Reissner_plate_brick
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...
Definition: getfem_linearized_plates.cc:28
getfem::add_enriched_Mindlin_Reissner_plate_brick
size_type add_enriched_Mindlin_Reissner_plate_brick(model &md, const mesh_im &mim, const mesh_im &mim_reduced1, const mesh_im &mim_reduced2, const std::string &ua, const std::string &Theta, const std::string &u3, const std::string &Theta3, const std::string &param_E, const std::string &param_nu, const std::string &param_epsilon, size_type variant=size_type(3), size_type region=size_type(-1))
Add a term corresponding to the enriched Reissner-Mindlin plate model for which varname_ua is the mem...
Definition: getfem_linearized_plates.cc:89
bgeot::size_type
size_t size_type
used as the common size type in the library
Definition: bgeot_poly.h:49
getfem::add_P0_projection
void add_P0_projection(model &md, std::string name)
Add the elementary transformation corresponding to the projection on P0 element.
Definition: getfem_linearized_plates.cc:395
getfem
GEneric Tool for Finite Element Methods.
Definition: getfem_accumulated_distro.h:46
getfem_models.h
Model representation in Getfem.
getfem::add_2D_rotated_RT0_projection
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...
Definition: getfem_linearized_plates.cc:288