Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F84558094
test_hertz.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
Mon, Sep 23, 14:44
Size
1 KB
Mime Type
text/x-python
Expires
Wed, Sep 25, 14:44 (2 d)
Engine
blob
Format
Raw Data
Handle
21006958
Attached To
rTAMAAS tamaas
test_hertz.py
View Options
# -*- coding: utf-8 -*-
# @file
# @section LICENSE
#
# Copyright (©) 2016-2020 EPFL (École Polytechnique Fédérale de Lausanne),
# Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Affero General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Affero General Public License for more details.
#
# You should have received a copy of the GNU Affero General Public License
# along with this program. If not, see <https://www.gnu.org/licenses/>.
from
__future__
import
division
,
print_function
import
numpy
as
np
from
numpy.linalg
import
norm
def
test_contact_area
(
pkr
):
# Testing contact area against Hertz solution for solids of revolution
model
,
hertz
=
pkr
contact_area
=
np
.
mean
(
model
.
getTraction
()
>
0
)
hertz_area
=
np
.
mean
(
hertz
.
pressure
>
0
)
area_error
=
np
.
abs
(
hertz_area
-
contact_area
)
/
hertz_area
assert
area_error
<
1e-2
def
test_pressure
(
pkr
):
model
,
hertz
=
pkr
error
=
norm
(
model
.
getTraction
()
-
hertz
.
pressure
)
/
norm
(
hertz
.
pressure
)
assert
error
<
5e-3
def
test_displacement
(
pkr
):
model
,
hertz
=
pkr
error
=
norm
(
model
.
getDisplacement
()
-
hertz
.
displacement
)
\
/
norm
(
hertz
.
displacement
)
assert
error
<
8e-3
Event Timeline
Log In to Comment