Constitutive models and ForwardDiff: Repalce Float64 by Real

src/PhysicalModels/ElectricalModels.jl

@@ -4,8 +4,8 @@
 # ===================
 
 struct IdealDielectric <: Electro
-  ε::Float64
-  function IdealDielectric(; ε::Float64)
+  ε::Real
+  function IdealDielectric(; ε::Real)
     new(ε)
   end
 end

src/PhysicalModels/ElectroMechanicalModels.jl

@@ -119,7 +119,7 @@ end
 
 struct FlexoElectroModel{EM<:ElectroMechano} <: FlexoElectro{EM}
   electromechano::EM
-  κ::Float64
+  κ::Real
 
   function FlexoElectroModel(electro::E, mechano::M; κ=1.0) where {E<:Electro,M<:Mechano}
     physmodel = ElectroMechModel(electro, mechano)

src/PhysicalModels/MagneticModels.jl

@@ -4,17 +4,20 @@
 # ===================
 
 
-struct Magnetic <: Magneto
-  μ::Float64
-  αr::Ref{Float64}
-  χe::Float64
+struct Magnetic{T<:Real} <: Magneto
+  μ::Real
+  αr::Base.RefValue{T}
+  χe::Real
 
-function Magnetic(; μ0::Float64, αr::Float64, χe::Float64=0.0)
-  new(μ0, Ref(αr), χe)
+  function Magnetic(; μ0::Real, αr::Real, χe::Real=0.0)
+    T = typeof(αr)
+    new{T}(μ0, Ref{T}(αr), χe)
+  end
 end
+
 function (obj::Magnetic)(Λ::Float64=1.0)
-  μ, αr, χe = obj.μ, obj.αr, obj.χe
-  ℋᵣ(N) = αr[] * Λ * N
+  (; μ, χe) = obj
+  ℋᵣ(N) = obj.αr * Λ * N
   # Energy #
   Ψmm(ℋ₀, N) = (-μ / 2.0) * ((ℋ₀ + ℋᵣ(N)) ⋅ (ℋ₀ + ℋᵣ(N))) * (1 + χe)
   ∂Ψmm_∂φ(ℋ₀, N) = (-μ) * (ℋ₀ + ℋᵣ(N)) * (1 + χe)
@@ -23,18 +26,15 @@ function (obj::Magnetic)(Λ::Float64=1.0)
 end
 
 
-end
-
-
 struct IdealMagnetic <: Magneto
-  μ::Float64
-  χe::Float64
-  function IdealMagnetic(; μ0::Float64, χe::Float64=0.0)
+  μ::Real
+  χe::Real
+  function IdealMagnetic(; μ0::Real, χe::Real=0.0)
     new(μ0, χe)
   end
   function (obj::IdealMagnetic)(Λ::Float64=1.0)
 
-    μ, χe = obj.μ, obj.χe
+    (; μ, χe) = obj
     J(F) = det(F)
     H(F) = J(F) * inv(F)'
 
@@ -72,15 +72,15 @@ end
 
 
 struct IdealMagnetic2D <: Magneto
-  μ::Float64
-  χe::Float64
-  function IdealMagnetic2D(; μ0::Float64, χe::Float64=0.0)
+  μ::Real
+  χe::Real
+  function IdealMagnetic2D(; μ0::Real, χe::Real=0.0)
     new(μ0, χe)
   end
 
   function (obj::IdealMagnetic2D)(Λ::Float64=1.0)
 
-    μ, χe = obj.μ, obj.χe
+    (; μ, χe) = obj
     J(F) = det(F)
     H(F) = J(F) * inv(F)'
 
@@ -113,20 +113,21 @@ end
 
 
 struct HardMagnetic <: Magneto
-  μ::Float64
-  αr::Float64
-  χe::Float64
-  χr::Float64
-  χt::Float64
-  βmok::Float64
-  βcoup::Float64
-  function HardMagnetic(; μ0::Float64, αr::Float64, χe::Float64=0.0, χr::Float64=8.0, χt::Union{Float64,Nothing}=nothing, βmok::Float64=0.0, βcoup::Float64=0.0)
-    χt_val = isnothing(χt) ? χe : χt
-    new(μ0, αr, χe, χr, χt_val, βmok, βcoup)
+  μ::Real
+  αr::Real
+  χe::Real
+  χr::Real
+  χt::Real
+  βmok::Real
+  βcoup::Real
+
+  function HardMagnetic(; μ0::Real, αr::Real, χe::Real=0.0, χr::Real=8.0, χt::Union{Real,Nothing}=nothing, βmok::Real=0.0, βcoup::Real=0.0)
+    χt = isnothing(χt) ? χe : χt
+    new(μ0, αr, χe, χr, χt, βmok, βcoup)
   end
 
   function (obj::HardMagnetic)(Λ::Float64=1.0)
-    μ, χe = obj.μ, obj.χe
+    (; μ, χe) = obj
     J(F) = det(F)
     H(F) = J(F) * inv(F)'
 
@@ -162,22 +163,24 @@ struct HardMagnetic <: Magneto
 end
 
 
-struct HardMagnetic2D <: Magneto
-  μ::Float64
-  αr::Ref{Float64}
-  χe::Float64
-  χr::Float64
-  χt::Float64
-  βmok::Float64
-  βcoup::Float64
-  function HardMagnetic2D(; μ0::Float64, αr::Float64, χe::Float64=0.0, χr::Float64=8.0, χt::Union{Float64,Nothing}=nothing, βmok::Float64=0.0, βcoup::Float64=0.0)
-    χt_val = isnothing(χt) ? χe : χt
-    new(μ0, Ref(αr), χe, χr, χt_val, βmok, βcoup)
+struct HardMagnetic2D{T<:Real} <: Magneto
+  μ::Real
+  αr::Base.RefValue{T}
+  χe::Real
+  χr::Real
+  χt::Real
+  βmok::Real
+  βcoup::Real
+
+  function HardMagnetic2D(; μ0::Real, αr::Real, χe::Real=0.0, χr::Real=8.0, χt::Union{Real,Nothing}=nothing, βmok::Real=0.0, βcoup::Real=0.0)
+    χt = isnothing(χt) ? χe : χt
+    T = typeof(αr)
+    new{T}(μ0, Ref{T}(αr), χe, χr, χt, βmok, βcoup)
   end
 
   function (obj::HardMagnetic2D)(Λ::Float64=1.0)
 
-    μ, χe = obj.μ, obj.χe
+    (; μ, χe) = obj
     J(F) = det(F)
     H(F) = J(F) * inv(F)'
 
@@ -208,3 +211,21 @@ struct HardMagnetic2D <: Magneto
     return (Ψmm, ∂Ψmm_∂u, ∂Ψmm_∂φ, ∂Ψmm_∂uu, ∂Ψmm_∂φu, ∂Ψmm_∂φφ)
   end
 end
+
+
+function Base.getproperty(obj::Union{Magnetic,HardMagnetic2D}, prop::Symbol)
+  if prop === :αr
+    return getfield(obj, :αr)[]
+  else
+    return getfield(obj, prop)
+  end
+end
+
+
+function Base.setproperty!(obj::Union{Magnetic,HardMagnetic2D}, prop::Symbol, val)
+  if prop === :αr
+    return getfield(obj, :αr)[] = val
+  else
+    return setfield!(obj, prop, val)
+  end
+end

src/PhysicalModels/MagnetoMechanicalModels.jl

@@ -68,12 +68,10 @@ function _getCoupling(mag::HardMagnetic, ::Mechano, Λ::Float64=1.0)
   # Hard magnetics in ultra-soft magnetorheological elastomers enhance fracture toughness and 
   # delay crack propagation, Journal of the Mechanics and Physics of Solids,
 
-
-  μ, αr, χr, χt, βcoup, βmok = mag.μ, mag.αr, mag.χr, mag.χt, mag.βcoup, mag.βmok
+  (; μ, χr, χt, βcoup, βmok) = mag
   J(F) = det(F)
   H(F) = det(F) * inv(F)'
 
-
   #-------------------------------------------------------------------------------------
   # FIRST TERM
   #-------------------------------------------------------------------------------------
@@ -85,7 +83,7 @@ function _getCoupling(mag::HardMagnetic, ::Mechano, Λ::Float64=1.0)
   Hℋ₀(F, ℋ₀) = H(F) * ℋ₀
   Hℋ₀Hℋ₀(F, ℋ₀) = Hℋ₀(F, ℋ₀) ⋅ Hℋ₀(F, ℋ₀)
 
-  ℋᵣ(N) = αr * Λ* N
+  ℋᵣ(N) = mag.αr * Λ* N
   Fℋᵣ(F, N) = F * ℋᵣ(N)
   Ψcoup(F, N) = (μ * J(F)) * (Fℋᵣ(F, N) ⋅ Fℋᵣ(F, N) - ℋᵣ(N) ⋅ ℋᵣ(N))
   ∂Ψcoup_∂F(F, N) = 2 * (μ * J(F)) * (Fℋᵣ(F, N) ⊗ ℋᵣ(N))
@@ -150,7 +148,7 @@ function _getCoupling(mag::HardMagnetic2D, ::Mechano, Λ::Float64=1.0)
   # Hard magnetics in ultra-soft magnetorheological elastomers enhance fracture toughness and 
   # delay crack propagation, Journal of the Mechanics and Physics of Solids,
 
-  μ, αr, χr, χt, βcoup, βmok = mag.μ, mag.αr, mag.χr, mag.χt, mag.βcoup, mag.βmok
+  (; μ, χr, χt, βcoup, βmok) = mag
   J(F) = det(F)
   H(F) = det(F) * inv(F)'
 
@@ -164,7 +162,7 @@ function _getCoupling(mag::HardMagnetic2D, ::Mechano, Λ::Float64=1.0)
   #-------------------------------------------------------------------------------------
   Hℋ₀(F, ℋ₀) = H(F) * ℋ₀
   Hℋ₀Hℋ₀(F, ℋ₀) = Hℋ₀(F, ℋ₀) ⋅ Hℋ₀(F, ℋ₀)
-  ℋᵣ(N) = αr[] * Λ* N
+  ℋᵣ(N) = mag.αr * Λ* N
   Fℋᵣ(F, N) = F * ℋᵣ(N)
   Ψcoup(F, N) = (μ * J(F)) * (Fℋᵣ(F, N) ⋅ Fℋᵣ(F, N) - ℋᵣ(N) ⋅ ℋᵣ(N))
   ∂Ψcoup_∂F(F, N) = 2 * (μ * J(F)) * (Fℋᵣ(F, N) ⊗ ℋᵣ(N))