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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
% 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.
% -*- mat-lab -*- (enables emacs matlab mode)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% parameters for program crack %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%% pde parameters : %%%%%
% MU = 77.0; % Lam coefficient.
% LAMBDA = 107.0; % Lam coefficient.
MU = 1.0; % Lam coefficient.
LAMBDA = 1.0; % Lam coefficient.
QUAD = 0;
BIMATERIAL = 1; % 1 : To enable the bimaterial case
if BIMATERIAL
LAMBDA_UP = 10.0;
LAMBDA_DOWN = 1.0;
MU_UP = 10.0;
MU_DOWN = 1.0;
end
REFINEMENT_RADIUS = 0.; %REF: 0.3; % 0 for no refinement
REFINEMENT_PROCESS = 3; % nb of desired refinement steps
REFERENCE_TEST = 0; % 1 to compute a reference solution to be exported into files
ERROR_TO_REF_SOL = 1; % 1 to compute errors with respect to the reference solution
%%%%% discretisation parameters : %%%%%
if (~QUAD)
MESH_TYPE = 'GT_PK(2,1)'; % linear triangles
else
% MESH_TYPE = 'GT_LINEAR_QK(2)';
MESH_TYPE = 'GT_QK(2, 1)';
end;
NX = 10; % space step.
MESH_NOISED = 0; % Set to one if you want to "shake" the mesh
if (~QUAD)
%FEM_TYPE = 'FEM_PK_WITH_CUBIC_BUBBLE(2, 2)';
FEM_TYPE = 'FEM_PK(2, 1)'; % PK element %REF: P3
DATA_FEM_TYPE = 'FEM_PK(2,1)';
INTEGRATION = 'IM_STRUCTURED_COMPOSITE(IM_TRIANGLE(6), 5)';
FEM_TYPE_P = 'FEM_PK(2,1)';
MORTAR_FEM_TYPE = FEM_TYPE;
else
FEM_TYPE = 'FEM_QK(2,1)'; % Q1 fem for quadrangles
DATA_FEM_TYPE = 'FEM_QK(2,1)';
INTEGRATION = 'IM_STRUCTURED_COMPOSITE(IM_GAUSS_PARALLELEPIPED(2, 1), 5)';
FEM_TYPE_P = 'FEM_QK(2,1)';
MORTAR_FEM_TYPE = FEM_TYPE;
end;
FEM_DISC = 'FEM_PK_DISCONTINUOUS(2,3,0.0001)'; % Discontinuous P1 for triangles
MIXED_PRESSURE=0; % Mixed version or not.
DIRICHLET_VERSION = 2;
% integration meth. for sub-simplexe of elements crossed by the level-set
SIMPLEX_INTEGRATION = 'IM_STRUCTURED_COMPOSITE(IM_TRIANGLE(6),3)';
% integration meth. for quasi-polar integration of sub-simplexes adjascent to the level-set
% (comment it to disable quasipolar integration). Should be a
% method defined on a square for 2D, or defined on a prism for 3D.
% SINGULAR_INTEGRATION = 'IM_GAUSS_PARALLELEPIPED(2, 10)';
SINGULAR_INTEGRATION = 'IM_STRUCTURED_COMPOSITE(IM_GAUSS_PARALLELEPIPED(2, 6), 9)';
ADDITIONAL_CRACK = 0;
%Enable the following two lines to use the precalculated solution as enrichement
%GLOBAL_FUNCTION_MF = "bimaterial_crack12.meshfem"
%GLOBAL_FUNCTION_U = "bimaterial_crack12.U"
ENRICHMENT_OPTION = 5; % 0 = Pas d'enrichissement
% 1 = standard XFEM on a fixed zone
% 2 = global functions with mortar junction
% 3 = global functions with cutoff
% 4 = spider fem alone
% 5 = spider fem enrichment
RADIUS_ENR_AREA = 0.2; % For XFEM and the mortar junction
CUTOFF_FUNC = 2; % 0 for the exponential cutoff.
% 1 for a 3rd degree polynomial cutoff
% 2 for a 5th degree polynomial cutoff
CUTOFF = 0.2;
CUTOFF1 = 0.01;
CUTOFF0 = 0.49;
SPIDER_RADIUS = 0.4;
SPIDER_NR = 10; % size of the cartesian mesh in r for spider fem
SPIDER_NTHETA = 10; % size of the cartesian mesh in theta for spider fem
SPIDER_K=1; % order of the spider fem
SPIDER_BIMAT_ENRICHMENT = 1; % 0 for a sqrt spider enrichment
% 1 for sqrt(r)cos(ln(r)) and sqrt(r)sin(ln(r)) spider enrichment
RESIDUAL = 1E-9; % residual for iterative methods if any.
%%%%% saving parameters %%%%%
ROOTFILENAME = 'crack'; % Root of data files.
VTK_EXPORT = 2 % export solution to a .vtk file ?
EXACT_SOL_NUM=1 % 1 for mode I, 2 for mode II, 6 for a combination
% of modes
|
