mirror of
https://github.com/JuliaFEM/JuliaFEM.jl.git
synced 2026-09-17 17:22:10 +00:00
Jukka Aho
57 lines
2.2 KiB
Julia
57 lines
2.2 KiB
Julia
# SPDX-FileCopyrightText: 2015-2026 Jukka Aho
|
||
# SPDX-License-Identifier: MIT
|
||
|
||
using Test
|
||
using JuliaFEM
|
||
using Tensors
|
||
|
||
@testset "update_material_cache! behavior branches" begin
|
||
mesh = create_test_mesh()
|
||
N = 8
|
||
NIP = 8
|
||
geometry_cache = JuliaFEM.create_geometry_cache(N, NIP)
|
||
elem_id = 1
|
||
Δt = 0.0
|
||
|
||
@testset "StatelessStrainDependent (NeoHookean)" begin
|
||
mat = NeoHookean(E_mod=210e9, nu=0.3)
|
||
kernel = ContinuumKernel(ContinuumFormulation{FullThreeD}(), mat)
|
||
element_cache = JuliaFEM.create_element_cache(mesh, kernel)
|
||
material_cache = JuliaFEM.create_material_cache(mat, NIP)
|
||
global_cache = JuliaFEM.create_global_material_cache(mat; n_ips=NIP, n_elems=1)
|
||
|
||
ug = [zero(Vec{3,Float64}) for _ in 1:nnodes_total(mesh)]
|
||
ug[2] = Vec((0.05, 0.0, 0.0))
|
||
|
||
JuliaFEM.update_geometry_cache!(geometry_cache, element_cache, elem_id, mesh)
|
||
JuliaFEM.update_element_cache!(element_cache, kernel, elem_id, mesh, ug)
|
||
JuliaFEM.update_material_cache!(
|
||
material_cache, geometry_cache, mat, element_cache, global_cache, elem_id, Δt,
|
||
)
|
||
|
||
σ = JuliaFEM.get_stress(material_cache, 1)
|
||
@test norm(σ) > 1.0
|
||
end
|
||
|
||
@testset "StatefulStrainDependent (PerfectPlasticity)" begin
|
||
mat = PerfectPlasticity(; E=210e9, ν=0.3, σ_y=200e6, H=0.0)
|
||
kernel = ContinuumKernel(ContinuumFormulation{FullThreeD}(), mat)
|
||
element_cache = JuliaFEM.create_element_cache(mesh, kernel)
|
||
material_cache = JuliaFEM.create_material_cache(mat, NIP)
|
||
global_cache = JuliaFEM.create_global_material_cache(mat; n_ips=NIP, n_elems=1)
|
||
|
||
ug = [zero(Vec{3,Float64}) for _ in 1:nnodes_total(mesh)]
|
||
ug[1] = Vec((0.0, 0.0, 0.0))
|
||
ug[2] = Vec((0.002, 0.0, 0.0))
|
||
|
||
JuliaFEM.update_geometry_cache!(geometry_cache, element_cache, elem_id, mesh)
|
||
JuliaFEM.update_element_cache!(element_cache, kernel, elem_id, mesh, ug)
|
||
JuliaFEM.update_material_cache!(
|
||
material_cache, geometry_cache, mat, element_cache, global_cache, elem_id, Δt,
|
||
)
|
||
|
||
σ = JuliaFEM.get_stress(material_cache, 1)
|
||
@test σ isa SymmetricTensor{2,3,Float64}
|
||
end
|
||
end
|