1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
% Copyright (C) 2017-2020 Yves Renard.
%
% This file is a part of GetFEM++
%
% GetFEM++ is free software; you can redistribute it and/or modify it
% under the terms of the GNU Lesser General Public License as published
% by the Free Software Foundation; either version 3 of the License, or
% (at your option) any later version along with the GCC Runtime Library
% Exception either version 3.1 or (at your option) any later version.
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
% or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
% License and GCC Runtime Library Exception for more details.
% You should have received a copy of the GNU Lesser General Public License
% along with this program; if not, write to the Free Software Foundation,
% Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
% -*- matlab -*- (enables emacs matlab mode)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% parameters for program static Coulomb friction problem %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MU = 0.83E+5; % Lamé coefficients
LAMBDA = 1.18E+5; % (in N/mm^2)
ROT_ANGLE = -1.5E-2; % Rotation angle of the first gear (in radians)
RESIDUAL = 1E-6; % residual for Newton
FRICTION_COEFFICIENT = 0.0E+0; % friction coefficient
MESHNAME_GEAR1 = 'gmsh:./gear1.msh'; % Meshes for gears 1 and 2
MESHNAME_GEAR2 = 'gmsh:./gear2.msh'; %
CONTACT_FACES_1 = [113];
CONTACT_FACES_2 = [113];
DIRICHLET_FACES_1 = [133,142,143,173,182,183];
DIRICHLET_FACES_2 = [133,142,143,173,182,183];
FEM_TYPE = 'FEM_QK(3, 1)'; % Main FEM (has to be Lagrangian)
INTEGRATION = 'IM_HEXAHEDRON(5)'; % Quadrature rule
%%%%% saving parameters %%%%%
ROOTFILENAME = 'static_contact_gears'; % Root of data files
% CONTACT_ALGO = 0 % nodal contact
% CONTACT_ALGO = 1 % integral contact (non-symmetric Alart-Curnier)
% CONTACT_ALGO = 2 % integral contact (symmetric one Alart-Curnier)
% CONTACT_ALGO = 3 % integral contact (non-symmetric Alart-Curnier method with an additional augmentation
% CONTACT_ALGO = 4 % integral contact (new unsymmetric method)
CONTACT_ALGO = 5 % integral large sliding contact
MULT_FEM_TYPE = 'FEM_QK(3, 1)';
|
