Commit 7e382500 authored by Andreas Zilian's avatar Andreas Zilian
Browse files

Revert to just consider elastic (material) stiffness in computing total stiffness.

parent 18b557c3
......@@ -42,8 +42,8 @@ classdef FeProblem < handle
self.postprocessElements(z)
end
%% assemble system K matrix
function A = assembleSystemStiffnessMatrix(self)
A = self.assembleSystemStiffnessMatrixElastic() + self.assembleSystemStiffnessMatrixGeometric();
function A = assembleSystemStiffnessMatrix(self)
A = self.assembleSystemStiffnessMatrixElastic() ;%+ self.assembleSystemStiffnessMatrixGeometric();
end
%% assemble system K matrix: elastic
function A = assembleSystemStiffnessMatrixElastic(self)
......@@ -59,7 +59,7 @@ classdef FeProblem < handle
% get element index vector
e_idx = elem.getDofSysIndices();
% get element matrix
e_mat = elem.stiffnessElastic();
e_mat = elem.stiffness();
% determine number of entries
e_size = length(e_idx)^2;
% add element contribution to system matrix
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment