JuliaFEM Матрица жесткости для уравнения Пуассона приносит нулевую матрицу - PullRequest
0 голосов
/ 08 февраля 2019

Есть ли проблемы с пакетом JuliaFEM для Julia версии 1.0?Я пытался запустить пример из документации пакета JuliaFEM для уравнения эллиптического тепла.Я использовал теплопроводность k = 1, что делает задачу проблемой Пуассона.Результатом матрицы жесткости является нулевая матрица?

Я также пробовал ее на другой квадратной сетке, которую я создал, и той же проблеме нулевой матрицы ... Ожидаемым результатом должна быть матрица Лапласа с граничным условием Неймана

Код из примера:

using JuliaFEM
el1 = Element(Quad4, [1, 2, 3, 4])
X = Dict(1 => [0.0,0.0], 2 => [1.0,0.0], 3 => [1.0,1.0], 4 => [0.0,1.0])
update!(el1, "geometry", X)
update!(el1, "thermal conductivity", 1.0)
time = 0.0
el2 = Element(Tri3, [3, 2, 5])
X[5] = [2.0, 1.0]
elements = [el1, el2]
update!(elements, "geometry", X)
update!(elements, "thermal conductivity", 1.0)
problem = Problem(Heat, "test problem", 1)
add_elements!(problem, elements)
assemble!(problem, time)
Matrix(problem.assembly.K)

Ссылка, по которой я взял код: http://www.juliafem.org/JuliaFEM.jl/latest/packages/FEMBase/problems.html

...