Homec4science
Diffusion tamaas (maxwell_viscoelastic-branch) Browse

tamaas/tests/test_patch_plasticity.py
b9cfa0546a2emaxwell_viscoelastic…

Last Change
Raw File

test_patch_plasticity.py
View Options

# -*- coding: utf-8 -*-
#
# Copyright (©) 2016-20 EPFL (École Polytechnique Fédérale de Lausanne),
# Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
# Copyright (©) 2020-2022 Lucas Frérot
#
# 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/>.
import pytest
from numpy.linalg import norm
import tamaas
from tamaas.nonlinear_solvers import \
DFSANESolver as PySolver, \
DFSANECXXSolver as CXXSolver, \
NewtonKrylovSolver
@pytest.mark.parametrize('solvers', [CXXSolver, NewtonKrylovSolver, PySolver])
def test_patch_plasticity(patch_isotropic_plasticity, solvers):
"Test analyitical solution of 1d plasticity"
tamaas.set_log_level(tamaas.LogLevel.info)
Solver = solvers
model = patch_isotropic_plasticity.model
residual = patch_isotropic_plasticity.residual
applied_pressure = 0.1
solver = Solver(residual)
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 < 3e-15
c4science · Help