<< plotsection FEMTruss

FEMTruss >> FEMTruss > rank_efficient_struct

rank_efficient_struct

Test the validity of a stiffness matrix

Calling Sequence

wrong_struct = rank_efficient_struct(K,p,e)

Parameters

K

the stiffness matrix of the structure (computed via a call to femtruss)

p

a matrix containing the coordinate of the points. It's a npoints x 2 matrix for a 2 dimensions structure and a npoints x 3 for a 3 dimensions structure.

e

a matrix containing the list of fixed degrees of freedom.

wrong_struct

%T if the structure produces a stiffness matrix which is not of full rank (this can mean that some bars are missing).

Description

Test the validity of a stiffness matrix

Examples

[t,p,e,A,E,rho,F] = build_fem_test('bridge2d');

[U,P,R,K,M]= femtruss(build_fem_test, %F, 'bridge2d');

res = rank_efficient_struct(K,p,e);

if (res) then
  printf('rank efficient structure\n');
else
  printf('not a rank efficient structure\n');
end

See Also

Authors


<< plotsection FEMTruss