28-getelementi…

test_patch_linear_elastic_explicit.py
View Options

			#!/usr/bin/env python3

			# ------------------------------------------------------------------------------
			__author__ = "Guillaume Anciaux"
			__copyright__ = "Copyright (C) 2016-2018, EPFL (Ecole Polytechnique Fédérale" \
			" de Lausanne) Laboratory (LSMS - Laboratoire de Simulation" \
			" en Mécanique des Solides)"
			__credits__ = ["Guillaume Anciaux"]
			__license__ = "L-GPLv3"
			__maintainer__ = "Guillaume Anciaux"
			__email__ = "guillaume.anciaux@epfl.ch"
			# ------------------------------------------------------------------------------
			from patch_test_linear_solid_mechanics_fixture import TestPatchTestSMMLinear
			import akantu
			import sys


			def foo(self):

			filename = "material_check_stress_plane_stress.dat"
			if self.plane_strain:
			filename = "material_check_stress_plane_strain.dat"

			self.initModel(akantu._explicit_lumped_mass, filename)

			coordinates = self.mesh.getNodes()
			displacement = self.model.getDisplacement()
			# set the position of all nodes to the static solution
			self.setLinearDOF(displacement, coordinates)

			for s in range(0, 100):
			self.model.solveStep()

			ekin = self.model.getEnergy("kinetic")
			self.assertAlmostEqual(0, ekin, delta=1e-16)
			self.checkAll()


			def test():
			TestPatchTestSMMLinear.TYPED_TEST(foo, "Explicit")


			if 'pytest' not in sys.modules:
			test()