linear_elasticity.py
4fa12654f76emaster

linear_elasticity.py
View Options

			# Copyright (C) 2018-2020 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 linear_elasticity_beam_poisson_ratio_engine import \
			LinearElasticityBeamPoissonRatioEngine
			from rod_3d import rod3Dsolver

			def test_elastic_beam():
			solver = LinearElasticityBeamPoissonRatioEngine(n = 10, rho_ = 1e3, g = 3,
			E = 1e6, nu0 = .45, length = 5, verbosity = 0)
			mu = .45
			uh = solver.solve(mu)[0]
			assert np.isclose(solver.norm(uh), 9.33957e-08, rtol = 1e-1)
			assert np.isclose(solver.norm(solver.residual(mu, uh)[0], dual = True),
			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], dual = True),
			1.2210129e-07, rtol = 1e-1)