Custom export of von Mises stress

Author: kaipartmann

src/auxiliary/io.jl

@@ -161,7 +161,7 @@ end
 
 function export_fields!(vtk, chunk, fields::Vector{Symbol}, t)
     for field in fields
-        point_data = export_field(Val(field), chunk, t)
+        point_data = export_field(Val(field), chunk.mat, chunk.system, chunk.storage, t)
         vtk[string(field), VTKPointData()] = point_data
     end
     return nothing
@@ -174,8 +174,8 @@ function export_fields!(vtk, chunk, fields_spec::Dict{Symbol,Vector{Symbol}}, t)
 end
 
 # this function can be specialized for each field, even custom export fields can be written!
-function export_field(::Val{field}, chunk, t) where {field}
-    return get_loc_point_data(chunk.storage, chunk.system, field)
+function export_field(::Val{field}, mat, system, storage, t) where {field}
+    return get_loc_point_data(storage, system, field)
 end
 
 @inline function get_loc_position(chunk::AbstractBodyChunk)

src/physics/correspondence.jl

@@ -182,15 +182,16 @@ end
     @pointfield damage::Vector{Float64}
     bond_active::Vector{Bool}
     @pointfield n_active_bonds::Vector{Int}
-    @pointfield stress::Matrix{Float64}
+    @pointfield cauchy_stress::Matrix{Float64}
     @pointfield von_mises_stress::Vector{Float64}
 end
 
 function init_field(::CMaterial, ::AbstractTimeSolver, system::BondSystem, ::Val{:b_int})
     return zeros(3, get_n_points(system))
 end
 
-function init_field(::CMaterial, ::AbstractTimeSolver, system::BondSystem, ::Val{:stress})
+function init_field(::CMaterial, ::AbstractTimeSolver, system::BondSystem,
+                    ::Val{:cauchy_stress})
     return zeros(9, get_n_loc_points(system))
 end
 
@@ -239,8 +240,7 @@ function calc_first_piola_kirchhoff!(storage::CStorage, mat::CMaterial,
     P = first_piola_kirchhoff(mat.constitutive_model, storage, params, F)
     PKinv = P * Kinv
     σ = cauchy_stress(P, F)
-    update_tensor!(storage.stress, i, σ)
-    storage.von_mises_stress[i] = von_mises_stress(σ)
+    update_tensor!(storage.cauchy_stress, i, σ)
     return PKinv
 end
 
@@ -315,3 +315,13 @@ function too_much_damage!(storage::AbstractStorage, system::AbstractSystem,
     end
     return false
 end
+
+# calculate the von Mises stress from the Cauchy stress tensor just when exporting
+function export_field(::Val{:von_mises_stress}, ::CMaterial, system::BondSystem,
+                      storage::AbstractStorage, t)
+    for i in each_point_idx(system)
+        σ = get_tensor(storage.cauchy_stress, i)
+        storage.von_mises_stress[i] = von_mises_stress(σ)
+    end
+    return storage.von_mises_stress
+end

src/physics/correspondence_rotated.jl

@@ -148,7 +148,7 @@ end
     @pointfield damage::Vector{Float64}
     bond_active::Vector{Bool}
     @pointfield n_active_bonds::Vector{Int}
-    @pointfield stress::Matrix{Float64}
+    @pointfield unrotated_stress::Matrix{Float64}
     @pointfield von_mises_stress::Vector{Float64}
     @pointfield left_stretch::Matrix{Float64}
     @pointfield rotation::Matrix{Float64}
@@ -164,7 +164,8 @@ function init_field(::CRMaterial, ::AbstractTimeSolver, system::BondSystem, ::Va
     return zeros(3, get_n_points(system))
 end
 
-function init_field(::CRMaterial, ::AbstractTimeSolver, system::BondSystem, ::Val{:stress})
+function init_field(::CRMaterial, ::AbstractTimeSolver, system::BondSystem,
+                    ::Val{:unrotated_stress})
     return zeros(9, get_n_loc_points(system))
 end
 
@@ -231,11 +232,10 @@ function calc_first_piola_kirchhoff!(storage::CRStorage, mat::CRMaterial,
     Δε = D * Δt
     Δθ = tr(Δε)
     Δεᵈᵉᵛ = Δε - Δθ / 3 * I
-    σ = get_tensor(storage.stress, i)
+    σ = get_tensor(storage.unrotated_stress, i)
     σₙ₊₁ = σ + 2 * params.G * Δεᵈᵉᵛ + params.K * Δθ * I
-    update_tensor!(storage.stress, i, σₙ₊₁)
+    update_tensor!(storage.unrotated_stress, i, σₙ₊₁)
     T = rotate_stress(storage, σₙ₊₁, i)
-    storage.von_mises_stress[i] = von_mises_stress(T)
     P = first_piola_kirchhoff(T, F)
     PKinv = P * Kinv
     return PKinv
@@ -249,3 +249,14 @@ function calc_zem_stiffness_tensor!(storage::CRStorage, system::BondSystem, mat:
     C_1 = calc_rotated_zem_stiffness_tensor!(zem_stiffness_rotated, params.C, Kinv, R)
     return C_1
 end
+
+# calculate the von Mises stress from the Cauchy stress tensor just when exporting
+function export_field(::Val{:von_mises_stress}, ::CRMaterial, system::BondSystem,
+                      storage::AbstractStorage, t)
+    for i in each_point_idx(system)
+        σ = get_tensor(storage.unrotated_stress, i)
+        T = rotate_stress(storage::AbstractStorage, σ, i)
+        storage.von_mises_stress[i] = von_mises_stress(T)
+    end
+    return storage.von_mises_stress
+end