issue-template

patch_test_linear_elastic_explicit.cc
View Options

			/* -------------------------------------------------------------------------- */
			#include "patch_test_linear_solid_mechanics_fixture.hh"
			/* -------------------------------------------------------------------------- */

			TYPED_TEST(TestPatchTestSMMLinear, Explicit) {
			std::string filename = "material_check_stress_plane_stress.dat";
			if (this->plane_strain)
			filename = "material_check_stress_plane_strain.dat";

			this->initModel(_explicit_lumped_mass, filename);

			const auto & coordinates = this->mesh->getNodes();
			auto & displacement = this->model->getDisplacement();
			// set the position of all nodes to the static solution
			for (auto && tuple : zip(make_view(coordinates, this->dim),
			make_view(displacement, this->dim))) {
			this->setLinearDOF(std::get<1>(tuple), std::get<0>(tuple));
			}

			for (UInt s = 0; s < 100; ++s) {
			this->model->solveStep();
			}

			auto ekin = this->model->getEnergy("kinetic");
			EXPECT_NEAR(0, ekin, 1e-16);

			this->checkAll();
			}