declarative

test_patch_plasticity.py
View Options

			# -- coding: utf-8 --
			# @file
			# @section LICENSE
			#
			# Copyright (©) 2016-20 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 numpy.linalg import norm
			from tamaas.nonlinear_solvers import DFSANESolver


			def test_patch_plasticity(patch_isotropic_plasticity):
			model = patch_isotropic_plasticity.model
			residual = patch_isotropic_plasticity.residual

			applied_pressure = 0.1

			solver = DFSANESolver(residual, model)
			solver.tolerance = 1e-15
			pressure = model['traction'][..., 2]
			pressure[:] = applied_pressure

			solver.solve()
			solver.updateState()

			solution, normal = patch_isotropic_plasticity.solution(applied_pressure)

			for key in solution:
			error = norm(model[key] - solution[key]) / normal[key]
			assert error < 1e-15