Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F82352908
detection-explicit.py
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Tue, Sep 10, 22:33
Size
1 KB
Mime Type
text/x-python
Expires
Thu, Sep 12, 22:33 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
20687829
Attached To
rAKA akantu
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
)
Event Timeline
Log In to Comment