Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F78161110
stiffness_matrix.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
Sun, Aug 18, 20:20
Size
1 KB
Mime Type
text/x-python
Expires
Tue, Aug 20, 20:20 (2 d)
Engine
blob
Format
Raw Data
Handle
19989540
Attached To
rAKA akantu
stiffness_matrix.py
View Options
#!/usr/bin/env python3
__copyright__
=
(
"Copyright (©) 2019-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
def
getStiffnessMatrix
(
material_file
,
mesh_file
,
traction
):
aka
.
parseInput
(
material_file
)
spatial_dimension
=
2
# --------------------------------------------------------------------------
# Initialization
# --------------------------------------------------------------------------
mesh
=
aka
.
Mesh
(
spatial_dimension
)
mesh
.
read
(
mesh_file
)
model
=
aka
.
SolidMechanicsModel
(
mesh
)
model
.
initFull
(
_analysis_method
=
aka
.
_static
)
model
.
assembleStiffnessMatrix
()
K
=
model
.
getDOFManager
()
.
getMatrix
(
'K'
)
stiff
=
aka
.
AkantuSparseMatrix
(
K
)
.
toarray
()
return
stiff
# --------------------------------------------------------------------------
# main
# --------------------------------------------------------------------------
def
main
():
mesh_file
=
'plate.msh'
material_file
=
'material.dat'
traction
=
1.
mat
=
getStiffnessMatrix
(
material_file
,
mesh_file
,
traction
)
print
(
mat
)
# --------------------------------------------------------------------------
if
__name__
==
"__main__"
:
main
()
Event Timeline
Log In to Comment