Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F88670488
test_model_solver.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, Oct 20, 02:17
Size
1 KB
Mime Type
text/x-python
Expires
Tue, Oct 22, 02:17 (1 d, 20 h)
Engine
blob
Format
Raw Data
Handle
21793475
Attached To
rAKA akantu
test_model_solver.py
View Options
#!/usr/bin/env python3
# ------------------------------------------------------------------------------
__author__
=
"Nicolas Richart"
__copyright__
=
"Copyright (C) 2016-2018, EPFL (Ecole Polytechnique Fédérale"
\
" de Lausanne) Laboratory (LSMS - Laboratoire de Simulation"
\
" en Mécanique des Solides)"
__credits__
=
[
"Nicolas Richart"
]
__license__
=
"L-GPLv3"
__maintainer__
=
"Nicolas Richart"
__email__
=
"nicolas.richart@epfl.ch"
# ------------------------------------------------------------------------------
import
numpy
as
np
import
numpy.linalg
as
npl
import
python_fe
as
pfe
sim_u
=
np
.
genfromtxt
(
"disp.csv"
,
delimiter
=
","
,
names
=
True
)
Ne
=
int
(
sim_u
[
'node'
][
-
1
])
F
=
sim_u
[
'force'
]
blocked
=
sim_u
[
'blocked'
]
u
=
sim_u
[
'disp'
]
trusses
=
pfe
.
TrussFE
(
Ne
=
Ne
,
F
=
{
f
:
[
np
.
where
(
F
==
f
)]
for
f
in
np
.
unique
(
F
)},
blocked
=
(
np
.
where
(
blocked
==
1
),
u
[
blocked
==
1
]))
solver
=
pfe
.
StaticSolver
(
trusses
)
solver
.
solveStep
()
upy
=
trusses
.
u
n
=
npl
.
norm
(
upy
-
u
)
/
npl
.
norm
(
upy
)
if
n
>
1e-14
:
raise
ValueError
(
'Something went wrong in the simulation'
)
Event Timeline
Log In to Comment