detection-explicit.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/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)