Files
JuliaFEM.jl/test/test_api.jl
T

66 lines
1.6 KiB
Julia
Raw Normal View History

2015-11-30 19:11:00 +02:00
# This file is a part of JuliaFEM.
# License is MIT: see https://github.com/JuliaFEM/JuliaFEM.jl/blob/master/LICENSE.md
module APITests
using JuliaFEM
using JuliaFEM.Test
2015-11-30 21:50:09 +02:00
using JuliaFEM: Model, Node, Element, ElementSet,
Material, Quad4, Seg2
2015-11-30 19:11:00 +02:00
function test_basic()
# basic workflow, copied from test_solver.jl
model = Model("my calculation model")
# create nodes
2015-11-30 21:50:09 +02:00
model.nodes[1] = [0.0, 0.0]
model.nodes[2] = [1.0, 0.0]
model.nodes[3] = [1.0, 1.0]
model.nodes[4] = [0.0, 1.0]
2015-11-30 20:03:28 +02:00
2015-11-30 19:11:00 +02:00
# create elements
2015-11-30 21:50:09 +02:00
e1 = Element{Quad4}([1, 2, 3, 4])
e2 = Element{Seg2}([1, 2])
2015-11-30 19:11:00 +02:00
2015-11-30 20:03:28 +02:00
# element set
elset = ElementSet("body", e1, e2)
elset2 = ElementSet("boundary", e2)
2015-11-30 21:50:09 +02:00
# push!(model, elset, elset2)
2015-11-30 20:03:28 +02:00
2015-11-30 19:11:00 +02:00
# material properties
2015-11-30 21:50:09 +02:00
# material = Material("my material")
# material["temperature thermal conductivity"] = 6.0
# material["density"] = 36.0
# set_material!(model, material, "Es")
2015-11-30 19:11:00 +02:00
2015-11-30 20:03:28 +02:00
# Create problem
2015-11-30 21:50:09 +02:00
# field_problem = HeatProblem()
2015-11-30 19:11:00 +02:00
# boundary conditions
2015-11-30 21:50:09 +02:00
# bc = DirichletBC("body", "temperature" => 0.0)
2015-11-30 20:03:28 +02:00
# bc = DirichletBC("Es", 0.0)
2015-11-30 21:50:09 +02:00
# ne = NeumannBC("boundary", "temperature flux" => ((0.0 => 0.0),
# (1.0=>600.0)))
2015-11-30 20:03:28 +02:00
# LoadCase
2015-11-30 21:50:09 +02:00
# case = LoadCase("Heat problem")
# push!(case, field_problem)
# push!(case, bc, ne)
2015-11-30 20:03:28 +02:00
# Get solver
2015-11-30 21:50:09 +02:00
# solver = LinearSolver
# push!(case, solver)
2015-11-30 20:03:28 +02:00
# add case
2015-11-30 21:50:09 +02:00
# push!(model, case)
2015-11-30 20:03:28 +02:00
# Solve problem
2015-11-30 21:50:09 +02:00
# solve!(model, 1.0, "Heat problem")
# xi = [0.0, -1.0]
# T = el1("temperature", xi, 1.0)
2015-11-30 20:03:28 +02:00
# T = model("temperature", [0.5, 0.0], 1)
2015-11-30 19:11:00 +02:00
end
2015-11-30 21:50:09 +02:00
test_basic()
2015-11-30 19:11:00 +02:00
end