Skip to content

Commit

Permalink
Update _storakers.py
Browse files Browse the repository at this point in the history
  • Loading branch information
adtzlr committed Nov 16, 2024
1 parent 2bf245d commit 1000c54
Showing 1 changed file with 5 additions and 3 deletions.
8 changes: 5 additions & 3 deletions src/felupe/constitution/jax/models/hyperelastic/_storakers.py
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,9 @@
along with FElupe. If not, see <http://www.gnu.org/licenses/>.
"""
from functools import wraps
from jax.numpy import array, sum as asum, sqrt

from jax.numpy import array, sqrt
from jax.numpy import sum as asum
from jax.numpy.linalg import eigvalsh

from ....tensortrax.models.hyperelastic import storakers as storakers_docstring
Expand All @@ -26,9 +28,9 @@
def storakers(C, mu, alpha, beta):
λ1, λ2, λ3 = sqrt(eigvalsh(C))
J = λ1 * λ2 * λ3

μ = array(mu)
α = array(alpha)
β = array(beta)

return asum(2 * μ / α**2 * (λ1**α + λ2**α + λ3**α - 3 + (J ** (-α * β) - 1) / β))

0 comments on commit 1000c54

Please sign in to comment.