stiffnessmatrixequationsolutionwithelasticrestraint.m

平面刚架（带弹簧支座）有限元程序

function [d,r]=StiffnessMatrixEquationSolutionWithElasticRestraint(n,ne,nr,s,ac,sr,crl)
s(n+1:n+ne,n+1:n+ne)=s(n+1:n+ne,n+1:n+ne)+diag(sr(n+1:n+ne));%把支座刚度系数累加到总节点刚度矩阵之中
d(1:n+ne+nr,1:1)=0;
r(1:n+ne+nr,1:1)=0;
d(1:n+ne)=inv(s(1:n+ne,1:n+ne))*ac(1:n+ne);%计算非约束位移
r(n+ne+1:n+ne+nr)=-ac(n+ne+1:n+ne+nr)+s(n+ne+1:n+ne+nr,1:n+ne)*d(1:n+ne);%计算支座反力
r(n+1:n+ne)=-diag(sr(n+1:n+ne))*d(n+1:n+ne);%计算弹性支座反力
dx=d(crl(1:n+ne+nr));d=dx;%构建总节点位移向量
dx=r(crl(1:n+ne+nr));r=dx;%够将总约束反力向量