materials

test_patch_westergaard.py
View Options

			# -- coding: utf-8 --
			# @file
			# @section LICENSE
			#
			# Copyright (©) 2016-2021 EPFL (École Polytechnique Fédérale de Lausanne),
			# Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
			#
			# This program is free software: you can redistribute it and/or modify
			# it under the terms of the GNU Affero General Public License as published
			# by the Free Software Foundation, either version 3 of the License, or
			# (at your option) any later version.
			#
			# This program is distributed in the hope that it will be useful,
			# but WITHOUT ANY WARRANTY; without even the implied warranty of
			# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
			# GNU Affero General Public License for more details.
			#
			# You should have received a copy of the GNU Affero General Public License
			# along with this program. If not, see <https://www.gnu.org/licenses/>.

			from __future__ import division, print_function

			from numpy.linalg import norm


			def test_patch_westergaard(patch_westergaard):
			model = patch_westergaard.model
			solution = patch_westergaard.solution
			pressure = patch_westergaard.pressure
			model['traction'][:] = pressure[:]
			model.solveNeumann()
			output_displ = model['displacement']

			error = norm(solution - output_displ) / norm(solution)

			assert error < 1e-15 and norm(solution) > 1e-15, \
			"Neumann error = {}".format(error)

			output_displ[:, :] = solution[:, :]
			model.solveDirichlet()
			error = norm(pressure - model['traction']) / norm(pressure)

			assert error < 1e-15, "Dirichlet error = {}".format(error)