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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
% Copyright (C) 2017-2017 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.
% This program (test_contact) is used for as a demo of
% Level set contact functionality in Getfem
% It simulates a contact betwen a maste and a slave
% the contact itself occures between the master surface
% and zero contour of level set function, defined on the slave
%
RESIDUAL = 1e-8;
N = 2;
NSTEP = 500;
APPLIED_DISP = -1;
LS_OFFSET = -0.12;
%
%
% Master contact body
DIVxM = 4;
DIVyM = 8;
DIVzM = 1;
xM = 0.3;
yM = 0.92;
zM = 0.1;
LxM = 0.4;
LyM = 0.5;
LzM = 0.1;
APPROX_ORDER_MASTER = 1;
LM_INT_TYPE = 'IM_STRUCTURED_COMPOSITE(IM_GAUSS1D(2),4)';
INT_ORDER_MASTER = 2;
MESH_TYPE_MASTER = 'QK';
LAMBDA_MASTER = 110.0;
MU_MASTER = 70.0;
%
%
% Slave contact body
DIVxS = 20;
DIVyS =
20;
DIVzS = 4;
xS = 0;
yS = 0;
zS = 0;
LxS = 1;
LyS = 1;
LzS = 0.3;
APPROX_ORDER_SLAVE = 1;
INT_ORDER_SLAVE = 2;
MESH_TYPE_SLAVE = 'QK';
LAMBDA_SLAVE = 0;
MU_SLAVE = 7;
|
