spring1d_elem.m

matlab实例

%***********************************************************************
%  function spring1d_elem
%
%  calculates the stiffness matrix of a 1D (co-linear) spring
%***********************************************************************
function [kelem, ke] = spring1d_elem(ElemData, Prop)

prop_i = ElemData.propID;

% NOTE: using propID as a property index.
% TO DO: check if Geom.Prop(prop_i).propID = prop_i

k = Prop(prop_i).properties(1);

kelem(1,1) =  k;
kelem(1,2) = -k;
kelem(2,1) = -k;
kelem(2,2) =  k;

ke = zeros(2,1);