Diffusion tamaas (features/netcdf-mpi) Browse

netcdf-mpi

rough_contact.py
View Options

			#!/usr/bin/env python3
			# @file
			# @section LICENSE
			#
			# Copyright (©) 2016-2020 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/>.

			import tamaas as tm
			import matplotlib.pyplot as plt

			# Initialize threads and fftw
			tm.initialize()
			tm.set_log_level(tm.LogLevel.info) # Show progression of solver

			# Surface size
			n = 512

			# Surface generator
			sg = tm.SurfaceGeneratorFilter2D()
			sg.setSizes([n, n])
			sg.random_seed = 0

			# Spectrum
			spectrum = tm.Isopowerlaw2D()
			sg.setFilter(spectrum)

			# Parameters
			spectrum.q0 = 16
			spectrum.q1 = 16
			spectrum.q2 = 64
			spectrum.hurst = 0.8

			# Generating surface
			surface = sg.buildSurface()
			# surface /= tm.SurfaceStatistics.computeSpectralRMSSlope(surface)
			surface /= n
			# print(spectrum.rmsSlopes())
			# print(tm.SurfaceStatistics.computeRMSSlope(surface))

			plt.imshow(surface)

			# Creating model
			model = tm.ModelFactory.createModel(tm.model_type_basic_2d, [1., 1.], [n, n])

			# Solver
			solver = tm.PolonskyKeerRey(model, surface, 1e-12,
			tm.PolonskyKeerRey.pressure,
			tm.PolonskyKeerRey.pressure)

			# Solve for target pressure
			p_target = 0.1
			solver.solve(p_target)

			plt.figure()
			plt.imshow(model.getTraction())
			print(model.getTraction().mean())

			# Cleanup threads
			tm.finalize()

			plt.show()

tamaas/examples/rough_contact.py274414d905f3features/netcdf-mpi

rough_contact.pyView Options

tamaas/examples/rough_contact.py
274414d905f3features/netcdf-mpi

rough_contact.py
View Options