patch_test_linear_elastic_explicit.cc
No OneTemporary
Actions

Subscribers

None

File Metadata

	/* -------------------------------------------------------------------------- */
	#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();
	}