mirror of
https://github.com/JuliaFEM/JuliaFEM.jl.git
synced 2026-09-19 01:48:47 +00:00
Jukka Aho
55 lines
1.3 KiB
Julia
55 lines
1.3 KiB
Julia
using JuliaFEM.elasticity_solver
|
|
|
|
using FactCheck
|
|
using Logging
|
|
@Logging.configure(level=DEBUG)
|
|
|
|
|
|
facts("test solve elasticity increment") do
|
|
X = [0 0; 10 0; 10 1; 0 1]'
|
|
elmap = [1; 2; 3; 4]
|
|
nodalloads = [0 0; 0 0; 0 -2; 0 0]'
|
|
@debug("nodal loads:\n", nodalloads)
|
|
dirichletbc = [0 0; NaN NaN; NaN NaN; 0 0]'
|
|
|
|
E = 90
|
|
nu = 0.25
|
|
mu = E/(2*(1+nu))
|
|
la = E*nu/((1+nu)*(1-2*nu))
|
|
la = 2*la*mu/(la + 2*mu)
|
|
|
|
la = la*ones(1, 4)
|
|
mu = mu*ones(1, 4)
|
|
u = zeros(2, 4)
|
|
du = zeros(2, 4)
|
|
R = zeros(2,4)
|
|
Kt = zeros(8,8)
|
|
|
|
N(xi) = [
|
|
(1-xi[1])*(1-xi[2])/4
|
|
(1+xi[1])*(1-xi[2])/4
|
|
(1+xi[1])*(1+xi[2])/4
|
|
(1-xi[1])*(1+xi[2])/4
|
|
]
|
|
|
|
dNdξ(ξ) = [-(1-ξ[2])/4.0 -(1-ξ[1])/4.0
|
|
(1-ξ[2])/4.0 -(1+ξ[1])/4.0
|
|
(1+ξ[2])/4.0 (1+ξ[1])/4.0
|
|
-(1+ξ[2])/4.0 (1-ξ[1])/4.0]
|
|
|
|
ipoints = 1/sqrt(3)*[-1 -1; 1 -1; 1 1; -1 1]
|
|
iweights = [1 1 1 1]
|
|
|
|
for i=1:10
|
|
JuliaFEM.elasticity_solver.solve_elasticity_increment!(
|
|
X, u, du, R, Kt,elmap, nodalloads, dirichletbc, la, mu, N, dNdξ, ipoints, iweights)
|
|
@debug("increment:\n",du)
|
|
u += du
|
|
if norm(du) < 1.0e-9
|
|
break
|
|
end
|
|
end
|
|
@debug("solution\n",u)
|
|
@fact u[2, 3] => roughly(-2.222244754401764) # Tested against Elmer solution
|
|
end
|