detection-explicit.py
No OneTemporary
Actions

Subscribers

None

File Metadata

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