v0.14.0
crack_propagation.cpp
Go to the documentation of this file.
1 /* This file is part of MoFEM.
2  * MoFEM is free software: you can redistribute it and/or modify it under
3  * the terms of the GNU Lesser General Public License as published by the
4  * Free Software Foundation, either version 3 of the License, or (at your
5  * option) any later version.
6  *
7  * MoFEM is distributed in the hope that it will be useful, but WITHOUT
8  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
9  * FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
10  * License for more details.
11  *
12  * You should have received a copy of the GNU Lesser General Public
13  * License along with MoFEM. If not, see <http://www.gnu.org/licenses/>. */
14 
15 #ifdef WITH_TETGEN
16 
17 #include <tetgen.h>
18 #ifdef REAL
19 #undef REAL
20 #endif
21 
22 #endif
23 
24 #include <MoFEM.hpp>
25 using namespace MoFEM;
26 #include <BasicFiniteElements.hpp>
27 #include <Mortar.hpp>
28 #include <NeoHookean.hpp>
29 #include <Hooke.hpp>
30 
31 #define BOOST_UBLAS_NDEBUG 1
32 
33 #include <CrackFrontElement.hpp>
34 #include <ComplexConstArea.hpp>
35 #include <MWLS.hpp>
36 #include <GriffithForceElement.hpp>
37 #include <VolumeLengthQuality.hpp>
38 #include <CrackPropagation.hpp>
39 #include <CPSolvers.hpp>
40 #include <CPMeshCut.hpp>
41 #include <AnalyticalFun.hpp>
42 
43 static char help[] = "Calculate crack release energy and crack propagation"
44  "\n\n";
45 
46 using namespace FractureMechanics;
47 
48 bool CrackPropagation::debug = false;
49 
50 bool CrackPropagation::parallelReadAndBroadcast =
51  false; ///< That is unstable development, for
52  ///< some meshses (propably generated by
53  ///< cubit) this is not working. Error can
54  ///< be attributed to bug in MOAB.
55 
56 int main(int argc, char *argv[]) {
57 
58  const char param_file[] = "param_file.petsc";
59  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
60 
61  try {
62 
63  {
64  PetscBool flg = PETSC_FALSE;
65  CHKERR PetscOptionsGetBool(PETSC_NULL, "", "-debug", &flg, PETSC_NULL);
66  if (flg == PETSC_TRUE)
67  CrackPropagation::debug = true;
68  }
69 
70  {
71  PetscBool flg = PETSC_FALSE;
72  CHKERR PetscOptionsGetBool(PETSC_NULL, "", "-read_and_brodcast", &flg,
73  PETSC_NULL);
74  if (flg == PETSC_TRUE)
75  CrackPropagation::parallelReadAndBroadcast = true;
76  }
77 
78  PetscBool flg = PETSC_FALSE;
79  char mesh_file_name[255];
80  CHKERR PetscOptionsGetString(PETSC_NULL, "", "-my_file", mesh_file_name,
81  255, &flg);
82 
83  const char *tan_list[] = {"tangent_internal_stress", "density_tangent",
84  "griffith_tangent"};
85  int choice_test = MWLS_STRESS_TAN;
86  PetscBool isTangentTest;
87  CHKERR PetscOptionsGetEList(PETSC_NULL, NULL, "-testing_mwls", tan_list,
88  LASTTAN, &choice_test, &isTangentTest);
89  moab::Core mb_instance;
90  moab::Interface &moab = mb_instance;
91 
92  ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);
93  auto moab_comm_wrap =
94  boost::make_shared<WrapMPIComm>(PETSC_COMM_WORLD, false);
95  if (pcomm == NULL)
96  pcomm = new ParallelComm(&moab, moab_comm_wrap->get_comm());
97 
98  if (pcomm->size() == 1 || !CrackPropagation::parallelReadAndBroadcast) {
99  if (flg) {
100  const char *option = "";
101  CHKERR moab.load_file(mesh_file_name, 0, option);
102  }
103  } else {
104  if (pcomm->rank() == 0) {
105  // Read mesh file
106  if (flg) {
107  moab::Core mb_instance_read;
108  moab::Interface &moab_read = mb_instance_read;
109  ParallelComm *pcomm =
110  ParallelComm::get_pcomm(&moab_read, MYPCOMM_INDEX);
111  if (pcomm == NULL)
112  pcomm = new ParallelComm(&moab_read, moab_comm_wrap->get_comm());
113 
114  const char *option = "";
115  CHKERR moab_read.load_file(mesh_file_name, 0, option);
116  Range ents;
117  CHKERR moab_read.get_entities_by_handle(0, ents, false);
118  CHKERR moab_read.get_entities_by_type(0, MBENTITYSET, ents, false);
119  CHKERR FractureMechanics::broadcast_entities(moab, moab_read, 0, ents,
120  false, true);
121  }
122  } else {
123 
124  Range ents;
125  CHKERR FractureMechanics::broadcast_entities(moab, moab, 0, ents, false,
126  true);
127  }
128  CHKERR FractureMechanics::clean_pcomms(moab, moab_comm_wrap);
129  }
130 
131  // Create mofem database
132  MeshsetsManager::brodcastMeshsets = false;
133  MoFEM::Core core(moab, PETSC_COMM_WORLD, VERBOSE);
134  MoFEM::Interface &m_field = core;
135 
136  auto core_log = logging::core::get();
137  core_log->add_sink(
138  LogManager::createSink(LogManager::getStrmWorld(), "CP"));
139  LogManager::setLog("CP");
140  MOFEM_LOG_TAG("CP", "cp");
141  MOFEM_LOG_C("CP", Sev::inform, "### Input parameter: -my_file %s",
142  mesh_file_name);
143 
144  // Create data structure for crack propagation
145  CrackPropagation cp(m_field, 3, 2);
146  cp.moabCommWorld = moab_comm_wrap;
147 
148  CHKERR MoFEM::Interface::getFileVersion(moab, cp.fileVersion);
149  MOFEM_LOG("CP", Sev::inform)
150  << "File version " << cp.fileVersion.strVersion();
151 
152  // Configure meshes
153  CHKERR cp.getOptions();
154 
155  // Register MOFEM discrete data managers
156  {
157  DMType dm_name = "MOFEM";
158  CHKERR DMRegister_MoFEM(dm_name);
159  }
160 
161  auto add_bits_tets = [&]() {
162  MoFEMFunctionBegin;
163  CHKERR m_field.getInterface<BitRefManager>()
164  ->addToDatabaseBitRefLevelByType(MBTET, BitRefLevel().set(),
165  BitRefLevel().set());
166  MoFEMFunctionReturn(0);
167  };
168 
169  CHKERR add_bits_tets();
170 
171  if (CrackPropagation::debug) {
172  std::string fn =
173  "mesh_after_cut_and_broadcast_" +
174  boost::lexical_cast<std::string>(m_field.get_comm_rank()) + ".vtk";
175  CHKERR moab.write_file(fn.c_str(), "VTK", "");
176  }
177 
178  if (string(mesh_file_name).compare(0, 8, "restart_") == 0) {
179  auto first_number_string = [](std::string const &str) {
180  std::size_t const n = str.find_first_of("0123456789");
181  if (n != std::string::npos) {
182  std::size_t const m = str.find_first_not_of("0123456789", n);
183  return str.substr(n, m != std::string::npos ? m - n : m);
184  }
185  return std::string();
186  };
187  std::string str_number = first_number_string(std::string(mesh_file_name));
188 
189  if (!str_number.empty())
190  cp.startStep = atoi(str_number.c_str()) + 1;
191  }
192 
193  // create bit levels, fields and finite elements structures in database
194  if (cp.doCutMesh) {
195 
196  if (string(mesh_file_name).compare(0, 8, "restart_") == 0) {
197 
198  Range ents;
199  CHKERR m_field.getInterface<BitRefManager>()
200  ->getAllEntitiesNotInDatabase(ents);
201  Range meshsets;
202  CHKERR m_field.get_moab().get_entities_by_type(0, MBENTITYSET, meshsets,
203  false);
204  for (auto m : meshsets)
205  CHKERR m_field.get_moab().remove_entities(m, ents);
206  CHKERR m_field.get_moab().delete_entities(ents);
207 
208  string cutting_surface_name =
209  "cutting_surface_" +
210  boost::lexical_cast<std::string>(cp.startStep) + ".vtk";
211 
212  if (m_field.get_comm_rank() == 0)
213  CHKERR cp.getInterface<CPMeshCut>()->rebuildCrackSurface(
214  cp.crackAccelerationFactor, cutting_surface_name, QUIET,
215  CrackPropagation::debug);
216  else
217  CHKERR cp.getInterface<CPMeshCut>()->rebuildCrackSurface(
218  cp.crackAccelerationFactor, "", QUIET, CrackPropagation::debug);
219 
220  CHKERR cp.getInterface<CPMeshCut>()->refineMesh(
221  nullptr, false, QUIET, CrackPropagation::debug);
222 
223  } else {
224 
225  string cutting_surface_name =
226  "cutting_surface_" +
227  boost::lexical_cast<std::string>(cp.startStep) + ".vtk";
228 
229  if (!cp.getInterface<CPMeshCut>()->getCutSurfMeshName().empty())
230  CHKERR cp.getInterface<CPMeshCut>()->setCutSurfaceFromFile();
231  CHKERR cp.getInterface<CPMeshCut>()->copySurface(cutting_surface_name);
232  Range vol;
233  CHKERR moab.get_entities_by_type(0, MBTET, vol, false);
234  CHKERR cp.getInterface<CPMeshCut>()->refineMesh(
235  &vol, true, QUIET, CrackPropagation::debug);
236  }
237 
238  CHKERR cp.getInterface<CPMeshCut>()->cutRefineAndSplit(
239  QUIET, CrackPropagation::debug);
240 
241  } else {
242 
243  Range vol;
244  CHKERR moab.get_entities_by_type(0, MBTET, vol, false);
245  BitRefLevel bit0 = cp.mapBitLevel["mesh_cut"];
246  CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevel(
247  vol, bit0.set(BITREFLEVEL_SIZE - 1), true, CrackPropagation::debug);
248 
249  // intact mesh will be saved in "spatial_domain" mapBitLevel to solve
250  // the elastic problem in the absence of the crack (restart files won't
251  // work)
252  if (!cp.doCutMesh && !cp.propagateCrack && cp.doElasticWithoutCrack) {
253  BitRefLevel bit1 = cp.mapBitLevel["spatial_domain"];
254  CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevel(
255  vol, bit1.set(BITREFLEVEL_SIZE - 1), true, CrackPropagation::debug);
256  } else {
257  CHKERR cp.getInterface<CPMeshCut>()->refineAndSplit(
258  QUIET, CrackPropagation::debug);
259  }
260  }
261 
262  if (string(mesh_file_name).compare(0, 8, "restart_") == 0) {
263  // restart code from file
264  CHKERR m_field.build_field("LAMBDA_ARC_LENGTH");
265  CHKERR m_field.build_finite_elements("ARC_LENGTH");
266  }
267 
268  std::vector<int> surface_ids;
269  surface_ids.push_back(cp.getInterface<CPMeshCut>()->getSkinOfTheBodyId());
270  std::vector<std::string> fe_surf_proj_list;
271  fe_surf_proj_list.push_back("SURFACE");
272 
273  CHKERR cp.createProblemDataStructures(surface_ids, QUIET,
274  CrackPropagation::debug);
275 
276  if (isTangentTest) {
277 
278  MOFEM_LOG("CP", Sev::verbose) << "Testing tangent variant";
279 
280  // set inital coordinates
281  CHKERR cp.setMaterialPositionFromCoords();
282  CHKERR cp.setSpatialPositionFromCoords();
283  if (cp.addSingularity == PETSC_TRUE) {
284  CHKERR cp.setSingularDofs("SPATIAL_POSITION");
285  }
286 
287  SmartPetscObj<DM> dm_elastic;
288  SmartPetscObj<DM> dm_eigen_elastic;
289  SmartPetscObj<DM> dm_material;
290  SmartPetscObj<DM> dm_crack_propagation;
291  SmartPetscObj<DM> dm_material_forces;
292  SmartPetscObj<DM> dm_surface_projection;
293  SmartPetscObj<DM> dm_crack_srf_area;
294  CHKERR cp.createDMs(dm_elastic, dm_eigen_elastic, dm_material,
295  dm_crack_propagation, dm_material_forces,
296  dm_surface_projection, dm_crack_srf_area, surface_ids,
297  fe_surf_proj_list);
298 
299  // set finite element instances and user data operators on instances
300  CHKERR cp.assembleElasticDM(VERBOSE, false);
301  CHKERR cp.assembleMaterialForcesDM(PETSC_NULL);
302  CHKERR cp.assembleSmootherForcesDM(PETSC_NULL, surface_ids, VERBOSE,
303  false);
304  CHKERR cp.assembleCouplingForcesDM(PETSC_NULL, VERBOSE, false);
305 
306  // set inital coordinates
307  CHKERR cp.setMaterialPositionFromCoords();
308  CHKERR cp.setSpatialPositionFromCoords();
309  if (cp.addSingularity == PETSC_TRUE)
310  CHKERR cp.setSingularDofs("SPATIAL_POSITION");
311 
312  switch (choice_test) {
313  case MWLS_STRESS_TAN:
314  CHKERR cp.testJacobians(cp.mapBitLevel["material_domain"],
315  BitRefLevel().set(), MWLS_STRESS_TAN);
316  break;
317  case MWLS_DENSITY_TAN:
318  CHKERR cp.testJacobians(cp.mapBitLevel["material_domain"],
319  BitRefLevel().set(), MWLS_DENSITY_TAN);
320  break;
321  case MWLS_GRIFFITH_TAN:
322  CHKERR cp.testJacobians(cp.mapBitLevel["material_domain"],
323  BitRefLevel().set(), MWLS_GRIFFITH_TAN);
324  break;
325  default:
326  break;
327  }
328 
329  } else if (cp.propagateCrack == PETSC_TRUE) {
330 
331  MOFEM_LOG("CP", Sev::verbose) << "Crack propagation variant is used";
332 
333  // set inital coordinates
334  CHKERR cp.setMaterialPositionFromCoords();
335  CHKERR cp.setSpatialPositionFromCoords();
336  if (cp.addSingularity == PETSC_TRUE)
337  CHKERR cp.setSingularDofs("SPATIAL_POSITION");
338 
339  SmartPetscObj<DM> dm_elastic;
340  SmartPetscObj<DM> dm_eigen_elastic;
341  SmartPetscObj<DM> dm_material;
342  SmartPetscObj<DM> dm_crack_propagation;
343  SmartPetscObj<DM> dm_material_forces;
344  SmartPetscObj<DM> dm_surface_projection;
345  SmartPetscObj<DM> dm_crack_srf_area;
346  CHKERR cp.createDMs(dm_elastic, dm_eigen_elastic, dm_material,
347  dm_crack_propagation, dm_material_forces,
348  dm_surface_projection, dm_crack_srf_area, surface_ids,
349  fe_surf_proj_list);
350 
351  CHKERR cp.setSpatialPositionFromCoords();
352  if (cp.addSingularity == PETSC_TRUE) {
353  CHKERR cp.setSingularDofs("SPATIAL_POSITION");
354  }
355  if (string(mesh_file_name).compare(0, 8, "restart_") == 0)
356  cp.getLoadScale() = cp.getArcCtx()->getFieldData();
357 
358  // Solve eigen problem
359  const auto load_factor = std::abs(cp.getLoadScale());
360  MOFEM_LOG("CPSolverSync", Sev::noisy) << "Lambda factor " << load_factor;
361  MOFEM_LOG_SYNCHRONISE(m_field.get_comm());
362  if (cp.solveEigenStressProblem) {
363  MOFEM_LOG("CP", Sev::verbose) << "Solve Eigen ELASTIC problem";
364  CHKERR cp.getInterface<CPSolvers>()->solveElastic(
365  dm_eigen_elastic, cp.getEigenArcCtx()->x0, 1);
366  CHKERR cp.postProcessDM(dm_eigen_elastic, -2, "ELASTIC", true);
367  }
368  MOFEM_LOG("CP", Sev::noisy) << "Solve ELASTIC problem";
369  cp.getLoadScale() = load_factor;
370  MOFEM_LOG("CPSolverSync", Sev::noisy)
371  << "Reset lambda factor " << cp.getLoadScale();
372  MOFEM_LOG_SYNCHRONISE(m_field.get_comm());
373  CHKERR cp.getInterface<CPSolvers>()->solveElastic(
374  dm_elastic, cp.getArcCtx()->x0, load_factor);
375  CHKERR cp.postProcessDM(dm_elastic, -1, "ELASTIC", true);
376 
377  // set finite element instances and user data operators on instances
378  CHKERR cp.assembleElasticDM(VERBOSE, false);
379  CHKERR cp.assembleMaterialForcesDM(PETSC_NULL);
380  CHKERR cp.assembleSmootherForcesDM(PETSC_NULL, surface_ids, VERBOSE,
381  false);
382  CHKERR cp.assembleCouplingForcesDM(PETSC_NULL, VERBOSE, false);
383 
384  CHKERR cp.calculateElasticEnergy(dm_elastic);
385  CHKERR cp.getInterface<CPSolvers>()->calculateGc(
386  dm_material_forces, dm_surface_projection, dm_crack_srf_area,
387  surface_ids);
388 
389  if (cp.doCutMesh == PETSC_TRUE) {
390  MOFEM_LOG("CP", Sev::verbose) << "Propagate crack and cut mesh";
391  CHKERR cp.getInterface<CPSolvers>()->solveCutMesh(
392  dm_crack_propagation, dm_elastic, dm_eigen_elastic, dm_material,
393  dm_material_forces, dm_surface_projection, dm_crack_srf_area,
394  surface_ids, CrackPropagation::debug);
395  } else {
396  MOFEM_LOG("CP", Sev::verbose) << "Propagate crack (not cut mesh)";
397  CHKERR cp.getInterface<CPSolvers>()->solvePropagation(
398  dm_crack_propagation, dm_elastic, dm_material, dm_material_forces,
399  dm_surface_projection, dm_crack_srf_area, surface_ids);
400  }
401 
402  } else {
403 
404  MOFEM_LOG("CP", Sev::verbose) << "Do not propagate crack";
405 
406  // set inital coordinates
407  CHKERR cp.setMaterialPositionFromCoords();
408  CHKERR cp.setSpatialPositionFromCoords();
409  if (cp.addSingularity == PETSC_TRUE)
410  CHKERR cp.setSingularDofs("SPATIAL_POSITION");
411 
412  SmartPetscObj<DM> dm_elastic;
413  SmartPetscObj<DM> dm_eigen_elastic;
414  SmartPetscObj<DM> dm_material;
415  SmartPetscObj<DM> dm_crack_propagation;
416  SmartPetscObj<DM> dm_material_forces;
417  SmartPetscObj<DM> dm_surface_projection;
418  SmartPetscObj<DM> dm_crack_srf_area;
419  CHKERR cp.createDMs(dm_elastic, dm_eigen_elastic, dm_material,
420  dm_crack_propagation, dm_material_forces,
421  dm_surface_projection, dm_crack_srf_area, surface_ids,
422  fe_surf_proj_list);
423 
424  SmartPetscObj<DM> dm_bc;
425  if (m_field.getInterface<MeshsetsManager>()->checkMeshset(
426  "ANALITICAL_DISP"))
427  CHKERR cp.createBcDM(dm_bc, "BC_PROBLEM",
428  cp.mapBitLevel["spatial_domain"],
429  BitRefLevel().set());
430 
431  // Solve eigen problem
432  if (cp.solveEigenStressProblem) {
433  MOFEM_LOG("CP", Sev::verbose) << "Solve Eigen ELASTIC problem";
434  CHKERR cp.getInterface<CPSolvers>()->solveElastic(
435  dm_eigen_elastic, cp.getEigenArcCtx()->x0, 1);
436  CHKERR cp.postProcessDM(dm_eigen_elastic, -2, "ELASTIC", true);
437  }
438  // set finite element instances and user data operators on instances
439  MOFEM_LOG("CP", Sev::noisy) << "Solve ELASTIC problem";
440  CHKERR cp.getInterface<CPSolvers>()->solveElastic(
441  dm_elastic, cp.getArcCtx()->x0, cp.getLoadScale());
442 
443  // set finite element instances and user data operators on instances
444  CHKERR cp.assembleElasticDM(VERBOSE, false);
445  CHKERR cp.assembleMaterialForcesDM(PETSC_NULL);
446  if (cp.doCutMesh)
447  CHKERR cp.assembleSmootherForcesDM(PETSC_NULL, surface_ids, VERBOSE,
448  false);
449 
450  CHKERR cp.getInterface<CPSolvers>()->calculateGc(
451  dm_material_forces, dm_surface_projection, dm_crack_srf_area,
452  surface_ids);
453  CHKERR cp.calculateElasticEnergy(dm_elastic);
454 
455  // post-processing
456  CHKERR cp.postProcessDM(dm_elastic, 0, "ELASTIC", true);
457  // set vector values from field data
458  CHKERR cp.savePositionsOnCrackFrontDM(dm_elastic, PETSC_NULL, 2, false);
459  if (!cp.doElasticWithoutCrack)
460  CHKERR cp.writeCrackFont(cp.mapBitLevel["spatial_domain"], 0);
461 
462  // Run tests
463  CHKERR cp.tetsingReleaseEnergyCalculation();
464  }
465  }
466  CATCH_ERRORS;
467 
468  MoFEM::Core::Finalize();
469  return 0;
470 }
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface reference to pointer of interface.
Definition: UnknownInterface.hpp:93
FractureMechanics::CrackPropagation::writeCrackFont
MoFEMErrorCode writeCrackFont(const BitRefLevel &bit, const int step)
Definition: CrackPropagation.cpp:8093
FractureMechanics::CrackPropagation::tetsingReleaseEnergyCalculation
MoFEMErrorCode tetsingReleaseEnergyCalculation()
This is run with ctest.
Definition: CrackPropagation.cpp:929
FractureMechanics::CrackPropagation::createBcDM
MoFEMErrorCode createBcDM(SmartPetscObj< DM > &dm, const std::string prb_name, const BitRefLevel bit, const BitRefLevel mask=BitRefLevel().set())
Create problem to calculate boundary conditions.
Definition: CrackPropagation.cpp:3625
MYPCOMM_INDEX
#define MYPCOMM_INDEX
default communicator number PCOMM
Definition: definitions.h:228
MoFEM::CoreTmp< 0 >
Core (interface) class.
Definition: Core.hpp:82
MoFEM::Version::strVersion
std::string strVersion()
Definition: UnknownInterface.hpp:24
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0
FractureMechanics::CrackPropagation::mapBitLevel
std::map< std::string, BitRefLevel > mapBitLevel
Definition: CrackPropagation.hpp:180
MWLS.hpp
Mortar.hpp
FractureMechanics::CrackPropagation::setSingularDofs
MoFEMErrorCode setSingularDofs(const string field_name, const int verb=QUIET)
set singular dofs (on edges adjacent to crack front) from geometry
Definition: CrackPropagation.cpp:9334
FractureMechanics::CrackPropagation::solveEigenStressProblem
PetscBool solveEigenStressProblem
Solve eigen problem.
Definition: CrackPropagation.hpp:163
FractureMechanics::CrackPropagation::getArcCtx
boost::shared_ptr< ArcLengthCtx > & getArcCtx()
Definition: CrackPropagation.hpp:907
FractureMechanics::broadcast_entities
MoFEMErrorCode broadcast_entities(moab::Interface &moab, moab::Interface &moab_tmp, const int from_proc, Range &entities, const bool adjacencies, const bool tags)
Definition: CrackPropagation.cpp:149
FractureMechanics::CrackPropagation::assembleCouplingForcesDM
MoFEMErrorCode assembleCouplingForcesDM(DM dm, const int verb=QUIET, const bool debug=false)
assemble coupling element instances
Definition: CrackPropagation.cpp:6178
MoFEM::CoreInterface::get_comm_rank
virtual int get_comm_rank() const =0
MoFEM::MeshsetsManager::brodcastMeshsets
static bool brodcastMeshsets
Definition: MeshsetsManager.hpp:881
MoFEM.hpp
MoFEM::UnknownInterface::getFileVersion
static MoFEMErrorCode getFileVersion(moab::Interface &moab, Version &version)
Get database major version.
Definition: UnknownInterface.cpp:16
FractureMechanics::MWLS_DENSITY_TAN
@ MWLS_DENSITY_TAN
Definition: CrackPropagation.hpp:56
BasicFiniteElements.hpp
MoFEM::CoreTmp< 0 >::Finalize
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
Definition: Core.cpp:112
FractureMechanics::CrackPropagation::getEigenArcCtx
boost::shared_ptr< ArcLengthCtx > & getEigenArcCtx()
Definition: CrackPropagation.hpp:908
FractureMechanics::CrackPropagation::setSpatialPositionFromCoords
MoFEMErrorCode setSpatialPositionFromCoords()
set spatial field from nodes
Definition: CrackPropagation.cpp:9330
MoFEM::LogManager::createSink
static boost::shared_ptr< SinkType > createSink(boost::shared_ptr< std::ostream > stream_ptr, std::string comm_filter)
Create a sink object.
Definition: LogManager.cpp:298
FractureMechanics::CPMeshCut
Definition: CPMeshCut.hpp:24
MoFEM::DeprecatedCoreInterface
Deprecated interface functions.
Definition: DeprecatedCoreInterface.hpp:16
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:2010
CPSolvers.hpp
Solvers for crack propagation.
CPMeshCut.hpp
FractureMechanics::CPMeshCut::getCutSurfMeshName
const std::string & getCutSurfMeshName() const
Definition: CPMeshCut.hpp:157
FractureMechanics::CrackPropagation::setMaterialPositionFromCoords
MoFEMErrorCode setMaterialPositionFromCoords()
set material field from nodes
Definition: CrackPropagation.cpp:9326
VERBOSE
@ VERBOSE
Definition: definitions.h:222
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:548
MoFEM::CoreInterface::get_moab
virtual moab::Interface & get_moab()=0
MoFEM
implementation of Data Operators for Forces and Sources
Definition: Common.hpp:10
MoFEM::CoreInterface::build_finite_elements
virtual MoFEMErrorCode build_finite_elements(int verb=DEFAULT_VERBOSITY)=0
Build finite elements.
MOFEM_LOG_C
#define MOFEM_LOG_C(channel, severity, format,...)
Definition: LogManager.hpp:311
Hooke.hpp
Implementation of linear elastic material.
MOFEM_LOG_SYNCHRONISE
#define MOFEM_LOG_SYNCHRONISE(comm)
Synchronise "SYNC" channel.
Definition: LogManager.hpp:345
FractureMechanics::CrackPropagation::assembleSmootherForcesDM
MoFEMErrorCode assembleSmootherForcesDM(DM dm, const std::vector< int > ids, const int verb=QUIET, const bool debug=false)
create smoothing element instance
Definition: CrackPropagation.cpp:5975
MoFEM::LogManager::getStrmWorld
static boost::shared_ptr< std::ostream > getStrmWorld()
Get the strm world object.
Definition: LogManager.cpp:344
FractureMechanics::MWLS_GRIFFITH_TAN
@ MWLS_GRIFFITH_TAN
Definition: CrackPropagation.hpp:57
MoFEM::DMRegister_MoFEM
PetscErrorCode DMRegister_MoFEM(const char sname[])
Register MoFEM problem.
Definition: DMMoFEM.cpp:43
FractureMechanics::CrackPropagation::doCutMesh
PetscBool doCutMesh
Definition: CrackPropagation.hpp:1327
main
int main(int argc, char *argv[])
Definition: crack_propagation.cpp:56
help
static char help[]
Definition: crack_propagation.cpp:43
FractureMechanics::CrackPropagation::fileVersion
Version fileVersion
Definition: CrackPropagation.hpp:86
FractureMechanics::LASTTAN
@ LASTTAN
Definition: CrackPropagation.hpp:58
FractureMechanics::clean_pcomms
MoFEMErrorCode clean_pcomms(moab::Interface &moab, boost::shared_ptr< WrapMPIComm > moab_comm_wrap)
MOFEM_LOG_TAG
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
Definition: LogManager.hpp:339
FractureMechanics::CrackPropagation::propagateCrack
PetscBool propagateCrack
If true crack propagation is calculated.
Definition: CrackPropagation.hpp:119
CrackFrontElement.hpp
FractureMechanics::CrackPropagation::assembleElasticDM
MoFEMErrorCode assembleElasticDM(const std::string mwls_stress_tag_name, const int verb=QUIET, const bool debug=false)
create elastic finite element instance for spatial assembly
Definition: CrackPropagation.cpp:3795
FractureMechanics::CrackPropagation::calculateElasticEnergy
MoFEMErrorCode calculateElasticEnergy(DM dm, const std::string msg="")
assemble elastic part of matrix, by running elastic finite element instance
Definition: CrackPropagation.cpp:7280
FractureMechanics::CrackPropagation::crackAccelerationFactor
double crackAccelerationFactor
Definition: CrackPropagation.hpp:1328
debug
static const bool debug
Definition: dm_create_subdm.cpp:12
FractureMechanics::CrackPropagation::postProcessDM
MoFEMErrorCode postProcessDM(DM dm, const int step, const std::string fe_name, const bool approx_internal_stress)
post-process results for elastic solution
Definition: CrackPropagation.cpp:9010
FractureMechanics::CPMeshCut::getSkinOfTheBodyId
int getSkinOfTheBodyId() const
Definition: CPMeshCut.hpp:167
FractureMechanics::MWLS_STRESS_TAN
@ MWLS_STRESS_TAN
Definition: CrackPropagation.hpp:55
mesh_file_name
char mesh_file_name[255]
Definition: mesh_smoothing.cpp:22
convert.n
n
Definition: convert.py:82
FractureMechanics::CrackPropagation::savePositionsOnCrackFrontDM
MoFEMErrorCode savePositionsOnCrackFrontDM(DM dm, Vec q, const int verb=QUIET, const bool debug=false)
Definition: CrackPropagation.cpp:8064
Range
MoFEM::CoreTmp< 0 >::Initialize
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
Definition: Core.cpp:72
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
FractureMechanics::CPSolvers
Definition: CPSolvers.hpp:28
FractureMechanics::CrackPropagation::moabCommWorld
boost::shared_ptr< WrapMPIComm > moabCommWorld
Definition: CrackPropagation.hpp:104
CATCH_ERRORS
#define CATCH_ERRORS
Catch errors.
Definition: definitions.h:385
AnalyticalFun.hpp
MoFEM::Core
CoreTmp< 0 > Core
Definition: Core.hpp:1148
BITREFLEVEL_SIZE
#define BITREFLEVEL_SIZE
max number of refinements
Definition: definitions.h:232
FractureMechanics::CrackPropagation::testJacobians
MoFEMErrorCode testJacobians(const BitRefLevel bit, const BitRefLevel mask, tangent_tests test)
test LHS Jacobians
Definition: CrackPropagation.cpp:5279
FractureMechanics::CrackPropagation
Definition: CrackPropagation.hpp:77
FractureMechanics::CrackPropagation::getLoadScale
double & getLoadScale()
Definition: CrackPropagation.hpp:210
FractureMechanics::CrackPropagation::assembleMaterialForcesDM
MoFEMErrorCode assembleMaterialForcesDM(DM dm, const int verb=QUIET, const bool debug=false)
create material element instance
Definition: CrackPropagation.cpp:5531
ComplexConstArea.hpp
MoFEM::PetscOptionsGetString
PetscErrorCode PetscOptionsGetString(PetscOptions *, const char pre[], const char name[], char str[], size_t size, PetscBool *set)
Definition: DeprecatedPetsc.hpp:172
MoFEM::MeshsetsManager
Interface for managing meshsets containing materials and boundary conditions.
Definition: MeshsetsManager.hpp:104
MoFEM::PetscOptionsGetEList
PetscErrorCode PetscOptionsGetEList(PetscOptions *, const char pre[], const char name[], const char *const *list, PetscInt next, PetscInt *value, PetscBool *set)
Definition: DeprecatedPetsc.hpp:203
m
FTensor::Index< 'm', 3 > m
Definition: shallow_wave.cpp:80
MoFEM::BitRefManager
Managing BitRefLevels.
Definition: BitRefManager.hpp:21
CrackPropagation.hpp
Main class for crack propagation.
MoFEM::Types::BitRefLevel
std::bitset< BITREFLEVEL_SIZE > BitRefLevel
Bit structure attached to each entity identifying to what mesh entity is attached.
Definition: Types.hpp:40
MoFEM::CoreInterface::build_field
virtual MoFEMErrorCode build_field(const std::string field_name, int verb=DEFAULT_VERBOSITY)=0
build field by name
QUIET
@ QUIET
Definition: definitions.h:221
MoFEM::LogManager::setLog
static LoggerType & setLog(const std::string channel)
Set ans resset chanel logger.
Definition: LogManager.cpp:389
MoFEM::SmartPetscObj< DM >
NeoHookean.hpp
Implementation of Neo-Hookean elastic material.
FractureMechanics::CrackPropagation::createDMs
MoFEMErrorCode createDMs(SmartPetscObj< DM > &dm_elastic, SmartPetscObj< DM > &dm_eigen_elastic, SmartPetscObj< DM > &dm_material, SmartPetscObj< DM > &dm_crack_propagation, SmartPetscObj< DM > &dm_material_forces, SmartPetscObj< DM > &dm_surface_projection, SmartPetscObj< DM > &dm_crack_srf_area, std::vector< int > surface_ids, std::vector< std::string > fe_surf_proj_list)
Crate DMs for all problems and sub problems.
Definition: CrackPropagation.cpp:9571
MoFEM::MeshsetsManager::checkMeshset
bool checkMeshset(const int ms_id, const CubitBCType cubit_bc_type) const
check for CUBIT Id and CUBIT type
Definition: MeshsetsManager.cpp:360
GriffithForceElement.hpp
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
FractureMechanics::CrackPropagation::getOptions
MoFEMErrorCode getOptions()
Get options form command line.
Definition: CrackPropagation.cpp:588
FractureMechanics::CrackPropagation::createProblemDataStructures
MoFEMErrorCode createProblemDataStructures(const std::vector< int > surface_ids, const int verb=QUIET, const bool debug=false)
Construct problem data structures.
Definition: CrackPropagation.cpp:9525
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
FractureMechanics::CrackPropagation::doElasticWithoutCrack
PetscBool doElasticWithoutCrack
Definition: CrackPropagation.hpp:1329
FractureMechanics
Definition: AnalyticalFun.hpp:15
FractureMechanics::CrackPropagation::addSingularity
PetscBool addSingularity
Definition: CrackPropagation.hpp:122
FractureMechanics::CrackPropagation::startStep
int startStep
Definition: CrackPropagation.hpp:168
MoFEM::PetscOptionsGetBool
PetscErrorCode PetscOptionsGetBool(PetscOptions *, const char pre[], const char name[], PetscBool *bval, PetscBool *set)
Definition: DeprecatedPetsc.hpp:182
Generated by Doxygen 1.8.17 and hosted at University of Glasgow