38-damage…

detection-explicit.py
View Options

			#!/usr/bin/env python3
			""" detection-explicit.py: Python contact detection example"""

			__author__ = "Mohit Pundir"
			__credits__ = [
			"Mohit Pundir <mohit.pundir@epfl.ch>",
			]
			__copyright__ = "Copyright (©) 2018-2021 EPFL (Ecole Polytechnique Fédérale" \
			" de Lausanne) Laboratory (LSMS - Laboratoire de Simulation" \
			" en Mécanique des Solides)"
			__license__ = "LGPLv3"

			import akantu as aka
			import time

			spatial_dimension = 2
			aka.parseInput('detection-explicit.dat')

			mesh = aka.Mesh(spatial_dimension)
			mesh.read('detection-explicit.msh')

			model = aka.ContactMechanicsModel(mesh)

			model.initFull(_analysis_method=aka._explicit_lumped_mass)
			surface_selector = aka.PhysicalSurfaceSelector(mesh)
			model.getContactDetector().setSurfaceSelector(surface_selector)

			model.setBaseName("detection-explicit")
			model.addDumpFieldVector("normals")
			model.addDumpField("gaps")
			model.addDumpField("areas")

			start_time = time.time()
			model.search()
			finish_time = time.time()
			print('Search time = %s seconds', finish_time - start_time)

			model.dump()

			# by default the contact model creates a group named contact_surface
			contact_surface = mesh.getElementGroup("contact_surface")

			normals = model.getNormals()
			gaps = model.getGaps()
			contact_elements = list(model.getContactElements())

			print(normals)
			print(gaps.ravel())
			print(contact_elements)