Homec4science
Diffusion akantu (features/34-poisson-solver-with-multiple-dofs) Browse

rAKA/examples/phase_field/phase_field_notch.cc
a4e2455c217efeatures/34-poisson…

Last Change
Raw File

phase_field_notch.cc
View Options

/* -------------------------------------------------------------------------- */
#include "non_linear_solver.hh"
#include "phase_field_model.hh"
#include "solid_mechanics_model.hh"
#include "coupler_solid_phasefield.hh"
/* -------------------------------------------------------------------------- */
#include <iostream>
/* -------------------------------------------------------------------------- */
using namespace akantu;
/* -------------------------------------------------------------------------- */
int main(int argc, char *argv[]){
initialize("material_notch.dat", argc, argv);
// create mesh
Mesh mesh(spatial_dimension);
mesh.read("square_notch.msh");
SolidMechanicsModel model(mesh);
model.initFull(_analysis_method = _static);
auto && mat_selector = std::make_shared<MeshDataMaterialSelector<std::string>>(
"physical_names", model);
model.setMaterialSelector(mat_selector);
PhaseFieldModel phase(mesh);
auto && selector = std::make_shared<MeshDataPhaseFieldSelector<std::string>>(
"physical_names", phase);
phase.setPhaseFieldSelector(selector);
phase.initFull(_analysis_method = _static);
model.applyBC(BC::Dirichlet::FixedValue(0., _y), "bottom");
model.applyBC(BC::Dirichlet::FixedValue(0., _x), "left");
model.setBaseName("phase_notch");
model.addDumpField("stress");
model.addDumpField("grad_u");
model.addDumpFieldVector("displacement");
model.addDumpField("damage");
model.dump();
UInt nbSteps = 1500;
Real increment = 1e-5;
for (UInt s = 0; s < nbSteps; ++s) {
if (s >= 500) {
increment = 1.e-6;
}
model.applyBC(BC::Dirichlet::IncrementValue(increment, _y), "top");
model.solveStep();
model.assembleInternalForces();
computeStrainOnQuadPoints(model, phase, _not_ghost);
phase.solveStep();
computeDamageOnQuadPoints(model, phase, _not_ghost);
phase.dump();
}
finalize();
return EXIT_SUCCESS;
}
c4science · Help