37-does-the…

detection-explicit.py
View Options

			#!/usr/bin/env python3
			__copyright__ = (
			"Copyright (©) 2021-2023 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)