mirror of
https://github.com/JuliaFEM/JuliaFEM.jl.git
synced 2026-09-19 09:54:55 +00:00
Jukka Aho
38 lines
1.0 KiB
Julia
38 lines
1.0 KiB
Julia
using JuliaFEM
|
|
using Base.Test
|
|
|
|
function test_one_element_solution()
|
|
model = new_model()
|
|
|
|
# 1. create mesh structure
|
|
# add nodes to model
|
|
nodes = Dict(1 => [0.0, 0.0], 2 => [10.0, 0.0], 3 => [10.0, 1.0], 4 => [0.0, 1.0])
|
|
add_nodes(model, nodes)
|
|
|
|
# 2. add elements to model
|
|
const QUAD4 = 0x4
|
|
element = Dict("element_type" => QUAD4, "node_ids" => [1, 2, 3, 4])
|
|
elements = Dict()
|
|
elements[1] = element
|
|
add_elements(model, elements)
|
|
|
|
# 3. add boundary conditions
|
|
# add dirichlet boundary condition u=0
|
|
boundaries = Dict("SUPPORT" => [1, 4])
|
|
add_dirichlet_boundary_condition(model, boundaries)
|
|
|
|
# add nodal load to node 3
|
|
nodal_loads = Dict(3 => [-2, 0])
|
|
add_nodal_loads(mode, nodal_loads)
|
|
|
|
# 4. model is defined, solve it
|
|
JuliaFEM.solvers.solve_elasticity!(model; parallel=false)
|
|
|
|
# 5. extract results
|
|
#write_xdmf(m, "results", ["displacement"])
|
|
disp = get_field(model, "displacement")
|
|
@assert_eq disp[2, 4] == -2.22224475
|
|
end
|
|
|
|
test_one_element_solution()
|