test_patch_plasticity.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	# -- 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)
	pressure = model['traction'][..., 2]
	pressure[:] = applied_pressure

	solver.solve()
	solver.updateState()

	solution = patch_isotropic_plasticity.solution(applied_pressure)

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