linear_elasticity.py
2aa049ac9d03feature_greedy_coarseni…

linear_elasticity.py
View Options

			# Copyright (C) 2018 by the RROMPy authors
			#
			# This file is part of RROMPy.
			#
			# RROMPy is free software: you can redistribute it and/or modify
			# it under the terms of the GNU Lesser General Public License as published by
			# the Free Software Foundation, either version 3 of the License, or
			# (at your option) any later version.
			#
			# RROMPy 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 Lesser General Public License for more details.
			#
			# You should have received a copy of the GNU Lesser General Public License
			# along with RROMPy. If not, see <http://www.gnu.org/licenses/>.
			#

			import numpy as np
			from rrompy.hfengines.vector_linear_problem import (
			LinearElasticityBeamPoissonRatio)
			from rod_3d import rod3Dsolver

			def test_elastic_beam():
			solver = LinearElasticityBeamPoissonRatio(n = 10, rho_ = 1e3, g = 3,
			E = 1e6, nu0 = .45, length = 5, verbosity = 0)
			mu = .45
			uh = solver.solve(mu)
			assert np.isclose(solver.norm(uh), 58.54349189072907, rtol = 1e-5)
			assert np.isclose(solver.norm(solver.residual(uh, mu)), 8.4545952e-13,
			rtol = 1e-1)

			def test_elastic_rod():
			solver = rod3Dsolver()
			uh = solver.solve(0)
			assert np.isclose(solver.norm(uh), 0.15563476339534466, rtol = 1e-5)
			assert np.isclose(solver.norm(solver.residual(uh, 0)), 5.708389944e-08,
			rtol = 1e-1)