@@ -68,60 +68,6 @@ function jacobian(physicalmodel::ElectroMechano, kine::NTuple{2,KinematicModel}
6868 jacobian (physicalmodel, ElectroMechano, kine, (u, φ), (du, dφ), (v, vφ), dΩ, Λ)
6969end
7070
71-
72- # =================
73- # Visco-elasticity
74- # =================
75-
76- # -----------------
77- # Stagered residual
78- # -----------------
79-
80- # function residual(physicalmodel::ViscoElectricModel, ::Type{Mechano}, kine::NTuple{2,KinematicModel}, (u, φ), v, dΩ, Λ, Δt, un, A)
81- # DΨ = physicalmodel(Λ, Δt=Δt)
82- # F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
83- # E = get_Kinematics(kine[2]; Λ=Λ)
84- # ∂Ψu = DΨ[2]
85- # ∫(∇(v)' ⊙ (∂Ψu ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)))dΩ
86- # end
87-
88- # function residual(physicalmodel::ViscoElectricModel, ::Type{Electro}, kine::NTuple{2,KinematicModel}, (u, φ), vφ, dΩ, Λ, Δt, un, A)
89- # DΨ = physicalmodel(Λ, Δt=Δt)
90- # F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
91- # E = get_Kinematics(kine[2]; Λ=Λ)
92- # ∂Ψφ = DΨ[3]
93- # -1.0*∫(∇(vφ) ⋅ (∂Ψφ ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)))dΩ
94- # end
95-
96- # # -----------------
97- # # Stagered jacobian
98- # # -----------------
99-
100- # function jacobian(physicalmodel::ViscoElectricModel, ::Type{Mechano}, kine::NTuple{2,KinematicModel}, (u, φ), du, v, dΩ, Λ, Δt, un, A)
101- # DΨ = physicalmodel(Λ, Δt=Δt)
102- # F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
103- # E = get_Kinematics(kine[2]; Λ=Λ)
104- # ∂Ψuu = DΨ[4]
105- # ∫(∇(v)' ⊙ ((∂Ψuu ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)) ⊙ (∇(du)')))dΩ
106- # end
107-
108- # function jacobian(physicalmodel::ViscoElectricModel, ::Type{Electro}, kine::NTuple{2,KinematicModel}, (u, φ), dφ, vφ, dΩ, Λ, Δt, un, A)
109- # DΨ = physicalmodel(Λ, Δt=Δt)
110- # F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
111- # E = get_Kinematics(kine[2]; Λ=Λ)
112- # ∂Ψφφ = DΨ[6]
113- # ∫(∇(vφ)' ⋅ ((∂Ψφφ ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)) ⋅ ∇(dφ)))dΩ
114- # end
115-
116- # function jacobian(physicalmodel::ViscoElectricModel, ::Type{ElectroMechano}, kine::NTuple{2,KinematicModel}, (u, φ), (du, dφ), (v, vφ), dΩ, Λ, Δt, un, A)
117- # DΨ = physicalmodel(Λ, Δt=Δt)
118- # F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
119- # E = get_Kinematics(kine[2]; Λ=Λ)
120- # ∂Ψφu = DΨ[5]
121- # -1.0*∫(∇(dφ) ⋅ ((∂Ψφu ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)) ⊙ (∇(v)')))dΩ -
122- # ∫(∇(vφ) ⋅ ((∂Ψφu ∘ (F∘∇(u)', F∘∇(un)', E∘∇(φ), A...)) ⊙ (∇(du)')))dΩ
123- # end
124-
12571# -------------------
12672# Monolithic strategy
12773# -------------------
0 commit comments